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Editorial 


POLISH MARITIME RESEARCH is a scientific journal of worldwide circulation. The journal appears as 
a quarterly four times a year. The first issue of it was published in September 1994. Its main aim is to 
present original, innovative scientific ideas and Research & Development achievements in the field of : 


Engineering, Computing & Technology, Mechanical Engineering, 


which could find applications in the broad domain of maritime economy. Hence there are published papers 
which concern methods of the designing, manufacturing and operating processes of such technical objects 
and devices as : ships, port equipment, ocean engineering units, underwater vehicles and equipment as well 


as harbour facilities, with accounting for marine environment protection. 
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Editor's message 


This is a special issue of Polish Maritime Research quarterly, devoted to results of the 
research on design, manufacture and operation of a new generation of ships for Baltic 
Sea shipping. The research was realized in the frame of the EU-supported project aimed 
at creating multi-mode shipping routes in that region, and satisfying the conceptual as- 
sumptions of Baltic Sea status as a Sensitive Sea Area. For that reason great importance 
has been attached to ecological problems in this project. We hope that this initiative of 
the Editors and the Prinicpal Coordinator of the project will be met with kind acceptance. 
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ABSTRACT 


The paper outlines the origin of the E.2772! project from the EUREKA 
group of European projects. Justification of undertaking the subject is gi- 
ven as well as the structure of institutional performers constituting the pro- 
јес! executive consortium. Characterised is the final result in the form of 
preliminary designs of four modern ecological Baltic short shipping ves- 
sels: container carrier, product tanker, ro-ro ship and river-sea type ship. 
Comprehensive bibliography contains a list of all the project reports, which 


е 


may be made available to the interested parties after a contact with the project coordinator. The paper 
opens a series of topical papers selected from the different problems dealt with in the project. 


Keywords: EUREKA project Baltecologicalship, project execution, preliminary designs of four ships 


INTRODUCTION 


The structural changes and enormous acceleration of the 
world economy cause a necessity of developing new genera- 
tions of surface transport, including sea and inland waterway 
shipping. This applies to transportation of the investment and 
consumer goods as well as to passenger traffic. In the highly 
developed areas, strongly urbanized and inhabited by weal- 
thier societies, one can hardly imagine further development of 
the traditional transport of goods through the road or even rail- 
way network, to meet the growing demand. With cautious fo- 
recasts of the increase of freight transport in the EU territory 
by approx. 38% by 2010, concentrated in the main land trans- 
port corridors and around large urban areas, activation of the 
sea and inland water transport incorporated into multimodal 
transport networks becomes a problem of vital importance. 

The emphasis put on the development of sea transport has 
been for the last several years a clear policy of the European 
Union, where surface transport is to be divided into long di- 
stance transport between established and intensively develo- 
ped central reloading nodes and short distance feeder trans- 
port. In the case of sea transport, the latter is defined as short 
shipping. This problem is also characteristic ofthe Baltic, which 
has become practically an EU Mare Nostrum. 

The amount of cargoes transported annually in the Baltic 
area reaches 450 million tons. Approximately 5500 thousand 
trips are performed served by the 500 Baltic ports. In recent 
years, the annual increase of the volume of transported goods 
has been as big as 30%. It is estimated that some 5000 ships 
are operated in the Baltic, 40% of them are older than 20 years 
and will have to be replaced by new generation ships in the 
nearest future. 

The Baltic short shipping development trends include ma- 
inly container carriers (300% increase of the volume of trans- 
port in 10 years), passenger ferries and product tankers. The 
actual structure and size of that new generation fleet will de- 
pend mainly on the long term development forecasts of the 
volume and type of cargo and the passenger tourist traffic. 


Regardless of location of the Baltic multimodal transport 
corridors, ecological requirements will have to be met connec- 
ted with new conception ofthe Baltic Sea status as PARTICU- 
LARLY SENSITIVE SEA AREA, to be implemented by 2005. 
That involves a necessity of creating a fleet of specialised new 
generation Baltic feeder service ships fulfilling strict ecologi- 
cal criteria. Designs of such ships, based not only on the veri- 
fied design practice but also on the state-of-the-art knowledge, 
should be on offer in the shipyards, and shipping companies 
aspiring to participation in the European transport system sho- 
uld provide for operating such ships. 


PRO-ECOLOGICAL REQUIREMENTS 


Regardless of the type of ship, the following 
general pro-ecological requirements have to be met : 


У Significant reduction of: 


c» all kinds of emissions and pollution connected with ship 
operation 

effects of possible sea averages and disasters 

negative effects of the solid waste utilization processes 
effects of vibration and noise 

negative effects of ship construction and repair pro- 
cesses. 


0000 


2 Improved safety and reliability of ship operation (structure 
and systems). 

> Application of additional technologies, equipment and sys- 
tems for proper utilization of the operational waste pro- 
ducts. 


Sources of hazards to the sea environment 
from an operating ship are schematically presented in Fig.1. 


Meeting all the requirements of the new generation Baltic 
ships will certainly cause higher construction costs but it is 
estimated that they will be compensated for by e.g. lower insu- 
rance premiums, lower port charges and smaller hazards to the 
marine environment. 
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PROJECT EXECUTION 


In order to perform the task of establishing the Baltic eco- 
logical multimodal transport corridors with new generation 
ships, a European project E!2772 Baltecologicalship was set 
up in the EUREKA group and a consortium was formed of the 
following participants: 


Technical universities 


ж Gdansk University of Technology, Faculty of Oce- 
an Engineering and Ship Technology – chief exe- 
cutor and coordinator of the whole project; Faculty 
of Chemistry as responsible for preparations of se- 
lected shipbuilding sector companies for the ISO 
9000 and 14000 certification 

3* Szczecin University of Technology, Faculty of Ma- 
ritime Technology. 


Companies from the small 
and medium enterprise (SME) sector 


3* SINUS Limited Company, performing design work 
on the conceptions of 4 types of Baltic ships 

3* Toka Consulting Company, cooperating in the ISO 
reviews 

% Ship repair yard, which unfortunately withdrew 
from the project at an early stage. 


Foreign partners 


ж Boy&Partners-Loma Company, Sweden, co-perfor- 
ming the ISO reviews 

ж Björn Carlsson Ecoship Engineering AB, Malmö, 
Sweden, performing own design work on the pro- 
duct tanker conception. 


The consortium performed the work in two stages. 


At the definition stage, studies were undertaken on the ne- 
eded structure of the Baltic region sea transport as well as de- 
sign conceptions and requirements of the ships. 

At the design and testing stage, designs of the ships were 
carried out and, as they were being completed, laboratory tests 
and analyses of hull shapes were performed and also verifying 
calculations and analyses of machinery, equipment and system 
solutions were carried out. In effect, the designs were correc- 
ted as appropriate and submitted to final evaluation. Parallel 
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Fig. 1. Ecological hazards created by ship (acc. to DNV — Maritime Solutions — Clean Performance advertising materials) 


deliveries 
to-shore — — 


bilgewater 


work was carried out to recognize the production capabilities 
of selected Polish shipyards. An advanced modernization con- 
ception of a selected shipyard was developed to ensure effecti- 
ve realisation of the contract. 

The project participants, both the university teams and the 
small and medium enterprise (SME) sector companies, cre- 
ated an efficient and well functioning organism which has ge- 
nerated up-to-date designs forming a prospective development 
offer of the shipbuilding sector in its aspirations for the Euro- 
pean operating structural funds. 

At the same time, as a consequence of activities in the Swe- 
dish branch ofthe project, in cooperation with Polish specialist 
groups, 40 shipbuilding sector enterprises, including two me- 
dium size shipyards, have been prepared for beginning the ISO 
9000 and 14000 standard certification process. The Swedish 
conception of a product tanker has obtained approval of a Swe- 
dish owner and after many rounds of negotiations in different 
shipyards, final contract talks are in progress in two selected 
shipyards. 


THE SINE SERIES SHIPS 
The project comprised the following ship types : 


> SINE 202 universal container carrier 
> SINE 203 oil product tanker 

> SINE 204 ro-ro ship 

> SINE 205 river-sea ship. 


The designed ships had to fulfil the Det Norske Veritas 
(DNV) CLEAN class conditions, the first in the world com- 
plex set of classification requirements to guarantee ecological 
cleanness of a ship. A range of different solutions have been 
applied in the designed ships in order to meet the requirements. 

In general, the solutions comprise : 


œ In connection with the increased safety requirements of ship 
in operation and in case of disastrous collision, the follo- 
wing solutions have been applied : 


** diesel-electric propulsion system (better use of fuel, mini- 
mum engine starting and stoppage time, significant limita- 
tion of water pollution by oil from the propulsion system) 

** pod propulsion and steering system (effectiveness in- 
creased by at least 10%, better manoeuvrability, less vi- 
bration and noise, elimination of rudder, minimum ship 
backing time, elimination of the cooling system) 


Environment friendly ships for Baltic area 


ve two separate power plants (breakdown safety of pro- 
pulsion and hull) 

эх double shell in the fuel tank and hold area (elimination 
of spills). 


In connection with restrictions on the emissions of gases to 
atmosphere and of the operational and technological media 
to water, the following solutions have been applied : 


* use of low sulphur content fuels 

Ус use of exhaust-gas catalyzers (limitation of NO, and 
SOx) 

ус design constraints on oil spills and cooling medium 
spills, elimination of freon 

* minimalization of ballast water exchanges in order to 

avoid the hazards from propagation of microflora and 

microfauna 

use of exclusively natural media in the fire fighting sys- 

tems 

* limitation of the oil content in bilge waters to a 15 ppm 
level with an automatic STOP device switching the sys- 
tem off when the value is exceeded 

* use of ecologically safe internal and external ship pa- 
ints 

** use of special ship sewage treatment procedures. 


x 


In connection with solid waste treatment aboard ship, to be 
transferred to shore for utilisation, waste sorting machines 
have been applied. 


Main parameters of the designed ships are listed in the Table. 


Fig. 2 presents general arrangement plan and 3D visualisa- 


tion of the ro-ro ship stern showing elements of the systems 
necessary to meet the CLEAN class requirements. Figs. 3, 4 
and 5 present the product tanker, container carrier and river- 
-sea ship, respectively. The ship documentation has been pre- 
pared to the level of preliminary design of the feeder service 
means of transport. 


У 


У 


У 
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Table 


Ro-ro ship Product tanker Container carrier 
| Ly[m] — | 14775 | — 1320 — | — 1320 — 
| — B[m] | 02480  Á | 2250 | 2250 — 
| Htomaindeck[m] | — todeck419.60 | — — 12.80 | tomaindeck 11.20.— 
|___Тагјт] | | |— 60 — | 80 |70 O 
|___ T[n] | | | 650 |). 870 | X 855 O 
|  DWT |740 | 12800 | 8550 


Loa [M 


L 

B [m 

H to main deck [m 
Taes. 
T [m 
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18.05 
8000 6000 10000 


m 
m 


19500 5980 11200 
Generator [kW/V] | 690 | 1120/660 1120/660 


pod pod 
2x2200 2x4350 


pod 
2x6500 


Propulsion 
[kW] 


SUMMARY 


The Eureka group European projects are an element of sys- 
tem support to the development of complex innovative pro- 
ducts of small and medium enterprises obtaining their know- 
ledge from the resources of scientific and research institu- 
tions. They create one of the mechanisms of organizing 
a European research and development network. 

The project execution consortium, sponsored by the Com- 
mittee for Scientific Research on the Polish side and the 
VIANOVA and NUTEK government agencies on the Swe- 
dish side, brought together a research potential of 2 techni- 
cal universities (3 faculties), many small and medium en- 
terprise sector companies from the Polish side and 2 com- 
panies of that sector from the Swedish side and effectively 
combined science, knowledge, theory and design practice 
in realisation of a complex innovative task. 

The cooperation proved effective. Four ship designs have 
been created on the Polish side and three designs (product 
tanker) on the Swedish side. Forty Polish shipbuilding sec- 
tor enterprises have been prepared for the ISO 9000 and 
14000 standard series certification and also talks are in pro- 
gress about an order of 4 product tankers for a Swedish 
owner. These are immediate effects of the cooperation. 
The obtained and processed knowledge — a basis of the 
design work — has been gathered together in 159 reports 
and other source materials, which will also be used for pre- 
paring a monograph publication (see bibliography). 

The product tanker design, selected with the use of synthe- 
sis of knowledge from earlier European projects, has been 
introduced as a subject of continued studies and analyses to 
the next Eureka group projects (ASPIS). 

Selected interesting topics of ће analyses (in hydromecha- 
nics, design, structural design, ship systems, economy, pro- 
duction engineering) performed in the project will be pu- 
blished later on in this journal. 


Type of ship 
SINE 205 
River-sea ship 
89.45 
87.47 
11.40 
5.45 
4.40/2.80 


2900 
12.0/14 
2500 
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2x1120/690 


2х700 


Prop. © [m] | 
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7200 m? — trailer 
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С шаа 
> 


Holds – 54 
Deck - 36 
385 т? 


950 ТЕО 
ог 475 ЕЕ4 


+1A1 Tanker for Oil 
Products, ESP, EO 

ETC, ICE-!C, RPS, 
ICS, W1, CCO, LCS 
(SID) 


+1A1 Container 
Carrier, EO, W1, 
S.C., RP 


*1A1 General Cargo 
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Fig. 3. Product tanker (SINE 203) 


POLISH MARITIME RESEARCH, Special issue, 2004 9 


EUREKA 


ENVIRONMENT FRIENDLY SHIPS FOR BALTIC AREA 
EUREKA PROJECT - E!2772 BALTECOLOGICALSHIP 


i 
— пааша шш шт 


T пп 
D Elec ЕРГЕ 
ја Ров АП L | | 
EE - 7 A = ФФ 
a == 7 Ж \ ) 
ELI Е чы е/ Q 
A-DECK WHEELHOUSE FORECASTLE DECK 
B-DECK D-DECK E-DECK PLE 
c = : 
/ MD А 
X Ж LAN 
"P at 
[=>] 
MAINDECK MAINDECK 
| | | 
/ / " 
TN AN 7 TN v 
7 а c» SS 
/ NX || N N i у N 
) = CUR = = = ў . 
и | EÁZS FJA 12 77 PRINCIPAL DIMENSIONS : 
7 7 ; N 
/ NP N V NIA b / \ Ж 
Length overall Loa = ~ 13810 m 
Length between perpend. Lbp = 132,00m 
Breadth moulded В = 2250m 
Depth H = 1120m 
Scantling draught TS = 8,55 т 
DECK 9000 Design draught Td = 760m 
3 | | Service speed У = 18,5 kn 
Ш ЦЕ | sj 5 = Se Deadweight at 7,60 m draught ab.8550 DWT 
ГЕСЕР = tH fs d 2 - = = s e 
Ан] 
ж Ш m) бі CREW 14 *1 
С | JE 
CLASS: 
DNV + 1A1 Container Carrier, EO, W1, SC, RP 
CLEAN DESIGN 
FLOOR: DOUBLE BOTTOM 
WH, | ІҢ ^| 
Sy d М р SINE202 
fh [XP 
а | h в в, 


10 


Fig. 4. Container carrier (SINE 202) 
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Fig. 5. River-sea ship (SINE 205) 
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RIVER - SEA VESSEL 


PRINCIPAL DIMENSIONS : 


LENGTH OVERALL 89,45 т 
LENGTH B.PERPENDICULARS 87,47 т 
BREADTH MOULDED 11,40 т 
DEPTH 545m 
SEA DRAUGHT 440m 
RIVER DRAUGHT 2,80m 
DEADWEIGHT 2900 DWT 
° CREW 8 
SPEED 12,0 knots/14,0 km 
TEU 20' IN HOLDS 54 
ON DECK 36 
HOLD CAPACITY 73850 cub m 
CLASS: 


DNV + 1A1 GENERAL CARGO CARRIER 
HC, EO, ICE-1C, RP 


CLEAN DESIGN 


SINE205 
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ABSTRACT 
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Paper presents results of studies on a parametric method of predicting ship building costs — 
— useful in the preliminary design. Conception and theoretical basis of the method are 
^ | presented, devised are also approximation formulae for estimating the building costs of the 
ship hull, ship equipment and the power plant with propulsion system. Factors of importan- 
ce for the ship building costs are identified and a computational algorithm formulated. The 
useful character of the method is illustrated by examples of building cost predictions for 
four different ship types designed in the Eureka project E!2772, i.e.: SINE 202 universal 


container carrier, SINE 203 oil product tanker, SINE 204 ro-ro ship and SINE 205 river-sea ship. 


Keywords: ships design theory, ship building costs 


INTRODUCTION 


At the initial stages of designing a ship, its parameters are 
determined by the ship design theory methods allowing to find 
design solutions fulfilling the owner design requirements as 
well as the valid international conventions, e.g. freeboard rules 
or the ship damage stability requirements. In the optimum de- 
sign methodology it is also necessary to formulate the desig- 
ned ship evaluation criterion. 

The criterion measure should be a ship property — of a me- 
asurable economic effectiveness character — important for the 
owner's operations on the market and containing, in a direct or 
indirect form, ship investment costs, predicted operating costs, 
evaluation of the ship transporting and money earning capabi- 
lity as well as the investment and acting capital costs. The ship 
evaluation criterion measure must be expressible by a set of 
ship technical parameters x determined in the ship design pro- 
cess - decision variables of the design optimization mathema- 
tical model. Optimization models in the methods of ship de- 
sign theory may be particularly useful in the work of ship de- 
sign offices at the time of stronger competition on the world 
shipbuilding and shipping markets. 

This paper contains results of research studies on a para- 
metric method of preliminary prediction of the ship building 
costs, based on a small set of its main design parameters. The 
developed method of ship building cost prediction (a top-down 
method) may be used as an evaluation criterion in the prelimi- 
nary ship design optimization models. The presented method 
is a contribution to the development of ship design theory as 
well as a tool to be applied in practical ship design. 

The problem of methodology of ship building cost predic- 
tion at the preliminary design stage has been a subject of many 
publications, both those already classic, like the works of Ben- 
ford [1, 2, 3], Sójka [4], Buxton [5], Fetchko [6], Dart [7] and 
Fisher [8] and the more recent items, e.g. [9, 10, 11, 12, 13]. 
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AIM AND SCOPE OF THE WORK 


Collecting data on real costs of building ships is a difficult 
task — in particular in relation to the ships currently under con- 
struction — as the data are treated as confidential business in- 
formation of shipyards and owners. For evaluation of the eco- 
nomic value of a designed ship important are, apart from the 
production costs, also the predicted future operating costs, avai- 
lability of cargoes, freight rates etc. in an appropriate time pe- 
riod, e.g. 15 to 30 years of the expected ship operation. Acqui- 
ring a reliable prediction tool is a significant difficulty in the 
research work on the development of preliminary ship design 
optimization methods. 

The ship building and operating cost data are scarce in the 
literature, usually they are related to different time periods and 
are presented in a form of diagrams or tables of little use to the 
optimization method computational algorithms. 

The presented research was initially aimed at results consis- 
ting in the development of ship design theory through devi- 
sing an original parametric method of ship building cost eva- 
luation at the preliminary design stage, a method of a mathe- 
matical structure appropriate for the computer optimization 
algorithms. 

A useful effect of the research was demonstrated by using 
the method to perform the building and operating cost predic- 
tions for ships designed in the EUREKA project E!2772. 


REQUIREMENTS OF THE METHOD 


In the early phases of preliminary ship design the buil- 
ding cost may be predicted by means of parametric relations 
corresponding with the level of the identified technical para- 
meters. At that stage no data are available to predict the cost 
from the equipment and material pricelists or from the ma- 
ker's proposals. The parameters in question are main dimen- 
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sions, component weight estimates, speed, propulsion power 
etc.. The method is based on the available data related to ship 
parameters and building costs as well as information on the 
cost structure. 


In view of the non-homogeneity of data, as they come 
from different sources and are differentiated as regards : 


» territory — they pertain to shipyards in different countries, 
operating in different economic systems 

> monetary system — they are expressed in different currencies 

> time- they pertain to different periods in the past. 


The costs have been reduced to comparable values – diffe- 
rent currencies were converted to US dollars at the rates valid 
in the respective periods and then updated to the current value 
using an average US inflation rate (in the years 1990-2003) of 
3% (according to the Bank of America indices). 


It has been assumed that the total 
ship building cost including : 


material cost 
labour cost 
other shipyard costs 


consists of : 


hull construction 
ship equipment 
power plant and 
propulsion system. 


It has been also assumed that costs of each of these groups 
are related to the weight of a respective group. Approximating 
functions were determined from the collected data to express 
the unit cost in US dollars per ton of a group weight. The po- 
wer plant horse-power was converted to an equivalent power 
plant weight. An advantage of such approach is balancing the 
designed ship displacement with the sum of component wei- 
ghts, e.g. by the Normand method, already in the preliminary 
stages of the design process. 


Depending on a functional type of the ship, 
the hull cost is considered for two cases, as in [11]: 


+ 


* single-deck hulls, e.g. tankers or bulk carriers 
* multi-deck hulls, with well developed internal volume sub- 
division, e.g. passenger ships, car carriers or ro-ro ships. 


* 


D 


The ship equipment cost 15 considered in three classes, dif- 
fering in respect of the equipment quality and of saturation of 
the ship with equipment. The equipment unit prices and the 
assembly labour cost may differ considerably. The equipment 
class variant depends on the owner requirements as well as on : 


* functional type of the ship (e.g. classic general cargo ships 
usually have low equipment class whereas passenger ships 
have high equipment class, etc.) 

* necessity of installing special equipment, e.g. refrigerating 
or air conditioning plant, according to owner's requirements. 


The power plant cost predicting formulae 
are related to diesel engine plants, where : 


Ж lower range ofthe main engine weight (power) corresponds 
to medium-speed engine power plants 

Ж upper range of the main engine weight (power) corresponds 
to low-speed engine power plants. 


Cost shares of the individual ship building stages in the 
total cost, according to an empirical estimation quoted in [14], 
are the following: 


Building stage Cost Impact on total 
5 stag of the stage | building costs 


The estimation shows that the design stage, having itself 
approx. 10% share in the total ship building costs, determines 
85% of those costs. Expenses on the design quality — proper 
choice of the ship main parameters, production technology, 
structural materials, equipment types etc. — have a significant 
impact both on the shipyard's and owner's economic effects. 


UNIT COSTS 


Preliminary studies have shown that analytical relations 
between a unit cost qj of a ship technological group and its 
weight mj may be approximated with sufficient accuracy by 
means of power functions containing four constant structural 
coefficients ej; determined by the least squares method. An 
identical structure of the unit cost and weight binding formula 
has been assumed for all the technological groups : 

9; = сој t0,;:m; c5; mj?! (1) 
The unit cost determined by this formula comprises: 


** material costs 
** labour costs 
* other shipyard costs. 


Structure of the formula was determined from tests perfor- 
med with different test functions, using the following selection 
criteria: 


$ good approximation quality in the considered range of de- 
sign parameters 

$ easy determination of numerical values of the formula as 
well as derivative relations, which is an important advanta- 
ge when manual calculations have to be performed, e.g. 
with a calculator, and also when an optimum design ma- 
thematical model is formulated (a problem of accurate de- 
termination of the Jacobian or Hessian determinant). 


The structural coefficients of approximation formulae, gi- 
ven in Table 1, have been determined by means of the regres- 


sion analysis. 
Table 1. Values of structural coefficients 
of the ship building cost predicting formulae 


0.053 


-1.175 |12957.0] 0.0518 


standard 
equipment 


25230.2| 2.928 |-161.66 | 0.6343 
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The preliminary ship building cost predictions, confronted 
with actual market prices (as published by the Centromor fore- 
ign trade agency in the „Budownictwo Okretowe" monthly), 
were evidently on the high side. That may be attributed to the 
fact that prices for some types of ships have decreased during 
the last decade not only in a relative sense (allowing for infla- 
tion) but also in an absolute sense, as an effect of sharp compe- 
tition on the shipbuilding market and dumping prices quoted 
by the Far East shipyards. 

Correction of the predictions was performed by calibration 
of the formulae with the use of the mentioned published actual 
prices. New corrected structural constants of the approxima- 
tion formulae are given in Table 2. 


Table 2. Corrected structural coefficients 
of the ship building cost predicting formulae 


. Unit 
group 


Single 1994.663 | 0.015549 | -154.0222 0.2471932 
deck hull 
Multi 
3241.7105|0.0194715 -5926707 n 1637139 
deck hull 


Low 
3| standard 6528.5325|9.6026706 -12.56288 | 0.979517 
equipment 


Average 
4| standard Me -2599.825 | -1.395667 | 9146.288 |0.0213938 
equipment 


High 
5| standard Qn |9749.0427| 14.66748 | -16.71265 |0.9963722 
eqipment 
TOIAN 16720.374|0.7839685| -221.3641 0.510682 510682 
plant 


The unit construction cost qs [$/t] of a single-deck hull 
of weight mp € (1000 + 30000) [t] is the following : 


qs = Co, +С "р t C5, ang Q) 
and that of a multi-deck hull 
of weight mp € (1000 + 30000) [t] : 


= C39 
qm = 20, + C12 Шу *C55:m, (8) 


The unit cost де [$/t] of ship equipment of weight 
me € (100 + 3000) [t], depending on the equipment standard : 


low standard equipment : 
= 63,3 
qı = Соз + C1,3°Me t C55: Me (4) 
average standard equipment : 
= 63,4 
Яа = Сод FC, 4°M, t C5 4m, (5) 
high standard equipment : 
= 3,5 
Gn = Co,5 +С15°Ш * C55: m, (6) 


The unit cost qy [$/t] of the power plant 
of weight mp € (100 + 2500) [t] : 


= ©з 6 
qp = Сов t+ Cy 6 m, tC56 m, (7) 


The formula validity ranges correspond to the statistical 
samples used for their approximation. As the unit cost func- 
tions are smooth, a moderate extrapolation is admissible. 


THE SHIP BUILDING 
COST PREDICTION ALGORITHM 


In the preliminary ship design stages, the component we- 
ight groups are usually determined by methods based on the 
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Fig. 1. Ship hull construction costs 
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parent ship data, e.g. by the Normand or Bubnov method or 
from the empirical and statistical relations, e.g. those given by 
Watson [11]. In both cases the component weights are expres- 
sed by main parameters X = (x,,x;...., x,) ofthe designed ship 
in functional analytical (modular) formulae : 


m, E m, (x) (8) 


so selected that they approximate the actual relations in the 
best possible way. The modular formulae determine the we- 
ight increases due to the changes of main design parameters of 
the ship, described by the X = (x4, x5,..., x,) vector, in relation 
to the parent ship parameters X, = (X01; X02- X94), Which in 
the domain of linear relations, for the mj weight group, is de- 
scribed by the expression : 


m;(x) = m;(x9) + Am;(x) = m;(xo) + 


m,-m,(x);m,-m,(x); 


n Qm. (X) У) 
j 
+ by 3 | (x; = Хо) 
м 9% y x 
The individual group costs are described 
by an expression of the following structure : 
0,00 = 4;(%) m) (10) 


The total ship building cost Q; , depending on the functio- 
nal type, saturation with equipment and equipment quality, may 
be expressed as a sum of respective components : 


QE) =Q, mEt 1, 00+%60 (0D 
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Table 3. Main design parameters of ships and predicted building costs in US dollars 


Building cost prediction 
simulation parameters for the 
baltecologicalships project 
ships 
m 


ight 
Depth to maindeck 


Scantling draught 


Speed at design draught 


| 


Power at 1800 rpm 


Number of decks : 
1- single-deck ; 2 — multi-deck 
Equipment standard : 
1— low, 2 — average, 3 – high 


Equipment cost 


[5] 


BUILDING COST PREDICTIONS 
FOR SHIPS DESIGNED 
IN THE EUREKA PROJECT E!2772 


The table above contains design parameters of the analysed 
ships, used for predictions of building costs. The total building 
cost of a ship consists of the hull construction cost, equipment 
cost and the cost of power plant with propulsion system. 


SUMMARY AND CONCLUSIONS 


| | | | | jeee 


О А parametric (top-down) method has been presented of pre- 
dicting the building costs of ships, to be used at the preli- 
minary design stages. The method may be used to determi- 
ne a ship evaluation criterion measure in the design metho- 
dology based on formal mathematical optimization models. 
The optimum design models may prove useful in the work 
of ship design offices, particularly at the time of strong com- 
petition on the world shipbuilding and shipping markets. 


O The presented original parametric method of the estima- 
tion of ship building costs, with a mathematical structure 
applicable to computer optimization algorithms, is a con- 
tribution to the development of ship design theory. 


O It is also a tool for practical ship design and was used for 
performing predictions of the building and operating costs 
of ships designed in the EUREKA project E!2772. 


О It is expected that empirical verification of the prediction 
results will be the cost estimations prepared in the respecti- 
ve shipyards constructing the analysed ships. 


SINE 202 
Container 
carrier 


8 550 14 300 
15 600 19 922 
138.70 138.10 


SINE 205 
River-sea 
ship 


2 950 

3 938 

89.45 

87.47 

11.40 
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3 640 

4 000 
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Paper presents results of studies on a parametric method of predicting ship operating 
costs - useful in the preliminary ship design. Conception and theoretical basis of the me- 
thod are presented, identified are also factors of significant importance for the ship opera- 
ting costs, taking into account changes in the value of money. Approximation formulae for 
estimating the operating cost components have been developed as well as a computational 
algorithm based on a minimum Required Freight Rate (RFR). The useful character of the 
method is illustrated by examples of operating cost predictions for four different ship types 


designed in the Eureka project E!2772, i.e.: SINE 202 universal container carrier, SINE 203 oil product 
tanker, SINE 204 ro-ro ship and SINE 205 river-sea ship. 
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INTRODUCTION 


In the preliminary ship design methodology formulated as 
an optimization problem, the aim is to determine such design 
solution which extremises a selected objective function (cho- 
ice criterion) dependent on the sought ship parameters. The 
criterion measure may be a measurable ship property impor- 
tant for the economic effects of owner's business operations. 
An important element of a preliminary ship design methodolo- 
gy are the ship operating cost prediction methods with costs 
expressed as relations dependent on the design decision varia- 
bles and the design mathematical model parameters. 

Atthe time of formulation of ship design requirements, in- 
formation on the future ship operation market conditions is very 
limited, uncertain and difficult to evaluate and any attempt to 
use it to perform expensive and time consuming studies aimed 
at optimizing the design requirements appears inefficient. There- 
fore, it may be justified to use simplified methods based on 
a small set of easy-to-get data of an index character in order to 
obtain approximate estimations with simple algorithms. That 
type of approach inspired the work on the presented method of 
predicting the ship operating costs by balancing the discoun- 
ted costs and incomes. 

The main design requirements of cargo ships, apart from 
the ship function, are usually: deadweight capacity Pn or car- 
go capacity Pc, service speed v, holds volume Ve and autono- 
mous steaming range R. Results of studies on the methods of 
predicting ship operating costs in early stages of design, per- 
formed on different levels of accuracy of the analysed models, 
were frequently published, e.g. [1-11]. Taking into account 
a greater number of model parameters leads to problem formu- 
lations which may be solved only with complex computational 
systems with the use ofthe optimization algorithms, e.g. those 
described in [12—14]. 
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PROBLEM FORMULATION 
AND ASSUMPTIONS 


A problem is under consideration where a ship of given 
deadweight capacity Pn, to be operating on a route of given 
length R, calls during trip at s ports for loading or unloading 
and it should achieve an assumed rate of return e during m 
years of ship operation. The required freight rate RFR 15 sou- 
ght that would be a minimum rate in given technical and eco- 
nomic conditions of ship operation but would cover all the ship 
operating costs taken into account in the model. 

Cargo supply in ports is generally of a stochastic character, 
particularly in the tramp shipping. A simplifying assumption 
is adopted in the method, allowing to use a deterministic mo- 
del in the considered line shipping case. The stochastic charac- 
ter of cargo supply is provided for in an indirect way by intro- 
ducing a coefficient €, expressing an average use of ship cargo 
capacity, which with sufficiently long ship operation periods 
allows to approximate an actual mean cargo supply. An advan- 
tage of this approach is easy way of determining the coeffi- 
cient, e.g. by analysis of the log-book records of ships opera- 
ting on a given shipping line. 

The Required Freight Rate (RFR) is a rate that the owner 
should obtain in order to ensure the assumed rate of return e 
(with the borne ship investment and operating costs) in the ship 
operating period of m years. 

The choice of RFR as a measure of covering the ship ope- 
rating costs is based on reasoning that with an established level 
of actual freight rates on a given line, the best profitability will 
achieve a ship with the lowest required freight rate. If the futu- 
re actual freight rates will be higher than the determined mini- 
mum freight rate RFR then the real rate of return e* will be 
higher than the assumed rate e, and if the opposite occurs then 
investment will not bring the assumed rate of return e. Accor- 
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ding to Schneekluth [9], in real operating conditions a mini- 
mum freight rate ship offers the greatest probability of achie- 
ving the required rate of return on shipping investment. 


INVESTMENT COSTS 


We assume, as in [11], that the total ship building cost J 
consists of the propulsion system cost Js, dependent only on 
the ship speed, and of other investment costs Jp, dependent on 
the ship size expressed by the deadweight capacity: 


J(Ne, Pn) = Js(Ne) + Jp(Pn) (1) 


According to [3], the Js and Jp costs are functions depen- 
dent mainly on the installed horse power Ne and deadweight 
capacity Pn and the functions are increasing slower than linear 
functions : 


Ја(Ме) = Cs- Ne? 
Jp(Pn)= Cp. Pn?” 


In this study, the Cs and Cp constants in expressions (2) 
were determined by the method described in [15]. 


SHIP ANNUAL CARGO CAPACITY (ACC) 


Q) 


The Annual Cargo Capacity (ACC) of a ship operating, on 
average, Z days in a year, is determined as proportional to the 
ship cargo capacity and the number of trips in a year, with ave- 
raging capacity usage coefficient € < 1: 


ACC=n-e-A-Pn (3) 
where A is a deadweight usage coefficient : 
Pn -Pz 
А---- (4) 
Рп 


where Pz means the weight of stores. 


The number of trips in a year n, with calling 
at s ports and the route length R, is expressed by : 
к= Z 7 Zw 
T B.Tm P-R 
where T means duration of one trip 
and Tm is the actual sailing time in one trip. 


G) 


The B > 1 coefficient is a correction taking into account the 
Tj times of reloading operations and roadstead waiting time in 
one trip : 


The value of B coefficient depends on the conditions on a given 
shipping line, number of ports, cargo handling efficiency etc. 


SHIP ANNUAL OPERATING COSTS (AOC) 


Taking the ship operating cost analysis [1] and [9] as 
a basis, it has been assumed that the ship Annual Operating 
Cost (AOC) depends mainly on the fuel cost. In order to sim- 
plify the model, it was assumed that the average annual cost 
of lubricating oil and power plant repairs (dependent on the 
engine power) would be taken into account as a correction 
coefficient > 1; then the annual operating costs are expres- 
sed as : 


AOC = y: n: Tm: Cj: бј: Ne (7) 


Cj means unit fuel price [$/t] 
and Gj means unit fuel consumption [t/(kW-h)]. 


Therefore, the annual ship operation cost 1s : 
2/3 


AoC = реп пасја 22200) А 
1 


У 2/3 2 ш 
„= Кеп-К·Рп ^. v 
Ca 


where Ке means a respective product of parameters. 


DISCOUNTED AVERAGE 
ANNUAL COST (AAC) 


With the required annual investment rate of return e and 
assumed m years of ship operation, the discounted Average 
Annual Cost (AAC) – operating cost and investment cost — is 
expressed as follows : 


AAC AOC 


CRF(e, m) CRF(e, m) 
where CRF means Capital Recovery Factor. With the invest- 
ment rate of return e and m years of ship operation, the CRF is 
determined from the formula : 


(9) 


e(14 e)" 
(l+e)™—1 


If a predicted rate of inflation i and the owner's income tax 
(on the profit calculated as a difference between the freight in- 
comes and operating costs) rate t are to be taken into account in 
the model, then the Capital Recovery Factor formula will be : 


(e+i)(l+e +i)" 


[(й+е+)”—1]@-9 


In such case, CRF* should be inserted 
in all the formulae containing the CRF expression. 


CRF(e,m) = (10) 


CRF" (e, m,i,t)- (11) 


The discounted average annual cost is : 


AAC = J- CRF(e, m) + AOC 


MINIMUM REQUIRED 
FREIGHT RATE (RFR) 


(12) 


Ifthe AAC is covered by continuous incomes from the car- 
go transport freight rate during the year then the rate of return 
e, as assumed in the CRF, is achieved with that freight rate. 


The required freight rate may be calculated from 
the required freight income to cover the AAC 
and the Annual Cargo Capacity (ACC) : 

RFR = AAC 
ACC 

After respective substitutions and transformations, an expli- 

cit relation is given between the required freight rate and the 

constants and variable parameters of the model : 


КЕК = 


(13) 


Pn? 


B-CRF(e,m) E Cs 
T S pv 1#®. Ca? Pn?’ "Yn 
u-Cj-Gj ui 


1253.13.15. Са 
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As all the factors and components of this equation have 
positive values, their impact on the minimum required freight 
rate may be evaluated in a simple way. 


EXAMPLES OF THE RESULTS 
OF SHIP OPERATING COST PREDICTIONS 


The following table contains technical and economic para- 
meters of the ships designed in the Eureka project E!2772, 
which were used to perform the ship operating cost predictions 
with the method described here. 


EVALUATION 
OF RESULTS AND CONCLUSIONS 


An impact of some model parameters on the predicted ope- 
rating cost is intuitively obvious. However, the obtained analy- 
tical relations allow to analyse the impact determined by the 
index exponents of parameters in the model. The following 
observations and conclusions may be drawn from the analysis : 


» increased RFR means higher ship operating costs 
» increased required rate of return e — expressed by increased 
CRF - implies higher RFR 


Table. Design parameters and operating cost predictions of the project ships, in US dollars 


RESULTS OF SHIP 
OPERATING COST Symbol] Unit 
PREDICTION SIMULATIONS 


SINE 202 
950 TEU 
8 550 DWT 
Container 
carrier 


SINE 203 
14 300 DWT 
Product tanker 


SINE 204 
7 400 DWT 
Ro-ro ship 


SINE 205 
2950 DWT 
River-sea 
ship 


Assumed ship deadweight capacity 8 550 14 300 7 400 2 950 
Price of a reference ship 12708330 | 12170670 | 23063402 3 544 370 
Price of a reference power plant 5 228 790 4 330 630 7 666 479 1 447 350 


Power of a reference ship engine 
Displacement of a reference ship | D | | 
Speed of a reference ship 

[t] 


Deadweight capacity Pn 
of a reference ship 
Cargo capacity [-] 
efficiency coefficient = 
Deadweight efficiency coefficient 


10 000 19 520 


3 640 
3 938 


.85 


Trip time to sailing time ratio 


| Trptmetosslingtimeraio | b | H | 11 | 
Cost (fuel+repairs+lub. oil) [-] 14 
Cost (fuel ` 
sa| x | x [| ж [| ж 


rift consumption |i [їй 
пасата [е 


Number of trips in a year 


Investment cost 


Discounted annual cargo 
handling cost (10 $/t 


Annual Cargo Capacity 89 763 189 353 
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7 


0 
0 


0.18 


[as | «|» | » 


3 544 373 


104 986 


50 223 


Parametric method of prediction of the ship operating costs 


» technical progress, characterised by better ship resistance- 
-propulsion efficiency and expressed by greater Ca coeffi- 
cient, causes decrease of RER 

> impact of the ship deadweight capacity Pn on RFR is am- 
biguous. For specific conditions of a task, the final effect 
depends on relations of the remaining model parameters 
and constants 

» the increase of efficiency factors in the model causes re- 
duction of the required minimum freight rate, 1.e. reduces 
the operating costs 

> increased ship hull and equipment prices cause increase of 
RFR 

> increased fuel price Cj implies increased RFR 

» technical progress in the engine design, leading to reduc- 
tion of the unit fuel consumption Gj, decreases RFR. 


The presented ship operating cost prediction method ana- 
lyses the cost as dependent on a number of factors. The studies 
presented in this paper were performed with average values of 
the factors, as given in the Table summarizing the operating 
cost simulation calculations. The obtained results should there- 
fore be read together with the respective parameter values, 
which reflect the current and future technical and economic 
relations. 


NOMENCLATURE 


Ca - Admiralty coefficient 

Cp - other cost coefficient 

Cs - installed power plant cost 

e* - real rate of return 

Gj - unit fuel consumption 

Jp - other investment cost 

Pe - cargo capacity 

Pz - weight of stores 

s - number of ports for reloading in one trip 
T - duration time of one trip 

Т, - time of reloading operation in one trip 
Tm- duration of one way voyage 

Ve - volume of holds 

В - waiting time correction coefficient 

n - deadweight/displacement ratio 


À - deadweight usage coefficient 
u - correcting factor for operating costs 
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Methodology of assessing the production 
capacity of selected shipyards 


Ryszard Pyszko, M.Sc., Eng. 


ABSTRACT 


The paper deals with selected subjects іп the European EUREKA ,,Baltecologicalship " 
E!2772 project. Within the project, preliminary designs of four feeder (short shipping) ships 


have been developed based on the market trend research, design knowledge and experi- 
mental investigations. The following tasks were to be performed for one of the ships: defi- 
ning optimum building conditions (timescale, production technologies, ecological aspect), 
selecting shipyards meeting the minimum production capacity requirement. For that pur- 
pose an index method was developed (based on the design and production indices). Then 


an index analysis was performed of the production capacities of selected shipyards (and also different 
configurations of cooperation between them) and results were presented in a form of the Gantt diagrams. 
Final conclusions were drawn. 


Keywords : structural and production process coefficients, index method, labour demand, worstand - hours (Stg) 


PROBLEM DEFINITION 


An increasing competition for contracts to build new cargo 
ships requires continuous watching the shipbuilding market as 
well as improving production methods and capabilities. At pre- 
sent, Polish shipyards have an order book large enough to ensure 
survival of the industry but they are not free of the shortage 
of capital [2]. This is a temporary situation and it makes ship- 
yards seek solutions necessary to win the competition. An exam- 
ple of such action is implementing the conception of the so 
called fractal factories [1]. 

In view of a complex and multistage process of building 
a ship from conception to delivery, certain stages may be di- 
stinguished where solutions should be sought to improve eco- 
nomic competitiveness and to meet the ecological production 
requirements. 

In the current economic conditions in the country and capi- 
tal weakness of the large production enterprises, one can har- 
dly expect any significant investments in the shipyard techni- 
cal infrastructure in the nearest future. Therefore, solutions 
should be looked for in the production engineering and organi- 
zation. With so defined strategy of competitiveness, the follo- 
wing aims may be formulated : 

In the technological conditions of Polish shipyards, 
a broad collaboration between shipyards and the coopera- 
ting enterprises, particularly from the small and medium 
enterprise (SME) group, should be developed in order to 
use the production and innovative potential in an optimum 
way. Such cooperation should be based on mutual econo- 
mic benefits and long-term planning as well as on the ship- 
building market development trends. 


Such collaboration depends 
on fulfilling the following conditions : 


1. The manufacturers have a well prepared (small number of 
modifications) and easily accessible (e.g. in an electronic 
form) documentation for customers. The present level of 
computer software and hardware allows to prepare com- 
plete ship construction documentation. [3, 4, 5]. 

2. Ship structure production quality standards for the project 
in question are implemented at the leading shipyard and at 
subcontractor plants, together with agreements between the 
quality control units. 
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3. Specialist teams are shifted between the cooperating enter- 
prises to perform specific tasks, or production of some struc- 
tural sections is outsourced and then transported to the lea- 
ding shipyard for assembly in dock. 

4. An efficient and effective transport operates between the 
enterprises (between the Gdansk-Gdynia area and the Szcze- 
cin area shipyards the waterway transport will be most ap- 
propriate). 

5. Technologies not meeting the ecological production requi- 
rements are eliminated. 

6. А clear and unified ship construction process is developed 
(together with all the quality control instruments) and im- 
plemented in the project executing enterprises. 

7. Aunified algorithm of assessment of the subcontractor pro- 
duction capacity is used to subdivide the tasks in a rational 
way. 


The above listed conditions comprise a broad range of dif- 
ferent problems. The problems of optimum production pro- 
cess planning and subdivision of technical tasks among the 
subcontractors were chosen for further analysis. 

With those conditions, a process algorithm was proposed 
in order to perform the main task. In view of the broad range 
of conditions, items 1, 2, 3, 4, 5 were assumed fulfilled and 
the emphasis was placed on the production process condi- 
tions 6,7. 

The analysis was performed on the SINE-203 product tan- 
ker hull [6], Fig.1. A series of simplifying assumptions related 
to the structure and production process were adopted. The more 
significant ones include : 


€ the subject of analysis is steel hull ofthe ship, without out- 
fitting 

€ the hull structure was simplified, with panel and stiffener 
elements only (the plate thicknesses and diversity of stiffe- 
ners were maintained) and with subdivision into flat and 
bent elements 

€ real arrangement of external shell seams was replaced with 
a virtual arrangement of the same total length (without lo- 
cating seams on the shell). 


With the above listed assumptions, a material balance 
and weld quantification were prepared. 
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Main dimensions : 


Length overall : Loa = 138.100 m 


Length b. perpendiculars : Lpp = 132.000 m 
Width : B -22.50m 
Depth : H =12.80m 
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Speed : V, = 14.0 knots 
Hull weight : Мк -4250t 
Draught : T -8.70m 
Underdeck volume (at T) : V =19375 т) 


#152+200 


Bow section 


Fig. 1. Subdivision of the SINE 203 ship into structural areas 


UNIFICATION 
OF THE TECHNOLOGICAL PROCESS 
OF SHIP HULL CONSTRUCTION 


Unification (schematisation) of the technological process 
was prepared to make comparisons possible between techno- 
logical capabilities of different shipyards. The process sche- 
matic diagram is presented in Fig.2. The most important stages 
of a technological process are included. 


Production stage Brief description of work 


Preliminary machining Cutting of plates and profiles 


Bending of plates and profiles 


ovine of machining line 


Fabrication of frames, girders, 


Preliminary fabrication ötacketsete. 


Fabrication of 2D objects Other lines, e.g. bent profiles line 


T Fom sections line 
Machine group method 


Construction of 3D objects 


v 
Fabrication of 3D objects : É 
outside a slipway (or a dock) 


Construction of 3D objects 
on a slipway (or in dock) 
Fig. 2. Unified production process of the ship hull construction 


Assembly 


SELECTION OF SHIPYARDS 
FOR THE COOPERATION ANALYSIS 


Analysis was performed with reference shipyards A, B and 
C of different production capacities. A hull production cycle 
of six months was assumed. 

Shipyard A has much greater production capacity than that 
required by the SINE 203 oil product tanker. 

Shipyard B covers a variant of cooperation between Ship- 
yard A and a small production capacity shipyard. Cooperation 
includes the preliminary machining, preliminary fabrication 
and 2D object fabrication stages. 


Shipyard C covers a variant of cooperation between ship- 
yards of similar production capacities, where the preliminary 
machining, preliminary fabrication, 2D object fabrication and 
3D object fabrication is performed by one shipyard and assem- 
bly in dock is performed by the other shipyard. This is a pro- 
duction cycle often practiced by the SME sector shipyards. 


INDEX METHOD OF ASSESSING 
THE SHIPYARD PRODUCTION CAPACITY 


Traditional method 


Assessing the capability of constructing a hull in the exis- 
ting production conditions requires an individual approach, 
i.e. an optimum production process should be found for a gi- 
ven hull geometry in a given shipyard with its existing infra- 
structure. 


Index method 


The indices determine the structure and production process 
(coefficients) assigned to the characteristic operations of indi- 
vidual construction stages, which allows to express objective- 
ly (in workstand-hours - Stg) the labour demand of respective 
stages. By determining the indices of individual shipyards, we 
obtain a picture of production capacities, which helps to take 
cooperation decisions (Fig.3). 


The production process coefficients T; are determined 
in relation to the shipyard infrastructure and they describe : 


> productivity, expressed by [t/Stg], [m?/Stg], [m/Stg] 
> size and weight, expressed by dimensions - a, b, с [m] and 
structure weight [t]. 


The T; values may be determined in several ways from : 


> counting the basic production operation times - classic ap- 
proach [7] 

» anadopted system of integrated calculation coefficients [8] 

> experience, i.e. an average value of annual production ba- 
lance. 


The structural coefficients W; are determined in relation 
to the hull structure (quantity of material to be processed in the 
production stages in order to obtain a hull) and they are : 
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Hull d 
! A set of informations on: _ 
! - geometry 

| -Structure method 


- materials 1 


Determination 
of the design coefficients 


Hull 
Structural coefficients 
W; (1- indices) 


Index 
method 


Traditional 


Infrastructure 
‚ A set of informations on: | 

- layout of compartments 

í -equipment characteristics ! 
| - production process lines | 


Determination 


of the production process coefficients 


Infrastructure 
Production process coefficients 
Т; (1— indices) 


Fig. 3. Flow diagram of the methods of comparing known hull structural and material conditions with production capabilities of a given shipyard 


А productivity related, expressed by [t], [m°], [m], [pes] 
A size and weight related, expressed by dimensions - a, b, c 
[m] and structure weight [t]. 


Production process coefficients. Using the adopted uni- 
fied production process divided into five stages (Fig. 2), coef- 
ficients were assigned to the respective work scopes. Their va- 
lues were determined by evaluation of the production capaci- 
ties (from the integrated calculation indices [7]) for specific 
work scopes in each ofthe analysed shipyards. The coefficients 
are presented in Table 1. 

Structural coefficients. From the documentation and the 
prepared material and weld length balance, coefficients were 
assigned to the respective work scopes. The coefficients corres- 
pond with the quantities of material to be processed in the pro- 
duction stages. The coefficients are presented in Table 2. 


A total of 15 production process coefficients and correspon- 
ding structural coefficients were adopted. In assigning values 
to the structural coefficients, a different hull structure subdivi- 
sion was assumed for each shipyard. 


RESULT OF ANALYSIS 


By linking the production process coefficients and structu- 
ral coefficients, the operation times Rj, expressed in [Stg], were 
obtained for the respective workstands, in accordance with the 
work scope assumed (Table 3). 

Gantt diagrams were prepared from the tables for different 
variants. Examples of hull construction schedules are presen- 
ted here below for the proposed A, B, C shipyards. 

The schedules presented here for the selected shipyards are 
one of the feasible variants with the assumed coefficients. 


Table 1. Identification of production process coefficients 


Item 


(tmin) depending 


Field dimensions (a/b/c) and load lifting capacity 
on the structure subdivision 


Description / Unit / Productivit 


Preliminary machinin 


line 


Fabrication of 3D objects 


Construction of 3D objects 
EN Fillet weld length [m/Stg] related to the 2D stage outside the-slinway (òr dodi) 


Demand for the operating area and crane lifting 
capacity needed to fabricate the object 


The 1, 2, 12, 14 coefficients are of the productivity evaluation character 
The 13, 15 coefficients are of the size and weight character 
[Stg] - workstand-hours. Time needed to perform a given scope of work on the workstand. Number of workers is not taken into account 


Item 


Determination of the plate sheet area 


Table 2. Identification of structural coefficients 


Description / Unit / Productivit Quantity of material to be processed 


Preliminary machining 


Я : ; Pl h in th г 
1 ! in the structure [m°] (3x12) - assumed sheet size ше зпеегатеа шше SOURS 
Determination of the length of stiffener | | 
2. уу, ete bi EIS епо РВИ ЕУ Total length of stiffeners in the structure 
in the structure [m 
Fabrication of 3D objects 
Fillet weld length [m/man-hour] ; ; | А 
12. related to the 3D stage List of welds differentiated according to type 
Field dimensions (a/b/c) and load carrying capacity Por a given s ишеге Шиш, demand 
13. M for the operating area and carrying capacity 
(tmin) depending on the structure subdivision 
of the transport means 
The 1, 2, 12, 14 coefficients are of the productivity evaluation character 
The 13, 15 coefficients are of the size and weight character 
[Stg] - workstand-hours. Time needed to perform a given scope of work on the workstand. Number of workers is not taken into account 
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Table 3. List of reference coefficients for the SINE 203 ship 


Preliminary machining line 
УУ, [m] / T; [m"/Stg 


W- [m] / T; [m/Stg 
Fabrication of 3D objects 


Wi; [m] / Ti; [m/Stg 


Тіз- Constructed object (a/b/c), object weight [t] 


Wi; - Workstand dimensions (a/b/c), load са 


One working month has 178 hours, with 3-shift work it is 534 [h] 


К; is a coefficient with index i. If it is a productivity-type 
coefficient then it assumes an assigned value [Stg]. The size 
and weight coefficients are assigned logical values: 1 when the 
operation is possible or 0 when the operation is not possible. 
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Fig. 4. Three hull construction schedules 


In the analysed cases, the size and weight coefficients 

were assigned the 1 values. 

The Shipyard A schedule shows a correct hull construc- 
tion process with the assumed coefficients. 

The Shipyard B schedule presents an incorrect process. 
Too large percentage (scope) of the tasks were entrusted to 
a weaker partner, which caused lengthening of the process. 

The Shipyard C schedule presents a correct process of the 
cooperation of shipyards, none of which would be individual- 
ly able to build the hull in the planned time of 6 months. With 
such way of determining the labour demand, the question of 
calculating the values of [Stg] and man-hours in different ship- 
yards remains open. Negotiations between the parties may be 
expected. 


ECONOMIC ASPECT 
OF THE CONSTRUCTION 


With so determined construction schedules (Fig. 4), the hull 
construction cost analysis may be carried out. The simplifying 
assumptions are the following : 


ж the [Stg] cost is the same for each construction stage 
ж the [Stg] cost for different subcontractors is the same. 


The schedules may be presented in a form 
of the so called logistic curve (S curve) - Fig. 5. 


30 000 
(> Stg] —e—shipyard А 
25 000 === Shipyard В 
= Mm Shipyard C 
20 000 
15 000 
10 000 
5 000 


[Stg] 


5000 


0 1000 2000 3000 4000 
Fig. 5. The S curves for shipyard A, B, C 


Three S Curves, reflecting the work intensiveness, were 
obtained for three shipyards (the assumed construction condi- 
tions). 

Shipyard A - carries out the hull construction in a most 
productive way, i.e. very intensively and in a short time. 

Shipyard B - carries out the hull construction in a least 
productive way and in a long time. 

Shipyard C - carries out the hull construction in a slightly 
worse way than Shipyard A, but still within acceptable limits. 


If we assume that the whole sum needed for the construc- 
tion, e.g. 100 units, is borrowed with an e.g. 2096 per annum 
interest, then the total cost will be as shown in Fig.6. 
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~ m Total cost of hull construction [units] 
BTotal cost [Stg/unit] 
250 | 229.7 
200 Я) 
156.4 
150 4 —] 
1091 1179 meg 111.4 
100 
50 = 
0 
Shipyard A Shipyard B Shipyard C 


Fig. 6. Construction costs and [Stg] costs for the analysed shipyards 


If a [> Stg/total hull construction cost (units)] relation 15 
formed then the Stg cost will be obtained (Fig.6). The follo- 
wing results were obtained for the analysed cases: Shipyard A 
had the best and Shipyard B the worst result. Although Ship- 
yard C had only a little longer construction time than Ship- 
yard A, the increase of [Stg/unit] is evident. The reason of 
distinct differences in the analysis (Total cost [Stg/unit]) is si- 
multaneous cumulation of the longer construction time (Ship- 
yard B and C) and longer credit interest payments. The hull 
outfitting was not taken into account in the analysis but it may 
be expected with certainty that cost differences ([X Stg / È units], 
Fig.6) between shipyards will be significantly increased. 

With the use of parametric description of a greater number 
of ships (structural coefficients) and parametric description of 
shipyard (and subcontractor) infrastructure, the following stu- 
dies may be performed : 


** production program for the leading shipyard, which may 
be optimized in respect of the distribution of tasks between 
the shipyard itself and the subcontractors 


* methods of financing the individual tasks 


* determination ofa minimum profitable production treshold 
and maximum production capacity of the leading shipyard 


** from the above studies, determination of a long-term pro- 
duction schedule for the whole industry branch; such pre- 
paratory work will allow to survive the low demand period 
and to start production of new ships as soon as the demand 
for them returns. 


These studies may be performed by building a mathemati- 
cal model ofthe shipbuilding industry branch. This is an inter- 
disciplinary task for the shipbuilding, banking, environmental 
engineering etc. experts. 


SUMMARY 


> The index method gives an answer whether in the techni- 
cal conditions characteristic of a shipyard the construction 
is possible, how long time it will take (with the assumed 
parameters), what will be the construction process ,,effi- 
ciency". 

> The index-based analysis may be applied to any type of 
ship. The differences between ship hull types will consist 
in the values ofthe structural coefficients, particularly tho- 
se connected with the structural subdivision of hull, i.e. 
weights and sizes of sections and blocks. 

> Comparison of the indices characteristic of a hull structure 
with indices characteristic of a shipyard infrastructure al- 
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lows to assess the construction capability in a rational (ob- 
jective) way. 

> The accuracy (reliability) of the proposed approach depends 
on correct indication of important points in the production 
process and proper correlation with the indices (characteri- 
stics) of the constructed hull. 

2 The method allows to compare in a simplified way the pro- 
duction capabilities of several shipyards in relation to the 
construction of the same hull. That gives a technical basis 
for the choice of contractor(s). 


Advantages of the index method. 
with precise identification of the production 
capabilities of Shipyard A, B, C : 


3* the proposed index method allows to perform a variant pre- 
liminary assessment of the hull construction effectiveness 

ж an exact information is obtained, in the form of a number 
or a logical value, related to the index character 

** the index structure (method of calculation) reflects its sco- 
pe and role from the point of view of efficient execution of 
the process 

** the whole construction process may be covered with an ap- 
propriate number of indices 

ж by increasing the number of indices, the construction pro- 
cess may be more precisely modelled together with cost 
estimations of individual operations (if needed) 

ж by carrying out a similar index analysis, one can easily 
learn of the strong and weak points of another shipyard. 


Disadvantages of the index method are the following : 


* some operations of the construction process are left out 

* too big scope of work represented by one index may lower 
the priority of an important operation, which in turn may 
lead to a too low (or too high) index value 

* partial overlapping of the index ranges, which leads to too 
high end values 

* the indices must be calculated by the same method for all the 
compared shipyards, regardless of their production profiles. 


Change of the method of calculation of one index makes it 
necessary to recalculate the indices for other shipyards. 


Directions of development: some actions 
may be indicated to improve the proposed method : 


A appoint a group of potential partners in the region 
A determine the indices, by the proposed method, for the ap- 
pointed group of partners 
А analyse the so far achieved results in hull construction : 
option | - individual work of a shipyard, 
option 2 - in cooperation with other partners 
^ compare and summarize the results. 
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Service margin - solution of the problem 
or a problem waiting for solution ? 


Janusz Stasiak, D.Sc., Eng. 


ABSTRACT 


Subject of the paper is the problem of predicting sea-going ship resistance and/or horsepo- 
wer in real operating conditions. Attention is drawn first of all to the need and possibility of 
revising the very much outdated but still used "service margin" method, which consists in 
adding an arbitrary percentage margin to the value of resistance or horsepower, relatively 
precisely determined for the calm water conditions. A negative impact in this respect is 
shown of the generally used delivery-acceptance procedures, where particular importance 
is attached to the ship propulsion tests on the "measured mile". The need of revising the 


"service margins" is a consequence of an obvious need for most efficient ships from the technical as well as 
economic point of view. Secondly, it is a "must" of permanent improving the ship design quality - the adequ- 
acy and accuracy of design methods. The work presents the "wave service margin" coefficient models. It is 
assumed that they may contribute to the necessary rationalization of the procedure of real ship resistance 
and/or horsepower determination. The work is based on the results of resistance tests of a series of ships 
designed within the Baltecologicalship project. The tests were carried out in the Chair of Ship Theory and 
Design of the Faculty of Ocean Engineering and Ship Technology, Gdansk University of Technology. 


Keywords : sea-going ship, design, seagoing qualities, resistance, main propulsion horsepower, service margins 


INTRODUCTION 


The role of ship hydromechanics in ship design 


There is no doubt that the hydromechanical properties of 
a sea-going ship are a fundamental quality (a set of qualities) 
of the craft. They determine its physical substance and, to 
a large extent, its functional effectiveness. It seems therefore 
fully justified to think that ship design, aimed at creating the 
best possible ship, will be particularly receptive to all the exis- 
ting test models and will stimulate development of more and 
more adequate and precise models. 

However, in the confrontation with reality that may appear 
to be wishful thinking. In fact, the ship design methods have 
been based for decades on the same hydromechanical solu- 
tions developed at the end of the 19th and beginning of the 
20th century. The later developed solutions, quite useful and 
definitely of better quality, have not been appropriately used in 
design. As that pertains to the basic and vital ship characteris- 
tics, such as transverse stability, load-lines, unsinkability and 
the resistance-propulsion properties, it is difficult to qualify 
the ship design as a modern and professional practical specia- 
lity. One can hardly expect its products to be, as a rule and not 
only by accident, as good as they can objectively be. 

Therefore one has to agree, though regretfully, that : 

»The naval architecture is a science of vague assumptions, 
based on debatable figures taken from inconclusive instruments, 
performed with equipment of problematical accuracy by per- 
sons of doubtful reliability and questionable mentality" 9. 

That opinion has to be accepted because in the more and 
more institutionalized and bureaucratic shipbuilding and ship- 
ping system there are no signs of will and action aimed at chan- 
ging the situation. On the contrary, quite clear are mechanisms 
and tendencies of maintaining the status quo in this respect. 

The problem discussed here is also important as it 1s rela- 
ted to the theory of ships - a research and practical science, by 
definition aimed at creating new more accurate and practically 
useful models of hydromechanical phenomena. Lack of de- 
mand from the practical design causes decreasing interest in 
creating such models and concentration on the in fact ineffec- 


tive polishing" of the hydromechanical solutions used in ship 
design procedures for years. That particular pragmatism” of 
the theory of ships closes the „positive feedback” cycle, which 
effectively petrifies and even deepens the hydromechanical 
backwardness of shipbuilding. 


Service margin 


The above characterised situation is reflected in predicting 
the ship service resistance and/or horsepower. Now, like in the 
past decades, the main design and research effort is concentra- 
ted on the idealised part of the problem, related only to the 
calm water conditions. Ship resistance and/or horsepower in 
real sea conditions is treated ,,marginally": they are not used as 
the hull size and shape selection criteria and their values are 
roughly predicted by means of a „service margin". The ,,servi- 
ce margin" method is the simplest of all the possible ways of 
predicting the real values of resistance Ета (Vg) and/or the 
ship main propulsion horsepower Рта (Ус) typical of the nor- 
mal sea conditions and the required service speed Vy. It con- 
sists in determining those values as a product : 


Ела (Vg) = (1+kw ) : Rr (МЕ) 
or: (1) 


Prs (МЕ) = (1+kw ) : Pr (МЕ) 
where : 


= ЕкУр) or P4(Vg) are the ship calm water resistance or 
horsepower values, respectively, for the speed Vg 

= k, is a dimensionless resistance and/or horsepower re- 
serve (margin) coefficient, whose values are determined 
by a rather subjective selection from the kw € (0.15 ; 0.35) 
range. 


In the past, up to mid 1960s, the so defined ,,service mar- 
gin" method was a justified practical design approach. At that 
time it was the only practical solution of the complex and im- 
portant problem, not yet investigated and, needless to say, 
modelled. Sticking to that method later when much more ac- 
curate and useful real sea resistance (horsepower) prediction 
design tools became available? has been a mere ship design 
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„Shop floor neglect". More importantly, that has been a neglect 
with an evident effect - the resistance and horsepower „оуег- 
size" ships. 

As it was demonstrated (e.g. in [1, 2, 3]), ships designed in 
the mid-twentieth century did not use in service ca. 2096 of 
their installed horsepower, whereas their average (taken from 
several dozen sea trips) service speed was in the range of 60% 
to 103% of the design speed. In the course of years, when the 
cargo ship design speed and horsepower increased quickly, the 
oversizing problem was becoming (it had to become) more 
and more clear. It gained importance also because ofthe signi- 
ficant improvement of ship hull construction technology.? 
Therefore, practical economic indications could be seen to ra- 
tionalize the design methods - to include the seagoing qualities 
into the ship hull design criteria or at least to reduce the service 
margins. It could have been expected that such revision would 
soon become a fact. Nothing of the sort has ever happened. 

Still a lot is being done in the ship design process to find 
a hull shape with a still water resistance lower by so much as 
a fraction of percentage point; the seagoing qualities of such 
ship are generally a second-rate matter, most often totally ne- 
glected. And still the so scrupulously minimized resistance is 
then roughly increased by a 15 to 35 percent service margin. 
Therefore, it may and should be assumed that even modern 
sea-going ships ,,carry” at least a dozen or so percent reserve 
of never used main propulsion system horsepower. 


The ship delivery-acceptance procedure 


Itseems that the reasons of such a deformed situation in the 
ship design resistance and/or horsepower prediction question 
should be sought first of all in the firmly established and insti- 
tutionalized ship delivery-acceptance procedures, where a par- 
ticular role is assigned to the ship speed (propulsion capabili- 
ty) test on the measured mile. 

That test, intended to be an objective verification of real 
propulsive and speed qualities of the ship, in fact does not pro- 
vide such a verification. It does stimulate, however, working 
out such design resistance-propulsion characteristics which 
meet the respective contract requirements in the test (still wa- 
ter) conditions. For the shipyard, negative result of that test 
means a sure financial loss determined in the contract. On the 
other hand, designing a ship having the resistance-propulsion 
characteristics well suited to the typical operating conditions, 
regardless of the measured mile test results, does not bring the 
shipyard (designer) any additional financial profits. 

As long as such a ship acceptance procedure is in use, one 
should not expect any significant changes (rationalizing) in 
designing the ship resistance-propulsion properties. Taking into 
account : 


© strong adherence of the shipbuilding and shipping industry 
to the generally accepted solutions, particularly those of 
the administrative and legal character 


c» objectively illusory but individually irrefutable conviction 
on the part of the owner of the necessity of checking the 
ship performance characteristics at the moment of delivery 


c» an evident interest of the shipyard and designer in applying 
a sanctioned service margin with its broad range of values, 
in order to ,,hide" all the inaccuracies of the resistance and 
horsepower predictions, 


it will take long for any changes in this respect to occur. There- 
fore, it may be assumed that the service margin will still be 
used in the ship design procedure for a long time. 


In this situation, what can and should be reasonably done is 
constant updating of the service margin values — improving 
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the accuracy of their ky coefficients. This task belongs to the 
domain of the theory of ships. It is an immanent duty of that 
speciality to generate more and more adequate models of all 
those hydromechanical phenomena which determine the func- 
tioning of ships. That duty is absolutely indispensable in rela- 
tion to the phenomena and problems connected with the ship 
safety and operational efficiency. 


The purpose and subject of investigation 


Apart from the above presented remarks, the purpose of 
this paper is also to present some real possibilities of determi- 
ning more adequate values of the service margin coefficient 
ky. Those possibilities in general consist in making an effecti- 
ve use of the seakeeping ability methods. In particular, they 
consist in generalization of the long-term prognosis of addi- 
tional ship resistance from sea waves Ray. 

It seems that the so obtained expressions of the type : 


PUE f (Fn, K) Q) 


T 


where : V. 


С Је 1. is a di- 
mensionless design ship speed and K represents the ship hull 
governing geometric characteristics, may be correct measures, 
if not the total ky coefficient then at least its very important 
component k - the service margin wave part coefficient. 


: Р : , Fn 
Ry - is the ship still water resistance, Ё 


One way or another, such expressions may be an important 
contribution to the necessary rationalization of the service 
margin values. 

The paper is based on the results of ship resistance tests 
carried out in the Faculty of Ocean Engineering and Ship Tech- 
nology, Gdansk University of Technology, in the ,,Balteco- 
logicalship" project and is an additional summary of those te- 
sts. The test materials are presented and discussed in detail in 
[5. 6, 7, 8]. 


LONG-TERM PROGNOSIS OF THE 
INCREASED RESISTANCE OF A SHIP 
MOVING IN HEAD WAVES 


Definitions 


> The ship resistance increase (additional resistance) was pre- 
dicted (determined) only for the head sea conditions i.e. 
for wave incidence angle B equal to 180°. 


> Thelong-term prognosis of resistance increase (additio- 
nal resistance from sea waves) Ray of a ship moving with 
constant speed V were determined as probabilistic estima- 
tion of a mean value of a set of random short-term additio- 
nal resistance RAw(x) prognosis for that ship : 


Ray = SR aw(x)-p(x) (3) 


where : 


= x= {a,b,c} is a set of discrete random sailing con- 
ditions, which have assigned to them the resistance 
increase values Raw(x) determined as short-term 
prognosis. Elements of the x set are: ship loading 
conditions [a], sailing regions [b(b,,b2)]—seas (61) 
and their regions (6) and sea states — wave condi- 
tions (c) characteristic of those seas 

= p(x)-p(a):p(b): p(c]b) is a probabilistic model 
of conditions x, or the combined probability func- 
tion of discrete random variables a, b, and c. 
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> The short-term prognosis of resistance increase 


Raw (x) (prognosis of resistance increase in the stationary 
irregular wave conditions) were determined in accordance 
with the principle of superposition expressed as follows : 


В ду (х)= 2f taw (o,a)S(m,b,c)do (4) 


[0] 
where : 


= Yaw (@, a) is a dimensionless resistance increase 
coefficient (resistance increase operator) determined 
from ship model tests carried out in regular head 
waves of different frequencies 00 

= $ (0), Б, с) is a model of the spectral density function 
of stationary irregular waves, the form of which was : 


1 
8,(0,T,,Hy)3)= 58 (o. T. Hus) Р 
. 3 goo [-80 206901, 1р3] 


where : S (œ, T1,H1/3) is a model of the ISSC spectrum. 


> Inthe third stage, using expression (1), the Ray prognosis 


were calculated. 

The Ray values for each ship and for each selected speed V 
were calculated by means of a probabilistic model of all- 
-year sailing conditions, whose elements are : 


= ship loading condition probability function p(a) presen- 
ted in Tab. 1 


= sailing area probability function р) presented in Tab.2 

= sea state probability function р(еђу) = p(T1,H1/3|b1) 
for the NORTH SEA presented in Tab.3 

= sea state probability function p(c|b1) = p(T1,H1/3|b1) 


— assigned to a sea region b». For restricted 
(coastal) waters — the JONSWAP spectrum, 
for unrestricted (open) waters — the ISSC 


spectrum | 
— identified each time by the sea state charac- ЈОНИ BATES Er рина 1н Ded 
teristics c = с(Нуз, T1) : significant wave Table 1 
Loading condition p(a; 
a [кту | E 


height Нуз and characteristic period Ті. 


» Theresistance increase prognosis were prepared for : 


омо | озю | тою | 


Table 2 


= shipping routes in the North Sea and Baltic Sea 
regions, whose probabilistic model p(b) was based 
on the Author's knowledge 

* random stationary waves, whose probabilistic mo- 
del p(c) = р(Ну/з , Ту) was based on the information 


зак [orase | ви 2 
taken from the ,,Global Wave Statistics" atlas [4] 0.150 0.100 | 0.250 
* three ships with two loading conditions [light (ај) 


and heavy (аз)] described by the model p(a) assu- 0.450 0.300 | 0.750 


med by the Author 
= three speeds V of each ship. The speed values were 0.600 0.400 | 1.000 
in the range Ve < 0.65Vg ; 1.1Vg > where Vg is Table 3 


the nominal (design) service speed of each ship. 


Periods T; [s] 


ERESERERERENES 
ee prse a a Тан 
0.048 | 0.120 | 0.096 | 0.037 | 0.010 | 0.002 | 0.316 
крип nsn [оа 
> In the first stage the functional relations: raw = f (Or, ) 
were determined, where : © = w(L/g)!” is a dimensionless 5 | - | 0.006 | 0.027 | 0.038 | 0.024 0.003 | 0.108 
wave frequency determined in a stationary system. s | - | 0.002 0.006 | 0.002 | 0.053 
The relations were based on the results of model tests car- 
ried out with 1:85 scale ship models towed in calm water s ЕГІ 0.004 | 0.009 | 0.007 | 0.004 | 0.001 | 0.026 
and on ca. 15 regular waves, whose frequencies @ were in 6.5 0.002 | 0.004 | 0.004 | 0.002 
Bore ү. КІН [| fof in өши 
T= fos oss [os [nom ons [oo [omi [ose 
> In the second stage the Raw prognosis were calculated in 
accordance with expression (4). Tables 
The Raw values were calculated ; 
for stationary waves (sea states) c (T1 , Нуз) : РОЗИ 


= whose Ту and Ну/з characteristics were taken from the sise[7|s[»| 


Tye < 45; 10s > , Нузе < 0.5m ; 6.5m > ranges 0.093 | 0.140 | 0.075 | 0.020 | 0.003 EX 
= whose spectral density function S (0, Ту, Нуз), depen- 
0.033 | 0.130 | 0.130 | 0.060 | 0.014 | 0.003 


Determination procedure 
of the long-term К ду prognosis 


The determination procedure of a single (for a given ship 
and its constant speed V) Ray prediction, as described by 
expressions (3) and (4), was the following : 


ding on sea region, is modelled by : 
- for coastal waters – the ISSC spectrum in the form: 


Es тіз -exp[- B-o-*) 


Si Tae 


| 35 [ооо 0.011 | 0.020 | 0.015 | 0.007 | 0.002 
Sel 0.002 | 0.006 | 0.005 | 0.002 
0.133 | 0.328 | 0.299 | 0.140 | 0.039 | 0.009 | 0.002 | 0.950 
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0.006 | 0.045 | 0.068 | 0.040 0.003 ШШ 


— for open waters - the JONSWAP spectrum in the 
form proposed by the Author : 
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hi | . ы 
Ships The Ray values were determined in accordance with pro- 
The resistance increase predictions were determined for three cedure described at the page 31, and the Вт, and Етс resi- 
ships, with their main characteristics presented in Tab.5. stance values were determined from model test results recal- 
Table 5 


NCUNITTTTIUETE 
Symbol 


Ship characteristics 


light | heavy light heavy light | heavy 
loaded | loaded | loaded | loaded | loaded | loaded 
ship ship ship ship ship ship 
Length between perpendiculars 132.00 132.00 147.75 
131.77 | 132.18 | 12928 | 132.00 | 145.24 | 147.50 


| 800 | 870 | 760 | 855 | 6.00 | 650 | 
ШІ 
radius of gyration i 
It can be seen that in the NORTH SEA : culated for a real ship by the three-dimensional extrapolation 
method. 


= annual mean wave parameters аге: T, 726.355; Нз = 1.95 m 
= calm water occurs with probability p = 0.031 i.e. approx. 
11 days in a year. 


It can be seen that in the BALTIC SEA : 


| ; ; Ray 
= annual mean wave parameters are: Тү = 5.365 На = 1.43m Figure presents a coefficient function k = 


в 1 t ith probability p = 0.050 i.e. й : Т 
is олы МОМИР SORS from the respective values taken from Tab.6 : 


Long-term prognosis of the relative increase 
of resistance in waves - wave service margin 


— f (Fn) built 


= points (different for different ships) mark the k = Ray 

Results of the calculations values taken directly from Tab.6 Кт 
= full line marks a universal function for the ships in question 
and their speeds, ofthe following mathematical form iden- 


Table 6 presents the values of : 


= mean still water resistance Ry = 0.8R4, + О2Етс (where tified here : 
Rr and Ryc are still water resistance values in the light Ray 0.0635 
and heavy loading condition, respectively) k == - 0.157 
= long-term prognosis of resistance in waves Ray R4 Fn (5) 
= total mean resistance in waves Ету = Rr + Ray where: Fn e «0.12 ; 0.30» 
А Г К А А 
= respective relative increases : --АУ- and zx The expression (5), with the assumption that : 
T TV 
Fn= Fn, = Vg 
determined for 3 speeds V 4g: L 
of each ship characterised in Tab.5. may also be treated as a certain model of a dimensionless co- 
V $ efficient (К) of the wave service margin, i.e. coefficient be- 
where: C, = s Gu DM ing a component (part) of the coefficient kw [see expression 
B L-B-d ' N L.B (1)], which determines the ship resistance and/or horsepower 
margin balancing the expected (in the long-term prognosis 
Table 6 sense) increase RAy(Vg) [or PAy(Vg)] caused exclusively by 


the sea waves. 
0.35 
12.5 | 164 | 20.3 


0.3 
Ad 
0.169 | 0.222 025 
5 02 
donani 
60.15 
О 


ZEJEXEJEI 
а је 


shape 


0.12 0.17 0.22 0.27 
Froude number Fn 


Ray 
Fig. Expression k = = f (Fn) 
Rr 
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An alternative version of model (5) is the expression : 


R 
k, йн -0.50 


Т (6) 
for Cp €<0.50; 0.85 > 


where: Сь = УБ is the design block coefficient 
"dg 
(corresponding to the design values of Vg and dg ). 
Model (6) was identified here 


by substituting to model (5) an expression : 
0.070 
ng = —— 
Св- 0377 
devised by the Author and very close (qualitatively 
and quantitatively) to the well known Ayre's formula : 
Св = 1.075 – 1.68Fng 

In order to illustrate the obtained results, Tab.7 presents 


percentage values of the k and ky wave service margin coeffi- 
cients, calculated : 


» according to models (5) and (6) 

> for ships characterised in Tab.5, where values of the Cg 
coefficients for those ships were determined as their long- 
-term prognosis, i.e. as weighted averages of the Cg values 
corresponding to the light (80 %) and heavy (20 %) loa- 
ding condition. 


Table 7 


carrier 


0.740 0.675 0.636 


Ук [kn] 14.00 16.50 19.50 
0.200 0.236 0.263 


k 0, 0, 0, 
kı 0, 0, 0, 


CONCLUDING REMARKS 


It is fitting to underline the following aspects 
in the summary : 


О The intention of this paper is to draw attention to the needs 
and possibilities of : 


Ж constant improvement of the quality and quantitative 
accuracy of the ship theory and design methods 

Ж more frequent application, in ship theory and in practi- 
cal ship design, of the useful but still little used models 
of ship hydromechanics - the seagoing qualities in par- 
ticular. 


O It was particularly intended to remind that such an impor- 
tant design problem as prediction of optimum ship resi- 
stance and/or horsepower in real operating conditions still 
remains a problem waiting for solution. 


O An important objective of the work was also to indicate 
a possibility of identifying more justifiable service margin 
formulae based on generalisation of long-term prognosis 
of the relative additional ship resistance in waves. 


O The identified models (5) and (6) do not aspire to be com- 
plete and directly applicable design service margin formu- 
lae, as they : 


** are based on too scarce research material, and also 
* are modelling only the wave part of the problem. 


O The above mentioned limitation of models (5) and (6) is to 
some degree compensated by the fact that their results are 
an upper estimation ofthe wave service margin. They were 
identified : 


* only for the head sea conditions when the resistance 
increase 15 the greatest 

* without taking into account the voluntary speed (and 
horsepower) reduction and/or change of the ship co- 
urse in order to protect the ship against too intense (dan- 
gerous) motions - accelerations, vibrations, slamming, 
shipping of water etc. 


O The form of models (5) and (6) shows clearly that the k 
coefficient depends either on the ship design speed Ук ог 
on the design values of block coefficient Ср: 


A itis inversely proportional to Vg 
A itis directly proportional to Cg. 
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a) Author of the sentence unknown 
b) There is a strong conviction that the work of research and 
practical sciences may be appreciated only when it is 
concentrated on fulfilling the existing and not potential design 
needs 
с) - „superposition principle” published іп 1953 by St. Denis and 
Pierson and its application to the resistance-propulsion 
problems published in 1967 by Vassilopoulos 
- standard spectral (spectral density function) wave models, e.g. 
the Pierson-Moskowitz model recommended in 1969 by ITTC 
- atlases of annual wave distributions in the main shipping sea 
areas, e.g. the „Осеап Wave Statistics" atlas published in 
1967 by Hogben and Lamb 
d) The problem is that the service margin is meant to be a remedy 
not only for the increased ship resistance (horsepower) due to 
sea-way but also due to the uncontrolled and gradually 
increasing hull roughness. There is no doubt that with totally 
welded hull shell (and constantly improved welding methods) 
and with a constant progress in the quality, durability and anti- 
-fouling effects of paint coatings, also the ,,roughness" 
component of the service margin should be revised (reduced) 
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ABSTRACT 


In this paper the predictions of resistance and propulsion for 4 types of ships: container carrier, oil product 

tanker, ro-ro vessel and river-sea vessel included in the EUREKA ,, Baltecologicalship " project E!2772 are 

presented. All the above ships were designed by the SINUS Enterprise Ltd. The scope of the research 

covers: experimental measurement of ship resistance in calm water, determination of resistance increment 

in irregular waves and computational powering prediction for POD propulsors. The experiments were 

conducted in the Laboratory of the Department of Underwater Technology, Hydromechanics and Design of 
Ships, Faculty of Ocean Engineering and Ship Technology, Gdansk University of Technology. 


Keywords : model tests, ship resistance in calm water, ship resistance in irregular waves, POD propulsors 


INTRODUCTION 


Trustworthy determination ofthe resistance characteristics 
of the designed ships through model experiments is necessary 
for confirmation of the hydrodynamic quality of design. De- 
termination of the full scale resistance enables realistic predic- 
tion of the propulsive characteristics of ships, assessment of 
the delivered power and propulsive parameters. Minimization 
of the engine power and ensuring correct conditions for engi- 
ne-propeller interaction is an important problem both from the 
economic (operational cost) and ecological (level of pollution) 
points of view. The above mentioned problems refer to the task 
"Resistance Model Experiments for 4 new types of ship", in- 
cluded in the EUREKA project E!2772. The Project covers the 
following ship types : 


* Container carrier 
+ Oil product tanker 
* Ro-ro vessel 

* River — sea vessel. 


All the above ships were designed by the SINUS Enterprise. 
The scope of the research covers : 


A experimental measurement of ship resistance 
in calm water 

A experimental determination of resistance increment 
in irregular waves 

A analysis of the ship-induced wave system 

A computational powering prediction 
for POD propulsors. 


The experiments were conducted in the Laboratory of the 
Department of Underwater Technology, Hydromechanics and 
Design of Ships, Faculty of Ocean Engineering and Ship Tech- 
nology, Gdansk University of Technology. 


RESULTS OF EXPERIMENTS 


The ship hull models for experiments were built in the model 
workshop on the basis of the theoretical hull lines supplied by 
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SINUS. Resistance experiments were performed in the towing 
tank and visualisations of the flow on the hull were made in the 
circulating water channel. Ship models were tested for two dis- 
placements: design draught and scantlings draught. The resi- 
stance prediction for full scale was performed using three-di- 
mensional extrapolation with form factor (1+k) calculated ac- 
cording to the ITTC-72 formula. The approximate computa- 
tional powering prediction was based on experimentally deter- 
mined resistance, on the propeller-hull interaction parameters 
determined by the Holtrop method, on the design POD propul- 
sor data supplied by SINUS and on propeller design software 
based on lifting line/lifting surface theory. Prediction of ship 
resistance in irregular waves was done for the Baltic and North 
Sea wave spectrum. The results of all experiments and compu- 
tations were the basis for assessment of the hydrodynamic qua- 
lity of each ship design. 


Container carrier — SINE 202 


The basic parameters of the container ship 
SINE 202 are as follows : 


> Length overall 138.10 m 
> Length between perpendiculars 132.00 m 
» Breadth 22.50m 
» Depth 11.20 m 
» Design draught 7.60 m 
> Scantlings draught 8.55 m 


The hull geometry is presented in the form of body lines in 
Fig.l. The model of the container ship was manufactured to 
the scale 1: 84.5. The resistance characteristics obtained for 
design draught and scantlings draught, recalculated into full 
scale are shown in Fig.2. The calculated effective engine po- 
wer curves for both draughts are presented in Fig.3. 

On the basis of experimental results it may be concluded 
that the ship will reach the speed of 17.80 knots on trials at 
design draught, 17.20 knots in service at design draught and 
16.25 knots in service at scantlings draught. Further increase 
of ship speed by 0.2+0.3 knots may be possible by increasing 
the propeller diameter by about 10%. 


The predictions of resistance and propulsion for 4 types of ships 
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Fig. 1. Body plan of the container carrier SINE 202 
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Fig. 2. Total resistance of the container carrier SINE 202 
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Fig. 3. Effective engine power of the container carrier SINE 202 


It is assessed that the container ship SINE 202 has good 
sea-keeping quality. The increment of resistance in waves is 
relatively low both in the sense of long-term predictions (Fig.4) 
and in the sense of short-term predictions (Fig.5). It has been 


determined that the sea margin of 15%, applied in propulsive 
prediction, 15 correct for typical sea states in the Baltic and 
North Sea. Additionally, it has been determined that the resis- 
tance increment in waves is identical for the design and scan- 
tlings draught. 
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Fig. 4. Resistance characteristics of the container carrier SINE 202 
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The detailed results of experiments 
and computational analyses are included in the report [1 , 2]. 


Oil product tanker - SINE 203 


The basic parameters of the oil product tanker 
are as follows : 


> Length overall 138.10 m 
» Length between perpendiculars 132.00 m 
» Breadth 22.50m 
» Depth 12.80 m 
» Design draught 8.00 m 
» Scantlings draught 8.70m 
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Fig. 6. Body plan of the ой product tanker SINE 203 


The hull geometry is presented in the form of body lines in 
Fig.6. The ship model was manufactured to the scale 1: 84.5. 
The resistance characteristics obtained experimentally for de- 
sign and scantlings draughts respectively are shown in Fig.7. 
The calculated effective engine power curves for both drau- 
ghts are presented in Fig.8. 

On the basis of experimental results it may be concluded 
that the ship will reach the speed of 12.80 knots in trials condi- 
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Fig. 7. Total resistance of the oil product tanker SINE 203 
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Fig. 8. Effective engine power of the oil product tanker SINE 203 


36 POLISH MARITIME RESEARCH, Special issue, 2004 


tion, 12.25 knots in service condition at design draught and 
12.05 knots in service at scantlings draught. 

Model experiments of resistance in waves have shown that 
the SINE 203 model has a small increase of resistance in wa- 
ves both in the sense of a long-term prediction Fig.9 and in the 
sense of a short-term prediction Fig.10 for sea states typical of 
the Baltic and North Sea. 

The detailed results of experiments 
and computational analyses are included in the report [3]. 
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Ro-ro vessel - SINE 204 


The basic parameters of the Ro-Ro vessel are as follows : 


2 Length overall 156.72m 
2 Length between perpendiculars 147.75 m 
2 Breadth 24.80 m 


The predictions of resistance and propulsion for 4 types of ships 
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The hull geometry is presented in the form of body lines in 
Fig.11. The ship model was manufactured to the scale 1:100. 
The resistance characteristics obtained experimentally for de- 
sign and scantlings draught respectively are shown in Fig.12. 
The calculated effective engine power curves for both drau- 
ghts are presented in Fig.13. 
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Fig. 12. Total resistance of the ro — ro vessel SINE 204 
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On the basis of experimental results it may be concluded 
that the ship has favourable resistance characteristics. The plan- 
ned value of delivered power is sufficient for attaining the plan- 
ned design speed of 20.0 knots. A minor power deficit for scan- 
tlings draught with sea margin is meaningless from the practi- 
cal point of view. The ship will reach the speed of 20.92 knots 
in trials condition, 20.45 knots in service condition at design 
draught and 19.90 knots in service at scantlings draught. 

Model experiments of resistance in waves have shown that 
the resistance increment in waves is relatively low both in the 
sense of a long term prediction Fig.14 and in the sense of 
a short term prediction (Fig.15). It has been concluded that the 
sea margin of 15%, adopted in the propulsion analysis, is cor- 
rect for sea states typical of the Baltic and North Sea. 
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The detailed results of experiments 
and computational analyses are included in the report [4]. 
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The basic parameters of the River-sea vessel are as follows : 
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Fig. 18. Effective engine power of the river - sea vessel SINE 205 
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River-sea vessel - SINE 205 
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The hull geometry is presented in the form of body lines in 
Fig.16. The ship model was manufactured to the scale 1: 50. 
The resistance characteristics obtained experimentally for sea 
and river draughts respectively are shown in Fig.17. The cal- 
culated effective engine power curves for both draughts are 
presented in Fig.18. 

On the basis of experimental results it may be concluded 
that the ship will reach the speed of 11.40 knots in trials condi- 
tion at sea draught, 10.90 knots in service condition at sea drau- 
ght and 12.00 knots in trials at river pon 
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Fig. 16. Body plan of the river - sea vessel SINE 205 : 
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The detailed results of experiments 
and computational analyses are included in the report [5]. 


CONCLUSIONS 


The above presented predictions of resistance and propul- 
sion for 4 types of ship obtained through model experiments 
enable realistic prediction of the power and propulsive charac- 
teristics. In the cases where the performance required by the 
EUREKA project was not achieved, the directions for neces- 
sary modification were indicated. The positive results of expe- 
rimental investigations demonstrate that the shape of the hull 
geometry and basic parameters were well selected and they 
create a possibility of ship design with favourable hydrodyna- 
mic characteristics. 


NOMENCLATURE 
Нуз - significant wave height [m] 
Р. - effective engine power [kW] 


Кач  - total resistance of ship [kN] 
Кту - long-term prediction of seaway resistance [kN] 


Т - characteristic wave period [s] 
Vs - Ship speed [knots] 
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Stability criteria as constraints 
in a fleet of ships optimisation problem 


Boguslaw Oleksiewicz, D.Sc.,Eng. 


ABSTRACT 


The paper has been written within the European EUREKA Project E!2772, initiated and 
completed at the Faculty of Ocean Engineering & Ship Technology, Gdansk University of 
Technology in the years 2001-2003. A problem has been solved concerning mathematical 
optimisation of a fleet of multipurpose sea-river vessels for European short-shipping regu- 


M 


lar lines, in the area of The North and Baltic Seas, on the level of marine transportation 
task, by the non-linear programming methods with constraints. A method is proposed which 


enables existing criteria of stability to be included as constraints in the optimisation model 
ofa fleet. In the numerical examples, three typical criteria of intact stability: by IMO, PRS, and HSMB have 
been selected to demonstrate a post-optimisation feasibility analysis of principal parameters of ships. 


Keywords : maritime transportation, computer-aided ship design, optimisation, intact stability criteria 


1. INTRODUCTION 


Computer aided ship design methods used at present, whi- 
le offering automation of the design process, require its ratio- 
nalisation and formalisation. In consequence, adequate mathe- 
matical models of the design object must be created which af- 
fect the design process by introducing a structure and termino- 
logy which unavoidably bounds reasoning to the terms of the 
model. 

In this case a fleet of ships at the stage of owner's study 15 
assumed to be an object and the task of optimising its main 
parameters is an objective of the adequate mathematical mo- 
del. In consequence, the global structure of the model (further 
called an *optimisation model") corresponds to that proposed 
by operational research methods in general and non-linear pro- 
gramming methods (NPM) in particular [3]. Within this struc- 
ture, optimisation models consist ofa set of sub-models of par- 
ticular properties of the object which have been recognised as 
significant to the predictive features of the model. Optimisa- 
tion models applied to fleet/ship design are definitely synthetic 
in nature. This feature requires the analytical representation of 
particular sub-models to be relatively simple. In consequence, 
sub-models usually neither become isomorphic with, nor con- 
form to the physical structure of that part of object to which 
they are related. Such type of models is sometimes referred to 
as “non-structural” [19]. NPM require for all the sub-models 
concerned to be formulated as constraints. Among them there 
are always those concerning safety of an object. In ship design, 
a special interest in this group is focused on the stability of 
ships. In naval architecture today, the stability requirements 
are imposed in the form of legal regulations by such institu- 
tions as IMO, classification societies, governmental organisa- 
tions and other bodies. An essential part of stability regula- 
tions are stability criteria. 

The paper deals with the problem of incorporating stability 
criteria as constraints in the optimisation model of fleet/ship 
design. At the initial stages of the design the principal difficul- 
ty is that the full geometry of a hull, necessary for the stability 
criteria to be applied, is usually unknown. A standard solution 
was to take into account the initial stability only, represented 


by the initial metacentric height GMO [2], [4], [7], [12], [14]. 
The paper proposes an alternative approach, based on an idea 
introduced by Wisniewski [20] and developed by Kupras [10], 
[11]. In this concept the full stability of ship can be accounted 
for by using systematic standard series of hull forms, follo- 
wing the methodology developed in ship resistance and power 
prediction. 

In order to accomplish the task, an attempt has been made 
to define all the stability-related geometrical characteristics of 
a ship analytically, based on the Series 60 body forms [19]. In 
consequence, an arbitrary criterion of intact stability can also 
be defined in an analytical way and so incorporated into opti- 
misation model as a constraint. 

Stability aspects in the computer-aided modelling of ship 
design have been addressed on the background of the optimi- 
sation problem concerning a fleet of multipurpose sea-river 
vessels for European short-shipping regular lines, in the area 
of The North and Baltic Seas, on the level of marine transpor- 
tation task, by non-linear programming methods with constra- 
ints. The problem has been undertaken within the European 
EUREKA Project [13] based on predictions that a significant 
increase of cargo transportation in Europe over the next 10 
years (or probably after this period) will take place between 
Western Europe and the Central and East European countries. 

In the numerical examples, three criteria of intact stability: 
IMO [6], HSMB [5], and PRS [17] have been selected, as typi- 
cal of contemporary stability regulations, to demonstrate the 
method in a post-optimisation, feasibility analysis of principal 
parameters of ships. 


2. PROBLEM STATEMENT 


A fleet of ships consists of a number of homogeneous ships 
operating as a maritime transportation system in a certain envi- 
ronment. The transportation task for a fleet of ships is to carry 
goods between ports during a prescribed period of time. An 
optimum fleet to perform this task, given the particular (ow- 
ner's) data, is a general problem under discussion. A solution 
to this problem needs adequate functional and mathematical 
models. 
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2.1. Functional model of a fleet 


In the particular case (Tab.2.1), a (potential) shipping line connects the furthest Western and Eastern regions of Europe (a). 
A corresponding model of shipping (b) is called a multi-port route model linking two areas of operation A and B with the two 
groups of clustered sea and hinterland river ports. There are two streams of goods transportation in the model: from A to B (called 
OUT) and back, from B to A (called IN). Ports A-0 and B-0 are the home and destination ports. For more details about the 
functional model of a fleet - see [13]. 


Tab. 2.1. An example of a shipping line and its graphical model 
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2.2. Optimisation model of a fleet 


The mathematical model chosen for the fleet optimisation problem can be described as deterministic, static, continuous, 
single level and single objective model, imposed and solved by non-linear programming methods. A standard formulation of the 
model, within NPM, is as follows: given a vector c (or a set C) of constants, find such a vector of decision variables x that 
minimises a single valued objective function Q(x,c) subject to a set of inequality constraints. An adequate mathematical form of 
the problem is : 


min О (хе); Q:R'xR^ >R! 


(х;е) e О = (oco) :2 (хе) x0 j= 12,..m}# @ ы) 


It is generally assumed that Q(-) and gj(-) are all non-linear functions. 
The conditions for existence and uniqueness of the (optimum) solution to (2.1) can be found in [3]. 


Tab. 2.2. Elements of an optimisation model of a fleet of ships 


Itis obvious that the optimum solution x of the problem (2.1) is, in fact, parametrised by constants C that can be classified 
according to different criteria (Tab.2.2a). In particular, the group C3 delivers constants for the stability criteria of a ship (legal 
environment). A crucial element of the optimisation model is a feasible solution region (FSR) Q, because it eventually decides 
about an optimum solution to the problem. FSR is such a set of pairs (x,c) that all the constraints hold. 


Q = ((x,c):gi(x,) <0, je J}, J= {1,2,...m} (2.2) 


For further discussion, it is useful to classify FSR from the functional and formal points of view. In the functional classifica- 
tion, FSR defines, in fact, the notion of feasibility in ship design. Formally it can be thought of as a common part of the four 
groups of requirements imposed on a fleet/vessel by the environment (Tab.2.2b) and can then be written as the following product 
of four sets : 

4 


Q= f 


k=1 
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АП of them affect constraints of the model. For example, the SAFETY (Q4) group contains, among other things, design 
restrictions concerning stability regulations. 

In (2.2), an individual constraint is an inequality-type relation in which a function gj(-) expresses a balance between certain 
(dependent) parameters of an object. In order to keep the same standard relation of inequality (X) for all constraints, the general 
form of gj(-) has to be alternatively : 


(Се) - руке) (а) 
5 (0) = Е Ы кии ог jeJ (2.3) 
руске) - pjGce) (b) 


where : 


the p;(-) parameters are those predicted by the model and pj(-) are the corr esponding ones required by Оқ. 
In most cases the pj parameters аге just constant figures such that p;7c;eC. 


In the formal classification of FSR, attention will be focused on constraints forming the boundary of Q. 
Figure in Tab.2.2c illustrates the problem. Let us define an FSR corresponding to a single j j constraints : 


О = {(x:c):g (ху) = 0) , јеј 
Its boundary is then: 


dQ, = (хе): (ко 20]. jeJ 


For the whole feasible region one has, of course: Q = ПО, , jeJ 
Similar relation does not, however, hold for boundaries : 


09 z(]0Q;,2U0Q9; but IQ CUIQ,, je J (2.4) 


A је constraint is called a boundary constraint if its own boundary contributes to the boundary of Q, otherwise the constraint 
is a non-boundary constraint. It is obvious that, in fact, the boundary constraints are those which determine © and, as such, affect 
the optimum solution of the problem : 


Q= По; = ПО; (2.5) 
jeJ је JgcJ 


Relation (2.5) suggests an apparently equivalent formulation of (2.1) with the smaller number of constraints (boundary 
constraints only). Unfortunately, this is not the case because the set of boundary constraints Jp is not known in advance. Among 
the boundary constraints one can further distinguish active constraints and passive constraints. The classification refers to the 
location of the optimum x* solution in Q, which exclusively depends on the objective function Q(-) used in the optimisation 
process. Ај constraint is active in x* if x* є 0Qj (= је Jj, ), otherwise the constraint is passive (= je ЈЕ = Jg \ Jẹ ). If the 
optimum solution belongs to the interior of FSR (x* € intQ), all the boundary constraints are passive and no constraint affects x*. 

Let us now introduce another classification of Q and the corresponding constraints : stability constraints vs. non-stability 
(remaining) constraints. 

O-O,nQ, e 7-09), (2.6) 

where : Qs is a feasible solution region with regard to stability requirements, Js is a set of indices of stability constraints, Әр, 

Jg. correspond to the remaining constraints accordingly. The current classification is independent of the previous one; i.e. both 

stability and remaining constraints can become the boundary or non-boundary constraints. Let us investigate the status of stabi- 

lity constraints via the relations: Qs vs. Q, ОО; vs. 0Qs, and Js vs. Jp. It is obvious that the relation Qs > Q is always satisfied, 
but similar relation does not occur for boundaries. One can distinguish here four cases (see Tab.2.3). 


Tab. 2.3. Stability constraints (gs) vs. remaining constraints (gp) 


ad. a) Stability constraints do not contribute to the FSR, so, for given stability regulations and a set of constants C, they are 
totally insignificant, regardless of the objective function Q. 
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ad. b) Stability constraints do contribute to the FSR, so, for given stability regulations and a set of constants C, they are currently 
insignificant but potentially significant, dependent on the objective function Q. 

ad. c) Stability constraints do contribute to the FSR, so, for given stability regulations and a set of constants C, they are currently 
significant, but potentially insignificant, dependent on the objective function Q. 

ad. d) Stability constraints form the boundary of the FSR, so, for given stability regulations and a set of constants C, they are 
currently and potentially significant, regardless of the objective function Q. 


The foregoing discussion will be recalled in Chapter 5 to verify the status 
of stability constraints in the fleet optimisation model. 


3. IDENTIFICATION OF INTACT STABILITY CRITERIA 


In order to illustrate the thesis of the paper, three criteria of intact stability have been chosen, as recommended by the 
following institutions: International Maritime Organisation (IMO [6]), Polish Register of Shipping (PRS, [17]), and Hamburg 
Ship Model Basin (HSMB [5]). In the light of the discussion in [8], all these criteria can be regarded as a present day standard as 
far as intact stability regulations in naval architecture are concerned. 

Tab.3.1 introduces the coordinate systems as well as basic notions and terminology which will be in use throughout the paper. 
In particular, (Y-Z) and (y-z) are, respectively, a dimensional [m] and non-dimensional [-] coordinate systems of the hull and 
(w-t) [-] is a non-dimensional coordinate system rotated around the x axis. 


Tab. 3.1. Coordinate systems and basic terminology of intact stability of ships 


ЛЕ + ІМ (3.1a) | 

1171 +12 | 
НЕ соѕф “Н (3.1b) | .. | 
t -sinó соѕф | | 2 ЈА 


La($:Zg)-lg($:za).T [m] (32а) 


(фас) e Wa ()-Zg ‘sinp [H] G20 | 15 2 $ * * PRS 


A common feature of the criteria under discussion (in fact, all the contemporary criteria of intact stability of ships) is that, as 
far as the righting moment of a ship is concerned, they all are totally based on the righting lever curve calculated in a calm water 
as the function : Lg ($&:D, Za) (3.2a), where : ф is angle of inclination and D [t], 2с [m] аге a constant displacement and 
coordinate of the centre of gravity G, describing current loading condition of a ship. In further discussion, the only loading 
condition accounted for is the design condition. It follows from (3.2) that all the stability-related characteristics and parameters 
of a ship have been normalised with regard to the design draught T. A short overview of the criteria by IMO, PRS апа HSMB in 
an analytical form has been shown in Table 3.2. 

It can be seen that the criteria by IMO and PRS involve both righting Lp(-) and heeling Ly(-) lever curves (,,weather-type” 
criteria), whereas the criteria by HSMB are based on the righting lever curve exclusively (,,Rahola’s-type” criteria). 

It follows from Tab.3.2, that, from a mathematical viewpoint, all the criteria can be formulated as combinations of certain 
operations on the functions Lg(-) and Ly(-) : linear - such as calculating ordinates, first derivatives, and integrals or non-linear – 
— such as calculating characteristic angles, the weather parameter K, etc. Two examples illustrate the problem : 


. Cross curve 
. righting lever 
curve 


(i) Linear case - a dynamic righting lever curve Тло) (PRS) is defined by integration of Lp(-) : 


–ф -Фо 

[Ол аф + [La ф)0ф for ф<0 

L4, (ф) = n ' 
РОТ) for ф>0 
0 


(3.3) 


(ii) Non-linear case - angles of the first and second interception фу, ф (IMO, PRS) are defined 
as the smaller and greater roots of the following non-linear equation : 


F9) = Lu (9) - Lg (9) =0 < $, = тіп[Е '(0)], ф, = max[F (0)] (3.4) 


A formulation of the intact stability criteria in Tab. 3.2 differ from but are fully equivalent to those originally formulated in 
the referenced documents by IMO, PRS and HSMB. Among other things, a notation has been unified and sequences of the 
righting and heeling lever curves as well as characteristic angles have been introduced to emphasize the analytical aspects of the 
criteria. 

In the subsequent chapters of the paper, the intact stability criteria will be examined with a special attention to those parts of 
them which concern the righting lever curve of a ship. This is because determination of the righting lever curve lp(-), via cross 
curves wp(:) involves the full geometry of the hull and as such decides on the reliability of the whole model. 
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Tab. 3.2. Intact stability criteria by IMO [6], PRS [17] and HSMB [5] 


Righting / heeling lever curves Stability criteria 


1. Lg (0) 20.15 m А 
2. [ (25) > 0.20 m 6. [Lg (ф)аф 20.30 m 
30 


Lg (30) 20.20m пф 
30 

. IBEX) аф > 0.055 mrad 8. фк 20; 
0 


9. ф, > 60° 
bs 
. [La (ф)аф 20.090mrad 10. Комо 21.0 
0 


Dus ш 226 . Lp (0) 20.20 т ‚бф, 20.30 
J- IN Lan) 2020 т ‚фк <ф 
ТІНЕ у = T i Rm K С 
BM TED · ф, 260 - Ab, 
| S а bm, 225-0.5- A0, ZONE 


a 30 E DE 60 & ф,-Фүдесі 


0.020 
. Le (30) x: j fu (ф)аф > 


INO ns ЕТЕ 


0.035 


0.12 
3. fr (ф)аф BU Ты (ф)аф > 


Characteristic angles 
фо - initial heeling angle, 6, - rolling amplitude, фу - ship flooding angle, Фф, - deck immersion angle 
$, - turning angle, ф, - superstructure flooding angle, p. - capsizing angle, Qc = min (фң:,%,) - critical angle 


фу -the angle for calculating K, Om, ; mn, -1* and 2" maximum of Lg angles, ф,, maximum of L angle $,, = Pm, 


m 
if Lr (фу ) Z Lr Cm, ) otherwise Фф, = Pm, . фу. Q2- first and second interception angles, фу = min (фр ,40) 
фа =min(0.8- 4,16). ф5 = шіп(ф;,0;,50). фе = max(0,0,). фу = min (0.5 - 4,15) 

Weather criteria indices: Kiyo = Ag / Ан [-], Кр = Lu, /Ly, [-] 

HSMB form factor 


Св. 100 2b-B 2l, 
· |— , where H'= H+h- :—— - а corrected depth, h - height of a hatch above 
zg (Сұ L B Ӯ 


deck [m], b - breadth of a hatch (b > B/2) [m]. |, - length of hatches within 0.5 L 


4. ANALYTICAL FORM OF STABILITY-RELATED CHARACTERISTICS OF A SHIP 


Analytical definition of the stability constraints needs the analytical definition of righting and heeling lever curves. For the 
reasons given earlier, the attention will focus on the righting lever curve only. This will be done by approximation or interpola- 
tion of a series of stability-related characteristics of a ship. 

The geometrical data for the task were prepared in [1], based on the systematic calculations of cross curves for Series 60 (S60) 
[19] for three different block coefficients Св and four h = H/T ratios. Three parent models of S60 had been chosen and then 
modified by: (1) extrapolating sections to obtain different h ratios (beyond the basic h = 1.50), (11) adding standard superstructu- 
res (a poop and forecastle), (iii) extracting the data concerning merely the design draught T (displacement), and (iv) making them 
dimensionless. The ratio b = B/T =b —2.50 was kept constant for all the models. Table 4.1 presents the resulting twelve models 
and, as an example, geometry of the M7072 model. 

Table 4.2 shows a full list of ship characteristics that have to be defined analytically in order to form the stability constraints 
in the fleet optimisation model. All of them are dimensionless (related to T) and classified from the analytical point of view as 
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one - (1D), two - (2D), three - (3D), and four-dimensional (4D) characteristics, taking into account a number of variables of 
corresponding interpolating or approximating functions. 


Tab. 4.1. The Series 60 hull models used for the approximation of cross curves and geometry of the M7072 model 


|. дф оу | — 060 | мв | ме | ме | М6090 | 


Tab. 4.2. Stability-related geometrical characteristics of a hull form and their analytical representation 


MD eem [ште | | m _ 
No Characteristics Drocess Function 


Cross curves in the (w-t) system Wp = #;(ф; Cg, h, b) 
[e| — — Riding ever cue сава) 


Fig.4.1 presents three examples of 2D/3D geometrical characteristics 
of a hull form based on the S60 data defined by the interpolation method. 


ет. 


~ b) c) 


20 00 06 сәз 


Fig. 4.1. Graphs of analytically defined 2D / 3D geometrical characteristics of Series 60 by the interpolation method 


In the analytical definition of stability constraints, the basic characteristic is righting lever curve Ig(0;:). 
A fleet optimisation model requires for the righting lever curve to be defined as the following 4D function : 


lk (ФС h.b;za) = wa (0.Cs. hb) - 26 -5іпф (4.1) 
The process of generation of the righting lever function Ір(Ф;) (4.1) has been divided into five steps : 
(1) Approximation of cross curves based on the twelve models of S60, resulting in a 3D function Wh (6, C4, h) 


for b = b? = const. 


(п) Correction ofthe Wh (-) function for the regions of small h (below the available S60 data), resulting in the wR ( -) function 
for these special regions. This step was necessary to cover H/T ratios typical of sea-river vessels and expected when optimi- 
sing a fleet of such ships. 


(iii) Affine transformation of the function s (-)/ wp ( -) with regard to b = b, resulting in a 4D function wh (6, Cg, h, b) 

(iv) Definition of the righting lever curve as a 4D function: lp (ф, Cg, h, b) = wp (6, Cg, h, b) - 26 - sino 

(v) Approximation ofthe АП function lp (ф, Cg, h, b ; Zg )bya 1D function (of 6 only): lp (0; Cg, h, b, Zg )for Cg, h, b, Zg = const. 
In all the steps listed above, the only problem that appears twice is that of approximation of multivariable (1) or single- 


variable (v) functions. The problem will be then addressed first as a separate numerical problem formulated in a compact matrix 
notation. 
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4.1. Approximation problem 


Let f (u) be a function to be approximated and Ға) be its approximation, both assumed to be single valued, 
multi-variable functions of u = (u,,u,,....uy). The approximation problem is formulated ав: 


pf, f) — min (4.2) 


where : P(-,) is a metric ofa proximation. Let the f(-) function be given in a discrete manner 
by the following sequence of data : "cf E ., and the ҒО function has the following linear representation : 


f(u;a) = Уо, e,(u) - a! е(а) , n<N (4.3) 


where : @ = (O4, 0t, ,....0t, ) a vector of unknown coefficients, e(u) = le, (и), e; (и),...е, (а) - a basis. 
For given arguments one has: fh - f (u.a), k-212,.,N and f= (f, на fu)" 
Using the representation (4.3), one gets the following equation : 
f -El.a (4.4) 
A matrix E = је ЈЕ [е; (0 )],, Will be later called the basis matrix. As the metric 


of approximation f(- , -) a weighted norm | s |l of an Euclidean distance between f and f has been chosen, then : 


o(f,f) = |г-Ё | (4.5) 
2w 
In consequence, the approximation problem (4.2) can be formulated as the following quadratic optimisation problem : 
тіп Q(a) , Q:R^ К! (4.62) 
with the objective function : 
оо) = |е] = (f —E'-a)'-W-(f-E'-a) (4.6b) 


where : W = diag(w) , w = (W1, Wo, ..., we a vector of weighting factors (уу; 20, 1=1,2,....,№). 
Now let us assume that some linear interpolatory constraints are imposed on the f (-) function : 


LODS, ‚ и=1,2,....т<п=йтпа (4.7) 
where I,,(-) is a functional of f(-) in the general form : 
I, (f)-L,(f)u,;a) . u-12... m (4.82) 


апа ГС) is a linear aa pui instance, Ly(-) can be a differential operator, as below : 
E 
È «Ә-Ә C 5 Au =0,1,2..., iu 212.3..5 (4.8b) 
u 
диг 
u 
where the sequences of indices: {Ан | И liu , U=1,2,....,m have to be defined separately. 
Let ГТ) - LENLE Joal) be a vector of functionals of ШӨ, . 
The interpolatory constraints can then be written as the following matrix equation : 


I(f (a) =" Pa 
where : €= (C,,C5,....C,,)' is a vector of given values of interpolatory constraints. Thanks to the linearity of I(-), one can write : 
If) = Y" -e(u))- a" -1((u)) - o -C7 (4.10) 


A matrix C= lew |= [L (e, Yu, )] will be later called the constraint matrix. 
Putting (4.10) into (4.9), one obtains the interpolatory constraints as the following matrix equation : 
Н(о)- С-а-с-0 (4.11) 


The approximation problem with the additional interpolatory constraints 
can now be formulated as the following quadratic optimisation problem with constraints : 


min Q(a) = (f E! -a)'-W-(f -E! -a) 
a € О = la: Н(о) = 0, dim0 = m Сп) c R” 
It has been proven [3], that both problems (4.6) and (4.12), as the convex problems, do have unique solutions. Moreover, in 
both cases the solutions can be found analytically. In particular, for the problem (4.12), this can be done either by Lagrange's 


multipliers or by the penalty function method. The latter approach has been applied in the case. Let the penalty function РС) be 
defined as a square of the norm | Н(а)|, in the form : 


P(a) = H' (a): 0- H(o) (4.13) 


(4.12) 
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where: œ = diag(0), © = (0, 092, ..., Ор) - a vector of weighting factors for the penalty function, œ. 2 0., i=1,2,....,m. 
The optimisation problem with constraints (4.12) can now be replaced by the equivalent optimisation problem 
without constraints, with the modified objective function ЕС): 


min F(a) = О(а) + P(a) , P:R" 2 К! 


deR” (4.14) 
The solution to the problem (4.14) follows from the necessary condition of a stationary point in R” : 
У.Ео)-0 © У,О(а)-У,Р(о)-0 (4.15) 
Finally, one arrives at the equation : 

(E- W-E! «C! .9.C).a -(E- W.f 4C! -@-g)=0 (4.16) 

| BO — b —— 

or, in short : 

B-a-b=0 (4.17) 


Existence of the solution to the problem (4.12) assures that the matrix B is non-singular (det B = 0). 
So one eventually gets the desired definition of the approximating function f (-) : 


а= Bb (4.18) 


One can observe that: (1) the solution 0 is parametrized by the weighting coefficients: & = a(w;œ) and (ii) putting œ = 0, one 
gets the solution of an unconstrained problem (4.62). 

From the numerical point of view, solution à to the approximation problem (4.2) can be found either by solving the system 
of linear equations (4.17) or by solving a single matrix equation (4.18) (inverting B). The first approach has been applied in 
approximation of cross curves (Section 4.2) and the second one in approximation of righting lever curve (Section 4.4). 


4.2. Approximation of cross curves 
In this case the approximated function is a cross curves function (f = w9) determined in a discrete way for the twelve S60 


models listed in Tab. 4.1. Given 11 ordinates of the function per one model, one gets N = 132 ordinates to be approximated. The 
following polynomial three-linear form has been chosen, as a particular representation of the approximating function f (-). 


f (uo) =a"-e(u) = У are (и) = У, о ге (фСућ) = У, о (ht Ср“) = (а) 
m m а= 


ul e, (u) (4.19) 
LJ K 2: 

E УУ Уакф "у Св“ (b) 
i-0 j-0 k-0 


The equivalency of the two notations introduced follows from the correspondence of both indications. After some trials and 
errors, the upper limits of the indices in (4.195) have been established as follows : I2 5,J 2 2, К = 2. It gives the following 
dimension of the basis e (^) in (4.3) : dime () =n = (1+1): (J - D)- (K +1) = 54, so the condition n < N holds. The basis matrix E 
(4.4) can now be readily calculated according to the definition : 


E =le = [69.0.5]. tor 1812 » 
Шу v v B Jy v=1,2,....,N (4.20) 
The constraint imposed on the approximated function Wh (ф,:) follows from the known property of cross curves : 
d — 


dó 


It imposes the following interpolatory constraints on the approximating function f () : 


д ~ 
> Са] Lo Za C3) аз» 


Eq. (4.22) should hold for all the models given, so the number of imposed constraints is m — 12. 
The determination of the weighting coefficient matrices W and о in (4.12 and 4.13) was done by the „шта! and error" method. 


The best results obtained are as follows : 


| 000.0 for icorresponding to ф = 0 


i о, =1 for 1-1,2. 20 (4.23) 
1.0 — otherwise 


Results of the approximation of cross curves we) for the S60 models (Tab.4.1) is shown in Fig.4.2. 
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a) b) c) 
Wel 8 А 8 Wp ll 8 
12 = 1 | то 1 
Св=0.60 | Em 1 Св=0.80 1 
0.8 
9.8 ОЧЕ 
12 – hz1.610 (о/в) 12 — Ь=1.610 (o/a) 1,2 – ћ= 1.610 (о/а) 
34 – h- 1.786 (o/a) 34 — h- 1.726 (0/8) 34 — h- 1.726 (o/a) 
0.4 5,6 – h- 1.B13 (о/а) 5,8 — ћ=1.813 (о/а) 94 5,6 - h- 1.813 (о/в) 
7,8 – h-1.900 (o/a) 7,8 — ћ=1.900 (o/a) 7,8 — ћ=1.900 (о/а) 
92 0.2 
ф гаа] ф гад] Фф trad 


De 0.4 9.6 9.8 1.0 1.8 1.4 157 ов 0.4 0.6 0.8 10 T% 14 1.67 Oe 9.4 9.8 0.8 1.0 LE 140 ал? 


Fig. 4.2. Comparison of the original (о) and approximated (a) cross curves of S60 


4.3. Correction of the cross curves function Wh for small H/T ratios 


The function wb (ф,Сь ,h) approximating the cross curves of the S60 data is formally valid for the form parameters from the 
region : (Cy, h) e [0.60,0.80] x [1.61,1.90], b = b? = 2.5. Such a range of data evidently corresponds to the standard hull form 
parameters of sea-going ships but does not cover the regions typical of the sea-river vessels where the block coefficient Св 
reaches values as large as 0.90 and h as small as 1.1 (Tab.4.3). 


Tab. 4.3. Hull form of a sea-river vessel [18]. a - body lines, b - a sample of statistics [13] 


It turned out that a natural extrapolation of the cross curves function (4.19) obtained for S60 towards large Св and small h-H/T 
is acceptable as far as Св and not acceptable as far as h is concerned. Fig.4.3 illustrates the problem. So a special extrapolation 
was necessary for the function w$ to be used in an optimisation model of sea-river ships. As a result, a corrected cross curves 2D 
function wy,” (ф, h) has been defined (by an interpolation technique), based on the boundary data of S60 (for h = 1.61) and 
a cross curve data for the SINE-205 sea-river vessel [18] (for h = 1.239) in the following tensor product form : 


wg ($.h) 2 d' (9): g(b) =e" (0): D-g(h) (424) 
where : d(@) — interpolatory characteristics, D - interpolatory matrix, е(ф), о(ф) — bases, dim(d) = dim(g) = 3, dim(e) = 5. 
For details - see [13]. 


1. S60 - H/T = 1.60 " 3. 560 – H/T = 1.20. % | ман 


ВЛ 22.500 
1 -S60 -Н/Т=1.239 | 
2- SIN -HIT=1,239 
3 - S60 -HIT-1.615 
—— - 
il фиеді 
10 20 30 ас ES eu т ac LJ 


| 10 ГЕ 3c 40 so 60 70 LJ ж 


Fig. 4.3. А natural extrapolation of the cross curves function who of S60 towards the region of small H/T ratio 


After correction, a definition of the basic (B/T = b= 2.5) cross curves function is : 


0 
wR(,Ca,h) for hz1.610 
wa AE n ind (4.25) 
erch, h) fo ћ<1.610 
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Fig.4.4 show comparison of the basic cross curves functions before and after correction with regard to the small H/T ratio. 


А. Before correction 


CB = 0.870 
H/T = 1.224 


CB = 0.870 
B/T = 2.590 


CB = 0 870 
НТ = 1.224 


Яғ, 721790 


hired 


Fig. 4.4. Basic (b = b°=2.5) cross curves function before and after correction with regard to the small H/T ratio 


4.4. Affine transformation of the cross curves function w},(-) with regard to B/T 


The definition of cross curves function wal ), obtained as the result of approximation of the Series 60 data and then corrected 
for small H/T ratio is only valid for B/T =b = 2.5. Accounting for a new B/T ratio = b! # b’ can be done by the way of affine 
transformation of a hull form. The situation is illustrated in Fig.4.5. A system of hull coordinates is transformed when the hull is 
heeled. As a result, for equivalent waterlines (e g. CWL or deck immersion waterline), one obtains, among other things, 
a changed angle of heel (фо — 91), changed coordinates of B (Во — В) and, in consequence, the desired transformation of cross 
curves мо (ф,Сь,ћ ;Б=Ь®)— wh (фСв,ћђ). 


th 


ІШ A, =b'/b° 
Zo 


Fig. 4.5. Affine transformation of coordinate system 


In order to derive a particular form of this transformation, 
let us recall the inverse relation between the (y-z) and (w-t) systems (3.1b) : 


y cosó —5іпф | |w 


=. (4.26) 
Z sind  cosó| t 
Using the formula (4.26) for the points Bo and B, one has accordingly 
the following (w-t) -> (y-z) transformation of coordinates : 
0 т 0 
ув |_ соѕфо —sinóg | Wn (4.278) 
0 | Lu 0 | 
Zh SINdy coso | | th 
wh = у ‘COS, + 25 “Sind, = (Àp -y5):cosó, +20, -sing (4.27b) 
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After substitution of (4.27a) into (4.27b) one arrives at the formula : 


Wp = cos, cosh; [we ‘(Ay + фо tgo) + th - (А tgo +1001)] (4.28) 


Further algebraic manipulations lead to the final transformation function (the indices | were dropped) : 


кф М] | w? 
еы =| EN viles] (4.29) 
БА) | | бо] 
where the auxiliary functions involved in (4.29) are : 
fy (ФА) = Ap TCD, Àp) (a) 
5(0,,) = 05-4-24)-бі29-т(фҚ,2,) (b) 
1 
i) = анан (с) 
| cos? +A, sin? $ (4.30) 
$o(0) 
5%) = zp- |в) (d) 
0 
$09(6.4,) = arctg(A, :0ф) (e) 


Fig. 4.6 presents graphs of a heeling angle function (4.30e) for different A, е [0.25 , 2.00] 
and the functions wg(-) (4.29) and ibo (4.30 d). 


10 20 90 40 50 66 70 80 56 


Fig. 4.6. Affine transformation of cross curves with regard to the B/T ratio (Аъ = b/b?) 


4.5. Analytical definition of the righting lever curve Ip(-) 


An analytical definition of the righting lever function for the hull form parameters о = (Cp, h , b) : 


Іф. g:zG) = Ув(ф, g) - za: sino (4.31) 
can instantly be obtained based on the analytical definition of the cross curves function (4.29) : 
Wp (>, в) = ws (0. в) = WR (O, Cg. h) fi, (,b) +0500, Съ.) f, (9, b) (4.32) 


When looking at the formulae above it becomes obvious that the righting lever function so defined would be too complicated 
for performing various analytical operations (with regard to ф) when generating the stability constraints. In order to cope with 
this problem, a 4D function 1; (0, 8; Zg) has been approximated by an 1D function k (0;g,z,) (for g , zg = const.) using the 
known solution to the approximation problem (4.18). The idea is based on the observation that the matrix B depends only on the 
given nodes (matrices E, C), the types of constraints (matrix C), the weighting factors (matrices УУ, @) and does not depend on 
ordinates of the approximated function or values of constraints (c). It follows from the above that for all the mentioned elements 
kept fixed, the matrix B and its inversion В can be generated once and be used many times for different combinations of the g 
and zg parameters, when running the optimisation procedure. Bearing in mind the need for simplicity, a linear polynomial 
representation has been assumed for the approximating function : 


n 6 
md s 1—1 
Ё(ща) =а' -е(и)= Уо еи) ©  l&(6:9) 2 Уо: ф (4.33) 
1=1 i=l 2» 
where: а = (ae Sener - is a solution of the approximation problem and the definition of 1, (-). Note, that о = Q(g, Zg). 
In analogy to approximation of cross curves, the only constraint (m — 1) imposed on the lever righting function is : 


d ~ 
ag 8X9 = 0 = ав(Св) ro ые (4.34) 
The resulting righting lever function has the desired simple form : 
e= БОСНА а Бел) Уа а Билај“ (4.35) 


i-l 
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Fig. 4.7 presents some diagrams of the righting lever curve as a 4D function Ір(ф, Cg , h , b ; zg) according to the definition (4.31). 


9 


7 T - b 
0.0 $ia 10 1571 


Fig. 4.7. Righting lever curve as a multivariable function of $ and hull form parameters of a ship 


5. STABILITY CONSTRAINTS 
IN POST-OPTIMISATION STUDIES 


An example of post-optimisation stability-oriented study 
is presented on the background of the optimum solution of the 
fleet of ships optimisation problem. To perform the calcula- 
tions, the optimisation model of a fleet has been implemented 
in the EUROS computer program [13]. 


5.1. The data 


A typical sea-river vessel and a sample 
of the main fleet transportation data are shown in Table 5.1. 


gramming method of optimisation, a combined, double level, 
algorithm for constrained problems has been applied, inclu- 
ding a penalty function shifting method by Powell and Wierz- 
bicki, for the constrained part, and Rosebrock's method for the 
unconstrained part of the problem ([3],[15]). 


5.3. Numerical results 


Numerical examples of optimisation are presented in Ta- 
ble 5.2. The calculations have been carried out separately for 
three criteria of stability: IMO, PRS, and HSMB with the same 
data, the same objective function: (NPV) and the same star- 
ting point x . The optimisation process has been repeated twi- 


Tab. 5.1. Fleet transportation problem : 
Typical sea-river vessel for short shipping in the North and Baltic Seas [18] and main transportation data [13] 


Sea-river vessel 


Sea-river ships, as multi-purpose vessels, carry different 
kind of commodities, such as containers, general cargo and 
bulk cargo. In the transportation data, these are accounted for 
by fixing relative quotas of cargo types for the whole period of 
fleet operation. 

Referring to the terminology introduced in Table 2.2a, the 
overall quantities of numerical data used in the model (inclu- 
ding those in Tab.5.1), in the four categories of constants, are: 
VESSEL - 154, CARGO 10, ENVIRONMENT - 374, OPE- 
RATION - 326. This sums up to 864 constants in total. 


5.2. The model and NPM algorithms 


The implemented model for fleet optimisation problem 
contains 12 decision variables (Tab.5.2), 70 constraints in total 
(including 26 stability constraints) and 3 alternative objective 
functions: Net Present Value (NPV), Internal Rate of Return 
(IRR) and Required Freight Rate (RFR). As a non-linear pro- 
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A sample of main transportation data : 


> | Period of operation : М = 15 years 
No. of ports : OUT / INP / total = 10/11/20 


No. of river ports = 12 

Mass of cargo OUT - 110 000.0 tons / year 

Mass of cargo IN = 103 000.0 tons / year 
Distance : OUT / IN / total: 2827 / 2959 / 5786 [NM] 
Max. draught / velocity in river : 2.8 m/ 14 km/h 
Period of calling at ports : 6 - 14 days 

General cargo stowage factor : 1.0 - 2.0 [m°/t] 
Relative quotas of cargo type : 

асс / дос /qbc OUT = 0.8 / 0.2 / 0.0 [-] 

асс /дос / qbc IN =0.9/0.1/0.0 [-] 

qcc - containerised c., qgc - general c., qbc - bulk c. 


ce : the results ofthe first optimisation were used as a starting 
point to the second optimisation and then convergence was 
reached. It turned out that, from the viewpoint of the final va- 
lue of the objective function, such an approach improved the 
results remarkably. Apart from the NPV, the results in Table 
5.2 show the corresponding values of two remaining measu- 
res of merit: IRR and RFR. It can be seen from Table 5.2 that 
the results for the stability criteria by IMO and HSMB are 
identical. This will become obvious when one analyses the 
status of the stability constraints with regard to the feasible 
solutions region Q. 


5.4. Feasibility analysis - the case study 


Feasibility analysis is an element of post-optimisation stu- 
dies (together with parametric study and sensitivity analysis) 
that aims at investigating the status of a certain group of con- 
straints (in the case - the stability constraints) versus boundary 


Stability criteria as constraints in a fleet of ships optimisation problem 


Tab. 5.2. Results of optimisation of a fleet of multipurpose 'sea-river ' ships for three intact stability criteria. Objective function: NPV 


| x | Description — | Unit | Stat | МО | PRS | HSMB | 
| x |  Numberofshipinaflet — | [J | 4000 | 3.927 | 443 | 3927 | 
| 


: 


| Мо of container rows under deck — | [] | 3.00 | 445 | 3000 | 455 | 
Мо of container tiers under deck — |  [:] | 200 | 2302 | 229 | 2302 | 


Tab. 5.3. List of stability constraints 


г— 


3 
Ха 
Х5 
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Xg 
Xo 
X10 
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flt 
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-- 


[-3.[-— GM0atrivr — | мил | m | ___н рр 
(3| 5200 — [mi | јој | __ нмв __| — 40 | 60 | 
| 6 | — — GZrage — — [ min | [lg] | (— PRS | — 50 |) 46 | 
| 8 |  — Weathercriterion — | min | mad | мо | .— 3 | 06. | 
| 9 |  Turningheelangle — | max | [deg] | PRS |5 | — - | 
| 10 | — Statcheelange [| max | [dg] | мо | s“ | 6 | 
| 12 |___Агса џпдег GZ [0, 40] | min |[mrad]| —IMO.HSMB | _ 56 | 068 | 
| 13 |  AreaunderGZ [30. 40] | min ||ша)| —IMO.HSMB | 5 | 60. | 


constraints on the background of the remaining constraints, 
significant in a vicinity of the optimum solution. It can also be 
thought of as a preparatory, qualitative part of sensitivity ana- 
lysis. Table 5.3 presents a list of stability constraints together 
with their indices in the optimisation model implemented in 
the EUROS program. 

Sea-river ships, as multipurpose craft, operate in their life 
as different types of vessels, both from the viewpoint of their 
functionality and the corresponding stability regulations. To 
cope with the problem three extreme situations have been re- 
cognised in the model as crucial ones : 


(1) A ship operates as a ,,pure" container sea-going vessel. 
She then carries containers both in holds and on deck and 
is allowed to take water ballast into the double bottom and 
wing tanks. 


(ii) A ship operates as a ,,pure" general cargo sea-going vessel. 
As such, she does not load cargo on deck but, on the other 
hand, she is not allowed to take water ballast into her tanks. 


(111) A ship operates as а ,,pure” container, river vessel. She 
then carriers containers both in holds and on deck but, as 
the total mass of ship must be smaller than that at sea due 
to restricted draught, the ship is not allowed to take a wa- 
ter ballast. 


All the three situations have been accounted for by formu- 
lation ofa triple set of corresponding stability constraints „ас- 
ting" simultaneously. Such an approach assures that stability 
requirements will also be satisfied in all the intermediate lo- 
ading conditions which can occur in operation. A feasibility 
analysis has been done in a graphical form using the method 
and software worked out in [9]. The results are presented in 
Tab. 5.4. The constraints are presented via visualisation of the- 
ir contour lines in the vicinity of the optimum solution. The 
arguments ofthe graphs are ship's dimensions: B vs. H (1) and 
hull form ratios: B/T vs. H/T (2), the most essential ship para- 
meters as far as the intact stability 1s concerned. The descrip- 
tion in Tables 5.4 refers to the terminology introduced in Chapter 
2. The stability constraints recognise three potential functional 
ship categories : a sea-going container carrier, a sea-going ge- 
neral-cargo ship and a container river ship. The feasible solu- 
tion region Q is shown as а darkened area and the ,,tufts” on 
the contour lines point in the unfeasible direction. The parame- 
ters ofa fleet, which are not arguments of the graphs, are those 
in Table 5.2. The feasibility analysis is summarised in Table 
5.5 where the attributes of stability constraints refer to those in 
Tab. 2.3. 

It follows from Tab.5.5 that in the case under investigation 
the most restrictive stability criteria are those by PRS, where 
the criterion No. 50 (GZ range) appears as a boundary and 
active constraint. 
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Tab. 5.4. Feasibility analysis - the case study 
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Legend : a - active constraint, p - passive constraint 


cs - container ship at sea, gs - general cargo ship at sea, cr - container ship in river 


Annual mass of cargo (INP) 
Mass reserve factor 
Mass reserve factor 


Boundary constraints : 


(min) 
(min) 
(max) 
(min) 


POLISH MARITIME RESEARCH, Special issue, 2004 


Vertical container 
Freeboard 

Height of bow 
L/B ratio 


clearance [m] 


[m] 


Stability criteria as constraints in a fleet of ships optimisation problem 


Tab. 5.5. The status of stability constraints versus boundary constraints 


B/T vs. H/T 


6. SUMMARY AND CONCLUSIONS 


A post-optimisation study — stability-oriented feasibility 
analysis has been demonstrated on the background of a marine 
engineering design problem - optimisation of a fleet of sea- 
river ships for a regular shipping line in the area of the North 
and Baltic Seas. Two questions can be addressed based on the 
results obtained: (1) why the full intact stability criteria (such as 
IMO or similar) should be included in the optimisation model, 
and (ii) how to prepare such a model to accomplish the task. 


As far as the first question 1s concerned, 
the following justifications seem to be of importance : 


designers and potential customers (owners) express speci- 
fic and growing interest in the safety of ships 

the stability criteria do exist as formal and legal design re- 
strictions so they can not be neglected 

their presence in the model improves its quality, thereby 
making it more credible 

their incorporating into the model enables post-optimisa- 
tion studies to be undertaken, such as parametric study, fe- 
asibility analysis and their quantitative extension - sensiti- 
vity analysis, with special emphasis on the economic aspects 
of the solution (,,shadow prices” and risk assessment). 


У У У У 


As to the second question, an attempt has been made to 
define arbitrary stability criteria as constraints based on a com- 
plete analytical definition of all the necessary geometrical cha- 
racteristics of hull form. A method has been proposed to use 
a systematic series of ship body forms. As an example, appro- 
ximation of cross curves and accompanying characteristics of 
Series 60 has been demonstrated which makes it possible to 
determine a righting lever curve as a multi-variable function of 
the heeling angle and ship parameters — decision variables in 
the fleet optimisation model. 


Fora particular case under study — a fleet of sea-river ships, 
Series 60 turned out to be not suitable for the whole range of 
form parameters expected of such type of ships (very large Cp 
and very small H/T). Consequently some corrections have been 
made to fit the available characteristics of SINE-205 sea-river 
vessel. Such an approach, however, must be regarded as a tem- 
porary solution. A research project is under way to work out 
a stability-oriented, analytically defined series of hull forms 
for specific functional types of ships such as full cellular con- 
tainer ships, tankers and so on, to be used in the corresponding 
fleet optimisation models. 


NOMENCLATURE 

B/b - breadth of ship / B/T ratio 

е/С - vector / set of constants 

Св - block coefficient 

50) - function of a j constraint 

H/h - depth of ship / H/T ratio 

Jg / Ја / Jg - sets of indices of boundary / stability / remaining 
constraints 

Lu - heeling lever curve 


lk / Lg - non-dimensional / dimensional [m] righting lever 
curve 

- number of constraints 

- number of decision variables 

- objective function 


R: - n-dimensional arithmetic space 

To - non-dimensional transversal metacentric radius 
T - design draught of a ship 

Wp - cross curves function 

x - decision variable vector 

ZG - non-dimensional centre of gravity of a ship 

ф - heeling angle 

Q - feasible solution region 
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Hull strength assessment of the designed ships 
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ABSTRACT 


The paper presents the method and scope of analysis of the hull strength of four ships designed in the 

Eureka project. Criteria of the Det Norske Veritas rules for classification and construction of ships were 

used in the analysis. Specific features of the structure of each hull have been described and the resulting 
problems connected with ensuring adequate strength. 


Keywords : hull strength, FEM calculations 


INTRODUCTION 


The strength analysis was performed 
for the following four ships : 


product tanker (length overall Год = 138.10 m, breadth 
B = 22.50 m, depth D = 12.80 m, draught T = 8.70 m) 
ro-ro ship (Год = 156.72 m, B = 24.80 m, D = 19.60 m, 
T = 6.50 m) 

,Tiver-sea" ship (Год = 89.45 m, B = 11.40 m, D = 5.45 m, 
T = 4.40 m — at sea, T = 2.80m – in river) 

container carrier (Год = 138.10 m, В = 22.50 m, D = 11.20 m, 
Т = 8.55 m). 


The design requirement of those ships was their ecological 
cleanness, i.e. а CLEAN DESIGN symbol given in the Det 
Norske Veritas class description [1]. 

The ecological cleanness of a ship is first of all connected 
with its installations and equipment. It does not impose any 
specific requirements on the internal subdivision or structure 
of the hull. A sufficient precaution 15 the use of classic solu- 
tions, 1.e. double bottom and double side in the case of a pro- 
duct tanker. However, double bottom and double side have been 
used in all the designed ships as they are classic solutions for 
those ship types. They protect the sea against pollution in the 
emergency situations, e.g. grounding, collision with another 
ship etc. The ecological cleanness feature does not mean, how- 
ever, that the hull strength standards have to be raised. 


У У У У 


The structure strength assessment was performed 
in accordance with the rule requirements [1, 2, 5]. 


The following structure strength levels had to be checked : 


+ overall strength 
+ zone strength 
+ local strength. 


Additionally, the criteria of total normal stresses in the lon- 
gitudinal girders and in longitudinal stiffeners of side shell, 
decks etc. from the general, zone and local bending have to be 
fulfilled. 

The overall strength involves the longitudinal bending of 
the hull in the vertical plane. In the case of ships with wide 
hatch openings, also torsion and bending in the horizontal pla- 
ne should be analysed. The calculations are relatively simple 
as the rules require a bent beam model to be used in the analy- 
sis of bending in the vertical or horizontal plane. In the torsion 
analysis, the use of constrained torsion theory of thin-walled 
bars is admissible. However, it is recommended to use a shell- 
-bar FEM (Finite Element Method) model of the whole hull 
and then the calculations are time-consuming. 

The zone strength involves bending of the hull girder sys- 
tem. Calculations are often time-consuming as in general it is 
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required to use a shell-bar FEM model. The rules [2] generally 
allow to use also simpler models — space frames, grids, flat 
frames, continuous beams. Such calculations are still relative- 
ly time-consuming. 

The simplest are local strength calculations. The respective 
formulae for plate thickness or stiffener section modulus are 
given in the rules. 

In the above mentioned analyses the stress level is chec- 
ked, by the permissible stress criterion. The structure should 
also fulfil the criteria given in the rules [1] for the structure 
element stability and fatigue life. 


Specific hull strength problems of individual ships 
are described below. 


PRODUCT TANKER 


The most labour consuming problem of the hull strength 
assessment was the zone strength analysis. 

The cargo section of the hull is divided with transverse bulk- 
heads into 6 compartments. Compartments no.2 to 6 (counting 
from the bow) are additionally divided with a longitudinal bulk- 
head in the hull plane of symmetry. 

An FEM model in the form of space frame was created 
comprising a hull fragment of the length of + 1 + V? cargo 
compartment (Fig.1). The possibility of asymmetric liquid cargo 
distribution made it necessary to develop a full hull breadth 
FEM model. 


lox i^ w е 1 
АРЕ 
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Fig. 1. FEM model of a product tanker hull module frame 


Calculations revealed a high level of the reduced stresses 
in floors and transverse web frames coming from shear in the 
area of lightening and communication holes. It was necessary 
to correct the original version of the structure. The holes were 
made smaller and thicker plates were used. Also the original 
transverse corrugated bulkheads appeared not strong enough. 
They had to be made thicker. 

Difficulties were met with ensuring local strength in the 
specific product tanker conditions of dynamic loads in partial- 
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ly filled tanks (sloshing). The problem was connected with the 
deck structure and upper parts of bulkheads in the first fore 
cargo tank. It was necessary to correct the original plate thick- 
nesses and stiffener sizes as well as the main girders in that 
region. 

A characteristic feature 15 that the elements of a product 
tanker hull structure (total length of 138.1 m) with double ship 
side shell and double bottom, without hatch openings in the 
deck, dimensioned according to the local and zone strength cri- 
teria, ensure automatic fulfilment of the general longitudinal 
bending criteria. 


RO-RO SHIP 


Characteristic solutions of the ship are the following : 


& three loading levels (inner bottom, main deck, upper deck 
with the loading ramp hatch openings) 

no transverse bulkheads in the cargo section of the hull 
double bottom, double side from bottom to main deck and 
single side above. 


оо 


Assessment of the zone strength of such structure was рег- 
formed by means of simple FEM models in the form of flat 
frames and grids loaded with the container or trailer weight. 
Additionally a shell FEM model was created of a typical web 
frame (Fig.2) in order to verify the structure correctness in the 
bottom and side joint region with communication openings. 
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Fig. 2. Transverse web frame shell FEM model 


Calculations indicated the necessity of strengthening that 
structure region in relation to the original version. Also the 
upper deck web beams in the side region appeared too weak. 
Their section modulus had to be increased. 

In the case of a ro-ro ship, the deck shell and stiffener di- 
mensioning load is that of the vehicle wheel load. The use of 
80 t and 60 t roll-trailers and 30.5 t carrying capacity container 
stackers was assumed. That gives relatively great wheel load 
values. Increase of the originally proposed cargo deck shell 
thickness and beam strength appeared necessary. 

As in the case of the product tanker, the structure dimen- 
sioned according to the local and zone strength criteria auto- 
matically fulfilled the general overall strength criteria. Impor- 
tant 1s here the relatively large depth of the ship and lack of 
hatch openings in the middle part of the hull. 


RIVER-SEA SHIP 


This is the smallest 
of the four designed ships (Loa = 89.45 m). 


Characteristic features of the hull are: wide hatch opening 
as long as 62 m over the only hold of the ship, strong double 


bottom adjusted to relatively large local loads from cargo, do- 
uble ship side and high hatch coamings on the main (single) 


deck. | Ши 
The structure dimensioning loads 


are those during the operation at sea. 


It is interesting that the overall strength problems occurred 
in the case of a small ship. Normal stresses in the upper fibres 
of hatch coamings reached, in the general bending conditions, 
values near the admissible level. 

In the case of that ship type, important is the problem of 
hull torsion on a course oblique to sea waves. Torsion analysis 
was performed with the use of thin-walled bar constrained tor- 
sion theory [5]. Torsional displacements of the hull appeared 
small but stresses reached significant values. 

The greatest normal torsional stresses occurred in the hatch 
rear end area, in the upper edge of hatch longitudinal coamings 
(75 MPa). Total normal stresses in that place, from torsion, 
hull bending in the vertical and horizontal plane and from deck 
strake bending in the horizontal plane (due to bottom and side 
loads) reached the admissible stress level, i.e. 195 MPa. 

In the zone strength analysis, it was necessary to use a FEM 
model in the form ofa space frame, containing main structural 
elements in the whole cargo part of the hull and across the 
whole breadth of the hull, in order to analyse e.g. the structure 
stresses in the conditions of ship heel with containers on deck 


(Fig.3). 


Fig. 3. FEM model of a system 
of structure elements in the river-sea ship hold 


With the use of such computational model a relatively gre- 
at deformability of the hull was demonstrated. Convergence of 
the ship sides in certain loading conditions was estimated at 
65 mm, which may require special hatch covers to be used 
acting as hatchway beams. 


138 


Fig. 4. Shell FEM model of the bilge region of river-sea ship hull 


Computations showed also exceeded admissible stresses 
in the bilge region floors. Therefore, a more precise (shell) FEM 
model of that region was prepared (Fig.4). Calculations have 
shown that holes are inadmissible in that region. 

Standard local strength calculations were also performed 
in the hull strength analysis. The design requirement that the 
ship should be able to transport iron ore, with its approximate- 
ly 100 kPa static load on the inner bottom, resulted in a relati- 
vely massive structure of that hull fragment. 
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CONTAINER CARRIER 


A characteristic feature of that ship hull 
are wide hatch openings and narrow double sides. 


There are 5 hatch openings in the deck, connected with 
batten plates. Such hull configuration causes the cross section 
neutral axis to be positioned much closer to the base plane than 
to the deck. 

There were some problems with adequate general bending 
strength in the vertical plane. In the case of a small ship (Loa 
= 140,0 m), it appeared effective to provide high tensile steel 
(315 MPa yield point) girders in the middle part of the hull in 
the deck region. 

The rules [1] require that torsion analysis be performed for 
ship hulls ofthe above mentioned characteristics. The calcula- 
tions were carried out with a twisted bar model. Bending of 
a frame created by the side deck strakes and the batten plates 
(webs of the frame beams) was taken into account. The beam 
flanges are hatch coamings together with the internal and exter- 
nal side fragments or transverse bulkheads below deck. The cal- 
culations were performed in accordance with the procedure 
required in [6]. 

Normal stresses from bending of the deck strakes due to 
deformation of the above mentioned frame (hull cross section 
torsional deplanation) reached the 30 MPa level. This is a si- 
gnificant value. 

The total normal stresses from hull bending in the vertical 
and horizontal plane, torsion (due to constrainment of the de- 
planation), deck bending (as described above) and deck ben- 
ding from sea water pressure on ship sides, reached a 231.5 
MPa level, which is as much as 9596 of the admissible value. 
Those stresses occurred on the top of the no. 5 hatch rear end 
coaming. 

Container carrier torsion is therefore 
a significant problem, even in the case of a small ship. 


The zone strength analysis was performed with a space fra- 
me model of the structural elements of halves of two adjacent 
holds (Fig.5). 


Fig. 5. FEM model of the container carrier structure elements 


Also an FEM shell model of a single web frame was deve- 
loped (Fig.6) in order to check the stress level in the web com- 
munication opening areas. 

Standard local strength calculations were carried out in the 
same way as for the other ships. In the case of container car- 
rier, it was necessary to check the strength of the bow section 
bottom structure subjected to slamming loads. This is connec- 
ted with the relatively high service speed of those ships. 


The calculations showed that the bottom plating 
in the original version of the ship design was too thin. 
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Fig. 6. FEM shell model of a container carrier frame 


FINAL REMARKS AND CONCLUSIONS 


** Four different ships have been designed in the EUREKA 
project. They are intended for different service and there- 
fore their hulls are subject to specific loads. 

* In the case of a product tanker, this is the sloshing pressure 
having a significant impact on the shell plate thickness as 
well as on the size of stiffeners and structural elements of 
ship sides, bulkheads and deck, particularly in the fore part 
of the ship. 

* In the case of a ro-ro ship, characteristic are the vehicle 
wheel loads on decks. 

* In the case ofa relatively small river-sea ship, problems oc- 
curred with providing adequate general strength in the ben- 
ding and torsion condition, with bottom strength under consi- 
derable local loads from cargo, and with the hull stiffness. 

** The container carrier hull is significantly strained in the 
overall bending and torsion condition. In this case, the use 
of high tensile steel for girders in the upper part of the hull 
appeared effective. 

** A common feature of all the hulls is double bottom and 
double sides — at least below the main deck. 

* The ship hull design process proceeded in a rational way. 
The preliminary scantlings proposed by the designers were 
corrected ,,upwards” as a result of the performed strength 
analyses. 
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ABSTRACT 


Importance of reliability and safety in the operation phase of engine rooms has been pro- 
ven. A concept of designing engine rooms taking into consideration their reliability, safety 
of functioning as well as ecological features has been proposed. The need for as well as the 
main principles of carrying out empirical research in particular ship design stages have 
been introduced. Possibilities of the use of the semi-Markov theory in designing engine 
rooms have been justified. Possibilities of formulation and importance of hypothesis in 
scientific research related to engine rooms of sea going and inland ships depending on 


their particular specific features have been introduced. Examples of such hypothesis and their verification 
methods have been given. Possibilities of the use of semi-Markov processes in optimization of the ship 
operation have been signalled. 


Keywords : ship power plant, sea-going ship, reliability and safety, design, semi-Markovian process 


INTRODUCTION 


Ship engine rooms should be designed so that, assuming 
that the ship carries her tasks, she could be operated with in the 
highest in these conditions (if possible — optimum) rentability 
but keeping proper reliability, required safety of motion as well 
as meeting ecological measures. This is especially important 
for transport ships operating in relatively restricted sea areas 
such as the Baltic Sea. Thus problems of economic functio- 
ning of engine rooms of ships taking into consideration their 
reliability should be judged as the most important. 

Design solutions used in engine rooms can be different but 
the general way of forming their reliability (understood as the 
choice of types of their activities as well as the choice of possi- 
ble resources ensuring reaching needed level of reliability) can 
in principle be the same. Building the reliability of each ship 
engine room in the design phase is based on consecutive solu- 
tion of problems resulting from the need to obtain answers to 
the following questions : 


O What criteria of creation of reliability should be adopted ? 


O What is the risk of excessive increase of the cost of desig- 
ning in case of starting an original design of engine room ? 


O What could be the cost of obtaining the required level of 
reliability of an engine room similar to an already designed 
and operated engine room ? 


O How to divide the machines and equipment of the engine 
room in order to shorten the design phase and in the same 
time to obtain the required reliability ? 


O Against which requirements should the reliability be con- 
sidered ? 


This is why scientific research is indispensable in the de- 
sign stage of ship engine rooms. Such research is also indi- 
spensable in further stages of existing of an engine room, na- 
mely during its building period as well as operation. Such scien- 
tific research enables to create scientific knowledge, i.e. such 


a knowledge which is justified using scientific methods of its 
creation. In result of such research information of cognitive as 
well as useful value can be obtained. In the first case new know- 
ledge is created in process of verification of hypothesis and/or 
proving thesis explaining relations between new facts, ascerta- 
ined during research carried out during operation of engine 
rooms of particular ships, used by different shipwoners, with 
the existing scientific knowledge. The second case refers to 
the sphere of scientific knowledge which is used in practice as 
to designing engine rooms as well as other structural nods e.g. 
ship hulls. This way the number of true information creating 
useful knowledge but checked using scientific methods is in- 
creased and in the same time the number of so called “com- 
mon sense knowledge" information is decreased. This way 
useful knowledge becomes a more precise tool which can be 
used in the design processes and later on in the phases of buil- 
ding and operation of engine room [4,8,9,17,18 and 19]. 

The importance of scientific research in designing engine 
rooms is based on the fact that their introduction enables to 
explain facts ascertained during designing, building and ope- 
ration of already existing engine rooms and using desirable 
solutions in the phase of consecutive designing. The importan- 
ce of such research is also based on the fact that during repair 
of engine rooms they should be modernized using knowledge 
of designers who followed up the construction as well as the 
operation of the engine rooms. 


To show this importance one should first signal 
the contemporary design concept of ships and navy vessels. 


OUTLINE OF THE DESIGN CONCEPT 
OF MODERN SHIPS 
VERSUS SCIENTIFIC RESEARCH 


Deterministic design methods are used in contemporary 
design of engine rooms (as well as of ships treated as a whole) 
[1,8,15]. Although the use ofthese methods enables to prepare 
designs leading towards building engine rooms with a level of 
reliability judged intuitively as high (in principle the measure 
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of such assessment is subjective or psychological probability) 
but does not allow to present such reliability in form of reliabi- 
lity factors. Such indicators even in form of probabilities (logi- 
cal or statistic) are indispensable for rational planning and later 
on control of the process of operation of engine room. This is 
why there is a need to develop probabilistic methods shaping 
the reliability of engine rooms and to use them in the engine 
room design phase [6,7,17,18 and 19]. 

Probabilistic methods enabling shaping reliability and func- 
tional safety including proper ecological safety can be develo- 
ped using the contemporary theory of reliability and safety of 
complex technical systems (engine rooms obviously can be 
treated as such systems), probability calculation methods, ma- 
thematical statistics, the theory of machines (mainly thermal 
machines) installed in engine rooms as well as using the theory 
of semi-Markov processes. In turn in order to prepare a set of 
reliability conditions of engine rooms one should use technical 
diagnostics as well as the results of damages of machines and 
equipment used in similar engine rooms. The knowledge of 
reliability and safety indicators of particular engine rooms ena- 
bles to take rational decisions related to operation. 

During operation of sea going ships decisions are taken al- 
ways in stochastic decision situation i.e. in conditions of uncer- 
tainty (statistic risk conditions). This means that the probabili- 
ty and induction (mathematic) statistic rules are to be used when 
making decision. Decisions concerning operation are taken be- 
fore commencement of the operation as well as during opera- 
tion of the mentioned ships. This means that at least once such 
decisions are taken basing on initial information (obtained for 
example from research of reliability of engine rooms as well as 
of their particular machinery and equipment or basing on data- 
banks of similar engine rooms). Such decisions can be named 
“a priori” decisions. Following decisions are based on infor- 
mation obtained during operation ofthe technical systems (e.g. 
in result of diagnostics, not only of technical diagnostics) and 
can be described as “a posteriori" information. 

Decisions taken at the beginning ofthe operation are indis- 
pensable to plan the process of operation and maintenance of 
the said ship engine rooms. When taking decisions one should 
consider the statistic risk estimated as probability of taking 
wrong decision resulting from [2,7,9,19] : 


© impossibility of precise assessment of unknown parame- 
ters of the distribution of random variables, especially of 
such variables, which represent conditions of the operation 
process of the engine room and of particular machinery and 
equipment 


© lack of possibility to prepare full and/or reliable enough 
information indispensable to take the right decision. 


The first case generates stochastic type of errors being sub- 
ject of the so called statistic precision of conclusions while 
the second — both random mistakes as well as such which can- 
not be described as random (systematic). Determining the lat- 
ter errors is a subject, which I suggest to name as problems of 
precision or scrupulous way of concluding. In total determi- 
ning such errors is a question described as accuracy of stati- 
stical conclusions. Accuracy of conclusions results from the 
present level of scientific and applied knowledge, while scru- 
pulous way of concluding results from insufficient apprecia- 
tion or neglecting some information by the decision maker, 
who first could have checked that the information was really 
not so important. On the other hand decisions made during 
ship engine room operation are erroneous or irrational due to 
difficulties in preparation of full diagnosis and lack of proper 
credibility of information related to the technical conditions of 
these engine rooms, their machinery and equipment as well as 
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to the foreseen external conditions (atmospheric and sea con- 
ditions), which could be expected during the ship’s voyage 
[3,5,10,12]. 

In the presented decision situation taking rational deci- 
sions is possible in case of use of statistic theory of decision 
[4,5,18,19]. However determining the set of decisions, which 
can possibly be taken in accordance with the assumed crite- 
rions of optimization, requires identification of problems of 
assessment of reliability of the engine room, at least of its main 
energetic systems. 

Different reliability and safety indicators of engine rooms, 
its machinery and equipment can be considered in designing 
engine rooms. 


The most useful reliability indicators are [6,19] : 


probability of correct work until the first damage 
probability of correct work between two consecutive damages 
single-dimensional distribution of process (instantaneous 
distribution) the elements of which are P,(t)functions de- 
scribing the probability that in moment / the process would 
attain condition k 
* limit distribution of the process Р= lim P(Y(t) = j} 

tco 


+» + 


* conditional probability described as probability of transi- 
tion of the process from i to j condition, 

P; (t) = P(Y(t) = j/ Y (0) = i} (probability of transition) 

* distribution of the time of the first transition from i to sub- 
set of A [®;,(t)] conditions, if such subset is a one element 
set — to condition j , i.e. condition ®;;(t) 

* distribution Ф; (0) of the time of return of the process back 
to condition j 

Ж asymptotic distribution of the renewal process (У (1):120) 
generated by the distance in time of the return of the sto- 
chastic process (to condition j obtained from condition i , 
which in the moment ¢ takes the value equal to the number 
of "entries" of the process into condition j 

* approximate distribution of the total time of the presence 
of the process Y(t) in condition j provided that condition i 
18 the initial condition 

* expected value E(T;) ofthe time T; of i-duration ofthe pro- 

cess condition independently to which condition the transi- 

tion is being effectuated in moment Т, 

variance D*(T;) ofthe time T; ofthe duration of i-condition 

expected value E(Tj) of the time Ту of i-duration of the 

process condition provided that condition j will be next 
condition 

Ж expected value Е(Ө;) of the random variable Ө, descri- 
bing the time of return of the process to condition j 

Ж expected value ЕМ (0) of random variable V;;(t) descri- 

bing the number of "entries" of the process into condition j 

in the range [0, t] 

variance D? {V;(t)} of the random variable М, (6) 

average number of "entries" A,;(t) of the process into con- 

dition j occurring per time unit provided that the initial con- 

dition of the process is condition i (i.e. intensity of "en- 

tries” of the process into condition provided that Ү (0) = i) 

* limit intensity of "entries" of the process into condition / 
i.e. intensity j= limA4(t) : 

m 


жж 


жж 


Reaching numerical values mentioned in characteristics is 
possible only in case when the following conditions are met : 


ж Appropriate statistics have been gathered. The set includes 
values of assessed probabilities of transition ру, of the expec- 
ted value E(T;) etc. 


ж Construction of a stochastic model of the process of opera- 
tion of technical objects having small number of conditions 


Designing engine rooms of new generation ships 


and not much complicated number (in mathematical sense) 
of transitions (changes) from one condition into other con- 
dition. 


All the above mentioned types of reliability indicators are 
important for the operator but the most useful and in the same 
time the most easy to determine are the first two indicators. 


Among the most useful safety indicators one can mention : 


> expected value of random variable N,(t), i.e. the number of 
damages of the machinery and equipment of the engine 
room within the time range [0, t], causing a complicated 
situation describer with symbol E{N,(t)} 

> expected value of random variable N,(t), i.e. the number of 
damages of machinery and equipment of engine room in 
range [0, t], causing a dangerous situation described with 
symbol E{N,(t)} 

> expected value of random variable N,(t), i.e. the number of 
damages of machinery and equipment of engine room wi- 
thin the time range [0, t], causing a emergency situation 
described with symbol E{N,(t)} 

> expected value of random variable N;,(t), i.e. the number of 

damages of machinery and equipment of engine room wi- 

thin the time range [0, t], causing a catastrophic situation, 

described with symbol EfN,(t)] 

expected value of sail outs ofthe ship per one sea accident, 

described with symbol E МАО) 

expected value of sail outs of the ship per one sea catastro- 

phe, described with symbol E{N,(t)} 

expected value of death casualties per one sea catastrophe, 

described with symbol E{N,(t)} 

expected value of death casualties per one sea mile, descri- 

bed with symbol E{N,,,(t)} 

probability of safe ship (engine room) movement at sea P,(t), 

within the time range [0, t] 

probability of non-occurrence of a complicated situation of 

the engine room at sea P,(t), within the time range [0, t] 

probability of non-occurrence of a dangerous situation of 

the engine room at sea P,(t), within the time range [0, t] 

probability of non-occurrence of an emergency situation of 

the engine room at sea P,(t), within the time range [0, t] 

probability of non-occurrence of a catastrophic situation of 

the engine room at sea P,(t), within the time range [0, t]. 


У у у У У У У У У 


It is obvious that in order to estimate the mentioned expec- 
ted values of particular random variables №(1 = s, n, а, К) one 
should apply interval and not point estimation because only 
then the error of estimation can be determined. Determining 
the mentioned expected values is relatively simple. During in- 
vestigations within the time period [0, t] one can determine the 
average (arithmetic) value n; being the observed statistics Nj, 
which as it is well known has a normal asymptotic distribution 


С. 
N| m;—— 
vn 
where : 
mj - expected value (average) 
бі - average (standard) deviation of the random variable 


Nj, nj -number of registered events [1,4,16,18]. 


Reliability and safety requirements concerning functio- 
ning of the given engine room are as important as the techni- 
cal and economic requirements formulated by the shipowner 
ordering a new ship in order to obtain the highest possible 
rentability of the ship as well as efficiency of its propulsion. 
Considering reliability and safety factors of a ship's engine 
room at the stage of design enables to influence its readiness 


to start up in any moment of operation, which should be ma- 
intained using possibly lowest operation cost. Thus econo- 
mic calculations should be considered in the design stage in 
such a way as to control the process of engine room operation 
when in the same time maximizing the income. However this 
economic goal cannot be met without considering in the de- 
sign process all the conditions of construction as well as ope- 
ration and service of the engine room. This means that there 
are two staged which should be differentiated in the design of 
engine rooms [1,8,14] : 


* The first stage encompassing activities enabling designing 
of the process of fulfilling the need for production of ener- 
gy indispensable to ensure the movement of the vessel as 
well as social needs of the crew and passengers, proper or- 
ganization and control of work process proper for the de- 
sign of engine rooms. 


* The proper design stage i.e. the stage of creation of new 
values in result of which 15 created a more or less original 
design of the engine room 15 created. Said design stage is 
needed to realize the above mentioned process of meeting 
the need for energy production. 


The second design stage is finished with constructing, which 
18 leading to creation of the final form of the designed engine 
room, namely to the creation of the final ship structure. 

In both mentioned design stages one should consider such 
properties ofthe engine room as: functionality (functional cor- 
rectness), safety, efficiency, reliability, durability, easiness to 
make diagnosis, renewability, ability to be controlled, life expec- 
tancy, ergonomics, ecological aspects, low level of noise and 
vibrations and toxicity of exhaust gases. 

In the so understood design process one should assume that 
the basic information for design of engine room 15 the infor- 
mation which is related to the performed targets and condi- 
tions in which it can be realized. This means that one should 
first design the process ofthe engine room operation. The pro- 
cess will be a discreet process in conditions and continuous in 
time. It results from the hitherto considerations that a semi- 
-Markov process, especially of controlled semi-Markov can 
be a model for of such a process [3,4,5,6,7,17,19]. Therefore 
the theory of semi-Markov processes, can be used in the pro- 
cess of designing ship engine rooms. 


SEMI-MARKOV PROCESSES 
AS MODELS FOR REAL PROCESSES 
CONSIDERED IN THE PHASE 
OF ENGINE ROOM DESIGN 


The use of the theory of semi-Markov processes to create 
semi-Markov models of real processes as semi-Markov pro- 
cesses can take place only when [4,6,19] : 


* the Markov condition is met, what means that the evolu- 
tion of conditions of the investigated objects (e.g.engine 
rooms) in the future, for which the semi-Markov models 
have been built, should depend only on the present con- 
dition of the given object (from the condition in the pre- 
sent moment) and not from the functioning of the object 
in the "future", what means that the future of the object 
should not depend on its "past" and only on its "present 
condition". 

** random variables describing the time intervals when the 
subjects of investigations remain in particular conditions 
reveal different distributions than exponential. 


Thus when modelling, which has to lead to preparation of 
a semi-Markov model of the engine room operation of any 
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ship, one should consider the change of real conditions of the 
process of the given subject of investigations. 

In case of engine room of each vessel its operation pro- 
cess can be interpreted as a process of simultaneous changes of 
technical and service conditions [4,6,7,19]. Designing such 
process of engine room of any ship requires to solve using scien- 
tific methods the problem of forecast of particular conditions 
of the process. This requires to assume the following hypothe- 
sis (НІ): prognosis of the condition of the process of opera- 
tion of any ship or marine vessel in moment Ту + T, when it 
is known in the moment t, is known because its condition 
considered іп any moment t, (n—0,1,...,m; To € t1 <... <Tm) 
depends in an important way directly from the previous 
condition and not from the conditions which occurred ear- 
lier and time intervals when the condition occurred. 

It has to be noticed that the formulated hypothesis does not 
include any contradictions, which would classify it from the 
logical point of view even before its checking. 


The consequences of this hypothesis 
are as follows [4,6,19] : 


** probabilities (pi; i # j; 1] € М) of the transition of the pro- 
cess of operation of engine room or any of its machinery or 
equipment from any conditions i, in which the process is at 
the moment, to any other condition j do not depend from the 
fact in which conditions the process used to be in the past 


* intervals of the unconditional time of durations of particu- 
lar conditions of ship engine room operation process i are 
stochastically independent random variables (Т;; i€ М) 


* intervals ofthe duration of each possible occurrence of con- 
ditions i of the engine room process are random stochasti- 
cally independent variables (Т; i # j; ijj € М), provided 
that the next condition is one of the remaining processes j 
of the process. 


The mentioned consequences show the probabilistic law of 
changes of conditions ofthe mentioned operation process. They 
are not internally contradictory and their logic truth does not 
allow for any doubts. Thus the condition of no contradiction of 
consequences is met. Therefore there is no reason not to un- 
derstand the named consequences as one joint consequence K1 
and to use it to check empirically the introduced basic hypo- 
thesis H1. This means to verify it in order to accept it or classi- 
fy. The verification of the introduced hypothesis is based on 
experimental checking of the truth of the joint consequence 
KI. The verification of the introduced hypothesis H1 through 
experimental checking of the truth of the mentioned consequ- 
ence КІ requires to assume the truth of the following syntactic 
implication : 


HI KI (1) 


In such a case one can apply noninduction (induction) 
conclusion carried out according to the following scheme: 


(Kl, HI КІ) - Hl Q) 


The logical interpretation of the scheme (2) is as follows : 
ifthe experimental checking of the consequence K1 proved 
its rightness and if implication (1) is right, then hypothesis 
H1 is also right and can be accepted. Concluding made in 
accordance with scheme (2) is a means concluding of reduction 
type, which is one of the types of induction concluding and 
thus does not lead to sure conclusions and is only probable. 

In case when for any structural solution of ship's engine 
room, its machinery or equipment operated in any operation 
system it would result that empirical checking of the said con- 
sequence (K1) may contradict the formulated hypothesis (H1), 
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then the conclusion about the truth of hypothesis (H1) can take 
place in accordance with the "modus tollens" scheme : 


(~K1, HI > КІ) Р-НІ (3) 


The scheme of concluding defined 
by the relationship (3) has the following interpretation : 


if experiment does not confirm the correctness of the 
sequence K1, then in case of correctness of implication 
НІ = КІ, the hypothesis НІ is not correct and therefore 
cannot be accepted. 


Both ways of concluding (2) and (3) 
have to be supported with statistical conclusions [6,7,19]. 


In service practice it is important to plan the preventive 
maintenance of particular engine rooms as well as of particular 
machinery and equipment with regards to the undertaken tasks 
and the need for gathering appropriate resources to carry out 
indispensable operations. 

The application of semi-Markov processes as models of real 
operation processes of ship's engine rooms enables, already in 
its design stage, optimization of the time interval after elapsing 
of which one should carry out particular preventive operations. 
In this case one can consider two variants resulting from the 
necessity and possibility of considering two functions [4] : 


O function К, : T, > С, which describes the dependence be- 
tween the readiness factor of the investigated object (engi- 
ne room or any of its machinery or equipment) and time Т, 
elapsing from the moment of the start of the use of the ob- 
ject, being in condition of ability to the moment ofthe start 
of the preventive operation 


O function Ка: Tj > D, which describes the dependence be- 
tween the average income (or cost) per unit of time and 
time Тр. 


The first case results from the need for maximization of the 
readiness factor of the object (ship's engine room, its main en- 
gine or any other machinery or equipment), while the second – 
— from maximization of income and minimalization of cost. In 
these variants G and D mean in turn : readiness ofthe object to 
start carrying tasks and average income per unit of time. 

The correctness of hypothesis H1 can be confirmed in case 
of verification and acceptance of a more detailed hypothesis. 
Among such hypothesis is the hypothesis H2 explaining the 
fact observed in practice that wearing of the trybologic sliding 
systems of marine machinery (e.g. diesel engines) is weakly 
correlated with time [2,9,16]. This observation allowed to fo- 
recast the technical condition of machinery basing only on its 
current condition with omission of earlier conditions. Expla- 
nation of this fact can be presented in form of the following 
hypothesis (H2): the condition of any trybologic sliding sys- 
tem as well as the time of its duration considerably depend 
on the earlier condition and not on earlier conditions and 
time intervals of their duration because its load and both 
implied speed and increase of wear are processes having 
asymptotically independent values. 

The statement given in this hypothesis that load and both 
implied speed and increase of wear are processes having 
asymptotically independent values results from two obvious 
facts : 


€ there exists an exact dependence between load of the sli- 
ding trybologic systems and their wear [2,9] 


€ there is a lack of monotonous changes of the load of the 
sliding trybologic systems of machinery in longer periods 
of operation. Thus it can be assumed that the load of such 
systems is stationary [10,12,13]. 
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Verification of the introduced thesis H2 requires to deter- 
mine (foresee) the consequences, the occurrence of which can 
confirm the correctness ofthe formulated hypothesis. The con- 
sequences K2, which can be drawn from the hypothesis H2 are 
as follows [2,9,16] : 


О irregular wear of particular sliding trybologic systems 


О interlacing of realized wearing processes of the sliding try- 
bologic systems 

О the course of correlation function for a given sliding trybo- 
logic system such that with the increase of distance the func- 
tion first quickly decreases and then oscillates around zero 
with a relatively small but gradually smaller amplitude in 
function of distance 


О nearly normal distribution of increase of wear of the sli- 
ding trybologic systems in appropriately long time period 
of proper operation 


О linear dependence between the variance of the process of 
wear of sliding trybologic systems and the time of opera- 
tion. 


The mentioned consequences show the probabilistic law 
governing the wear of sliding trybologic systems. The verifi- 
cation of hypothesis H2 can be done in a similar way as in case 
of verification of hypothesis H1. However in this case one sho- 
uld accept the correctness of the following syntactic implica- 
tion : 


H2 => K2 (4) 


and then to apply noninductive (inductive) concluding 
in accordance with the following scheme : 


(K2, H2 > K2) H H2 (5) 


what is known as reduction conclusion. 
FINAL REMARKS AND CONCLUSIONS 


Carrying scientific research in order to prepare semi-Mar- 
kov models of engine rooms adequate for operation processes 
already in the design phase of these objects gives the possibili- 
ty to present to shipowners the concept of control of the above 
mentioned processes in accordance with optimization criterions, 
which are important for them. Among these criterions one can 
mention: maximum income, minimum operation investments, 
required (also maximum) technical readiness factor etc. To pre- 
pare this concept one can apply decision (controlled) semi- 
-Markov processes i.e. such semi-Markov processes in which 
realization depends on the decisions taken in particular mo- 
ments of changes of process conditions. The processes are sub- 
ject of considerations of semi-Markov decision (controlled) 
process theory. 


The control of the processes of engine room operation in 
the phase of operation enables maintaining proper durability, 
reliability and safety of functioning of the objects and of their 
ecologic features. 


Preparation of the above mentioned models of operation 
processes is equal to determination of a set of information, 
which should be at the disposal of designer, builder and opera- 
tor so that they could ensure rational control of real processes 
of such engineer rooms. Without doubt it will influence the 
progress in operation of engine rooms and thus the progress in 
shipbuilding technology. 
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A probabilistic model of environmental safety 
of ship power plant 


Roman Liberacki, M.Sc., Eng. 


ABSTRACT 


In the paper a probabilistic model of environmental risk of ship power plant is presented. 
This is a linear strategy model with an additional restraint. It can be used for assessing risk 
to marine environment, which results from ship power plant operation. The risk is a sum 
of component probabilities of exceedance of the limits assigned - by MARPOL convention - 
to particular kinds of pollution discharged from ships, multiplied by weight factors. The 
factors determine a degree of harmfulness of a discharge for the environment. The restraint 
concerns the case of exceedance of the limits assumed unacceptable by the convention. 


A risk value is be contained within (0,1) interval. Moreover a criterion for environmental safety of ship 
power plant was propsed. This is the criterion of the ALARP class, in which an intolerable risk level and 
acceptable one is distinguished. Suggestions concerning determination of the levels are submitted. The con- 
siderations are illustrated by results of example computer investigations of influence of reliability of techni- 
cal elements of the systems responsible for environmental safety of a hypothetical ship power plant on envi- 
ronmental risk value. Three qualification levels of ship 5 crew were accounted for : high, average and low. 


Keywords : environment protection, environment safety, risk, reliability, ship power plant 


INTRODUCTION 


Ships are technical objects which produce serious hazards 
to the seas and atmosphere. Importance of the problem 15 sup- 
ported by the fact of establishment of IMO MEPC Committee 
which deals with problems of pollution from ships to the seas. 
It has been estimated that maritime economy and shipping share 
in 20% in total pollution amount discharged to the marine envi- 
ronment [3]. 

When analyzing safety of a technical object one assumes 
that the object is a component ofthe man-technology-environ- 
ment system, and that losses which may be generated by the 
object within the system, are considered. The losses are un- 
avoidable, however one can influence a level of risk of their 
appearance. It should be reduced to a size which one would be 
willing to tolerate or even to accept [1]. Different safety crite- 
ria which define permissible risk levels for different categories 
of losses, are proposed. It also concerns the ships and risk in- 
volved to the marine environment. The risk mainly arises due 
to harmful substances released from ships to the seas and at- 
mosphere. It is obvious that environmental safety of a ship sho- 
uld be shaped when designing the ship. During its operation it 
can be only maintained at a certain level or corrected. 

International conventions adopted by IMO and ratified by 
member states, as well as other legal instruments of a lower 
rank impose definite limits upon harmful substances dischar- 
ged from ships. They impose necessity of installing special anti- 
pollution devices and complying with appropriate procedures 
of action to be undertaken by ship crews. The main legal in- 
strument in this respect is the 73/78 MARPOL Convention [11]. 

Norway is especially active in the area of marine environ- 
ment safety. This state elaborated a procedure called Environ- 
mental Indexing of Ships, submitted to IMO and implemented 
in the waters of its state responsibility [9]. The procedure ma- 
kes different kinds of dues paid by a ship, dependent on pollu- 
tion-preventing systems installed on it. 
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Both the requirements of the MARPOL Convention and 
the Norwegian procedure have one important drawback: they 
do not account for a random nature of risk created by ships to 
marine environment. The pollution-preventing systems are 
anthropotechnical, and their technical and human components 
have definite reliability. The fact that a ship 1s equipped with 
devices complying with the standards imposed by the Conven- 
tion does not show itself that the standards are also fulfilled 
during ship's operation because the devices can be unservica- 
ble, wrongly used or simply switched off. Therefore a propo- 
sal is to supplement the deterministic requirements with pro- 
babilistic methods for assessing reliability and environmental 
risk of ships. The probabilistic model of environmental safety 
of ship power plant, elaborated by this author and presented in 
the paper, may serve as such a tool. 


PROBABILISTIC MODEL 
OF ENVIRONMENTAL SAFETY 
OF SHIP POWER PLANT 


In the environmental safety model of ship power plant the 
kinds of pollutions which are dealt with by the MARPOL Con- 
vention standards, were taken into account. These are : oil pol- 
lutants, sewage, noxious components of exhaust gas emitted 
by piston diesel engines and waste incinerators, as well as co- 
oling media ( freons). 

The environmental safety of ship power plant 1s expressed 
by risk level it creates for the environment. The lower risk the 
greater safety. The risk measure associates probabilities of cau- 
sing hazards to environment and their consequences into a sin- 
gle numerical value. A caused hazard is the emission of such 
pollution amount which exceeds its limit value imposed by the 
MARPOL Convention. 

Noxious substances released to the environment during 
ship's power plant operation are of different forms of harmful 
influence on the environment. Due to complexity of the pro- 
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blem, to use weight factors for determining the consequences 
of released pollutions is most suitable. To each kind of pollu- 
tion a weight factor reflecting seriousness of environmental 
losses due to its discharge, was assigned. 

The probabilities and weight factors were joined together 
into a single value by means ofthe linear strategy model [5], to 
which a restraint was added. The strategy is of a compensating 
character. A high value of one component of the model is ca- 
pable to compensate a low value of the other. Whereas the re- 
straint has to select the ship power plants not complying with 
the Convention's requirements as those hazardous. An addi- 
tional advantage of the assumed linear strategy is the simplici- 
ty resulting from its additive character. The environmental risk 
measure of ship power plant can be expressed as follows : 


m m 
1=1 i=l (1) 


where : 


Xj — amount of discharge of i-th kind of pollu- 
tion (random variable) 
Xil — limit amount of discharge of i-th kind of 
pollution, imposed by the Convention 
m — number of kinds of pollution, taken into 
account 
Wi — weight factor for i-th kind of pollution 
Pi(Xi > xj) - probability that the discharge amount of 
i-th kind of pollution will exceed the year- 
ly limit amount of discharge permitted by 
the Convention. 
It is worth observing that : 
3i Р(Х, >x,)=1>R=1 
The first relationship corresponds with the situation when 
a ship does not comply with any of the Convention's require- 
ments, i.e. it is not fitted with appropriate devices and systems. 
The maximum risk value is then equal to R = 1. 
The second relationship corresponds with the situation when 
a ship and its crew are beyond reproach regarding environ- 
mental safety. The situation is obviously not real; then R = 0. 
Therefore the risk measure R is contained within the interval 
between 0 and 1. 
The probabilities РХ; > xij), 1 = 1,2,...,m, which appear in 
(1) can be determined on the basis of singular models of envi- 
ronmental safety of ship power plant. 12 such models were 
elaborated : 


for operation of fuel oil delivery to ship 

for operation of fuel oil pumping between tanks onboard 
the ship 

for operation of lubricating oil delivery to ship 

for operation of lubricating oil pumping between tanks on- 
board the ship 

for operation of transferring oily substances out of the ship 
for packing (sealing) of stern tube,controllable pitch pro- 
peller boss, and controllable pitch propeller boss of thru- 
ster 

for bilge water deoiling system 

for ship sewage system 

for ship diesel engines regarding their NO, emission 

for ship diesel engines regarding their SO, emission 

for ship incinerator 

for ship refrigerating system. 


+ + жж жж 


жж 


For every kind of pollution a review of design, features, 
structure and modes of operation of a relevant technical sys- 
tem, was performed with taking into account crew's actions. 
Next, appropriate unserviceability and event trees were elabo- 
rated. This way cause - effect models were obtained. Each of 
the models contains — in its structure comprising the above 
mentioned trees - description of a scenario of events leading to 
triggering a hazard. A triggered hazard is the emission of such 
pollution amount which exceeds its limiting value imposed by 
the MARPOL Convention. The specific models make it possi- 
ble to determine probabilities of violating the Convention con- 
ditions, i.e. oftriggering the considered hazards. In the models 
any conscious, deliberate, unlegal pollution discharges by ship 
crews were not taken into account. Input values to the singular 
models are the unserviceability events of technical elements 
and human errors. The model tree structures are compliant with 
recommendations of the FSA method being developed within 
IMO activity [2, 12, 13]. 

To determine numerical values of the weight factors the 
Norwegian report on the Environmental Indexing of Ships was 
used [9]. However the there given values of weight factors were 
not applied in this work, because such pollutions as bilge wa- 
ter, lubricating oils, and oil residues were not accounted for in 
it. The report's authors asssumed that 1f the discharging to the 
sea a given kind of pollution is prohibited then it is not dischar- 
ged in fact. They also assumed that the pollutions subjected to 
treatment (bilge water) are always cleaned up to a sufficient 
degree. They do not take also into account oil leakage from 
stern tube as it should be prevented by packing. This approach 
results from a deterministic viewpoint on the environment pro- 
tection systems, without accounting for reliability of technical 
elements and human relability. 

However in the Norwegian report in question, a very detail 
analysis of influence of different pollutants on the environment, 
was performed. Their impact to the environment, which re- 
sults from their destructive mechanisms to the environment, 
range of influence and potential amount of discharge, was the- 
re linguistically described. 

Values ofthe weight factors of detrimental influence of pol- 
lutants on the marine environment are given in Tab.1. The fac- 
tors were determined by changing the linguistic variables cha- 
racterizing the environmental impact into appropriate numeri- 
cal values (col. 4 and 5). Next, for each kind of pollution the 
assigned numerical values (col.5) were summed up for parti- 
cular destructive impacts as shown in col. 6. Finally, the values 
of col. 6 for particular kinds of pollution were divided by the 
total sum of the numerical values determined for all kinds of 
pollution (2: =72). This way the searched values of the weight 
factors, w; , (col. 7) were found. The sum of all weight factors 
equals one in accordance with the formula (1). 

Incinerators are only used at high sea, far from land. Hence 
their influence is very limited. Amount of their exhaust gases 
18 small relative to that emitted from the combustion engines. 


INVESTIGATION OF THE MODEL 


The investigating subject was a hypothetical ship power 
plant of a typical technical structure complying with contem- 
porary solutions. Its operational process also complied with 
that presently applied. To perform calculations the following 
assumptions were made : 


2 Regarding risk of NO, excessive emission : 
in the power plant 4 combustion piston engines are instal- 
led : 1 main engine, and 3 auxiliary engines driving inde- 
pendent electric generating sets; no systems for purifying 
exhaust gas from NO, are applied. 
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Tab.1. Values of weight factors of harmful influence of polluting agents on marine environment У = 72 Х-1 


Climatic changes 
Acid rains 
Production of ozon 


Destruction of ozone layer 
Climatic changes 
Exhaust 
gases from V 
incinerators 
V 


Numerical 
value of 
environmental 


Sum of Weight 
numerical factor 
values 


Heavy fuel Pollution of beaches High to very high | 4 to 5 (4.5) 
oil Pollution of sea High to very high | 4to 5 (4.5) 


Lubricating Pollution of beaches Moderate 2 
; VII j 
oil Influence on sea organisms Moderate 2 


Pollution of beaches 


Fuel oil VI Pollution of beaches Mean to high 3to4 (3.5) 
(diesel) Pollution of sea Mean to high 3to4 (3.5) 


Influence on sea organisms 
Pollution of beaches 
Influence on sea organisms 


Pollution of beaches 
Sewage X ; 
Influence on sea organisms 


Low to high 1104 (2.5) 
Low to high 1104 (2.5) 
Very high 5 


* Harmfulness of exhaust gases from incinerators was not determined by the Norwegian authors. 


Regarding risk of SO, excessive emission : 

the ship in question 24 times a year enters a protected re- 
gion in which using a fuel oil having sulphur mass content 
lower than 1,596 15 required; beyond the region a fuel oil of 
sulphur content exceeding 1,5% is used on the ship; no sys- 
tems for purifying exhaust gas from SO, are applied; ship's 
fuel system is common for all the engines (a change of fuel 
oil kind for one engine simultaneously applies to the remai- 
ning engines); hence the set of engines ofthe power plant in 
question should be considered as a single engine. 


Regarding risk of freon release to the atmosphere : 

in the ship's refrigerating system a freon of HCFC group is 
used; repair and maintenance operations are carried out once 
a year, during which a release of freon may happen. 


Regarding risk of of excessive emission of noxious sub- 
stances resulting from the use of a waste incinerator : 
waste combusting operations by means of the incinerator 
are performed 15 times a year. 


Regarding risk of heavy fuel oil release to the sea : 

the heavy fuel oil bunkering operations are carried out 12 
times a year on the ship; the heavy fuel oil pumping opera- 
tions between the ship's tanks are performed 300 times 
a year. 


Regarding risk of (diesel) fuel oil release to the sea : 
the fuel oil bunkering operations are carried out 2 times 
a year on the ship; the fuel oil pumping operations between 
the ship's tanks are performed 24 times a year. 


Regarding risk of lubricating oil leakage to the sea : 
the lubricating oil bunkering operations are carried out 4 
times a year by means of oil supply piping; the filled tank 
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18 fitted with an air-escape pipe reaching over the deck; the 
oil pumping operations between the ship's tanks do not in- 
volve any risk of oil leakage to the sea; the ship 15 equipped 
neither with a controllable pitch propeller nor thrusters; oil 
leakage may only occur through stern tube packing. 


2 Regarding risk of leakage of oil residues to the sea : 
the operations of transferring oil derivatives out ofthe ship 
are carried out 4 times a year; the ship is fitted with the oil/ 
bilge water separator NEPTUN made by Warma, Grudziadz 
(Poland); the separated oil collected in the separator may 
be released to the sea due to a failure of oil discharge sys- 
tem of the 274 chamber of the separator. 


> Regarding risk of release of bilge water containing more 
than 15 ppm of oil : 
the ship is fitted with the bilge water/oil separator NEP- 
TUN; in the engine room cleaning operations with the use 
of detergents are performed 6 times a year. 


2 Regarding risk of sewage pollution of coastal waters : 
the ship is fitted with LK biological-chemical sewage tre- 
atment plant, Warma, Grudziadz; the probability of ship's 
sailing over sea waters up to 12 NM from the nearest land 
amounts to 30 % ; the seawage treatment plant is operated 
in accordance with DTR recommendations. 


Way of performing the model investigation 


Three levels of crew qualification were assumed: high, ave- 
rage and low. Probabilities of making human errors were de- 
termined by using tables elaborated on the basis ofthe THERP 
method [8] as well as calculations performed by using ASEP – 
— HRAP method [4] and TESEO method [6], [14]. 
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The same reliability level Rye was assumed for all technical 
elements of the investigated systems. For elements of the sin- 
gular operation models of : bunkering, pumping, transferring 
out fuel oils, lubricating oils and oil derivatives, of combusting 
wastes as well as of entering a region protected against SOx 
emission, the reliability level concerns a single operation only. 
For technical elements of singular models of : packing, bilge 
water treatment system, sewage system, combustion engines 
regarding NOx emission, as well as refrigerating system, the 
yearly reliability level is assumed. For elements of automation 
and signalling systems the so called reliability on demand is 
assumed. 

The investigation consisted in simultaneous changing re- 
liability levels of all technical elements within the range from 
0,9 to 1. On the basis of relevant singular models as well as 
values of the weight factors of harmfulness of pollution impact 
onto the environment, values of yearly environmental risk were 
calculated in function of the reliabilities of technical elements 
and crew qualification levels. Results of the calculations are 
shown in Fig.1. 

0,45 - 


0.4 


R(per year) 


0 T T 

0.9 0.91 0.92 0.93 0.94 0.95 0.96 0.97 0.98 0.99 1 
Rie 

Fig. I. Influence of the reliability level of technical elements, Rie, and crew 

qualification levels on the yearly environmental risk of ship power plant, R. 


Conclusion derived 
from the example model investigation 


* Theenvironmental risk of ship power plant, defined by the 
expression (1), obviously depends on the reliability of rele- 
vant technical elements as well as on the crew qualification 
level - the reliability of human operator. 


* The greater reliability of technical elements and human re- 
liability the lower environmental risk from the side of ship 
power plant, i.e. the higher level of environmental safety. 


A PROPOSED CRITERION 
FOR ENVIRONMENTAL RISK 
OF SHIP POWER PLANT 


The proposed risk criterion follows the proposal elabora- 
ted in the frame of FSA method used by IMO. The entire risk 
range is divided into three categories : 


the intolerated risk 

negligible risk 

the so called ALARP (As Low As Reasonably Practicable) 
range. 


The intolerated risk cannot be accepted in any case. The 
negligible risk is acceptable and no way of its mitigation is 
searched for. The risk within the ALARP range is subjected to 
an analysis. Possible lowering the risk level and costs associa- 


ted with the possibility are considered. If lowering the risk le- 
vel is practically possible and its cost is reasonable then ac- 
tions aimed at its reduction are undertaken [7,13]. 

The division of the risk range into three categories requires 
two limit risk values : upper Ry and lower Еј , to be establi- 
shed. IMO carries out work on establishing such limits in the 
case of human life. They reflects risk produced by ships to 
human life. As far as the risk to the environment due to ship 
power plant is concerned such limits still have to be determi- 
ned. A graphical interpretation ofthe risk criterion 15 shown in 
Fig.2. 


Intolerated risk range 


R, - upper risk limit 
ALARP range 


R, - lower risk limit 
Negligible 
risk range 


Fig. 2. A graphical interpretation of the risk criterion 


Here it is proposed the upper risk limit Ry to be a value 
somewhat higher than the average risk according to the formu- 
la (1) for power plants of ships in service, under assumption of 
the average qualification level of ship crews. The proposal is 
based on the author's conviction that one should strive for im- 
proving present state of environmental safety. On the other hand 
it cannot be said that the risk level due to existing ship power 
plants complying with the international requirements, is not 
admissible. 

A value of the lower risk limit Rj can be also determined by 
using the formula (1), in this case for a ship power plant desi- 
gned with taking into account possibly proecological solutions 
and ways of its operation, following the Norwegian ship's class 
notation Clean Design [10]. The notation does not take into 
account reliability of technical elements which should be on 
the level as high as possible. The human reliability would be 
assumed for the high qualification level. 

In this case a proecological ship represents the ship the stern 
tube of which is filled with water instead of lubricating oil, on 
which oily bilge water and sewage are collected in tanks and 
transferred out to land — based stations, only diesel fuel oil 
containing less than 1.5% of sulphur is used, and freons are not 
applied to its refrigerating system. 


FINAL REMARKS AND CONCLUSIONS 


* Randomness of the event of inadmissible pollution release 
from ship power plant is obvious. It implicates applying 
probabilistic approach to its environmental safety. The sa- 
fety is a feature of the power plant which has to be assessed 
by using the notion of risk and safety criterion defined in 
safety engineering. 


** Application of the elaborated models is closely associated 
with data on reliability of elements of environment protec- 
tion systems; however they are unavailable in Poland. Per- 
haps they are saved in databases of Det Norske Veritas or 
Lloyds Register. Also experts — experienced ship engineers 
— have not sufficient experience to base on their opinion in 
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formulating credible assessments. Ship environmental pro- 
tection systems are considered by ship engineers as auxilia- 
ry ones only, moreover they subject to failures rather rarely. 


* The elaborated models make it possible to formulate relia- 
bility requirements for elements of environmental protec- 
tion systems. The requirements can be put forward to pro- 
ducers ofthe elements, which are able to appropriately shape 
reliability of the elements as well as to verify it by means of 
laboratory tests or investigations in service. The problem is 
covered by PN-IEC 60300-3-4 standard. It consists only in 
formulating appropriate requirements and, certainly, in as- 
sociated costs of their implementation. 


** The models make it also possible to design reliability of 
the power plant systems responsible for releasing noxious 
substances – to a given amount only, as well as the entire 
power plant — to a given level of environmental risk. 


% The environmental safety of ship power plant is a partial 
problem. The environmental safety ofthe entire ship is that 
mainly interesting. It seems that just the so defined models 
have a chance to be applied in the worldwide shipping. 
Hence the presented models exemplify only a methodical 
approach to solving the problem of environmental safety 
assessment. 
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A concept of the predicting 
of technical state of devices 
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ABSTRACT 


The paper presents possibilities and limitations associated with technical state predicting 
procedure of devices, as well as it indicates the necessity of improving the present schedu- 
ling methods of preventive maintenance operations of ship power plant devices. It has 
been also showed that by making transformation of the wear rate distributions of device 
elements into the distributions of time of their correct- to-limit-state operation it is possi- 
ble to predict their technical state and to suitably schedule their maintenance operations. 


Keywords : technical state, prediction, maintenance, mathematical transformations, probability distributions 


INTRODUCTION 


In the process of predicting technical state of a device, in 
a result of which an optimum prediction containing a date and 
range of the next testing of the state 15 formulated, the device's 
user generates information which makes it possible to control 
keeping the device in the state of serviceability. An essential 
element of such control is to know prediction methods of devi- 
ce's state, containing procedures for scheduling a date and ran- 
ge ofthe next testing and the possible scope of overhaul opera- 
tions resulting from it. 

Due to a complexity of ship devices, their profitability and 
operational safety, the preventive maintenance strategy realized 
on the basis of their wear value is the most important. In order 
to elaborate such strategy it is necessary to choose an appro- 
priate strategy for predicting working time resources of a given 
object with taking into account its service conditions [9,11]. 

When analyzing prediction methods it is not possible to 
unambiguously show a greater suitability of any of them aga- 
inst other ones because it depends on which kind of objects 15 
subjected to predicting investigations [1,6]. 

In the present phase of knowledge and practical experience 
dealing with predicting problems in diagnostics, it is today hard 
to select an algorithmic procedure which would contain clear- 
-cut recommendations and would exactly precise which me- 
thods have to be used for elaborating an optimum prediction in 
given conditions. Experience shows that in very many situ- 
ations there is no general method which could be applied as 
a rule in realizing the process of choosing right diagnostic pre- 
dictions. 

Ifthe ageing ofthe object is considered as a set of results of 
actions of processes limiting its service effectiveness, a question 
arises if it is possible to predict amount of wear of a repairable 
object and to assess cost-effectiveness of its further using with- 
out its replacement or carrying-out its major repair. 

For this reason the attempt to present possibilities and limi- 
tations which concern the predicting of technical state of devi- 
ces, as well as to indicate a necessity of improving the current 
methods of planning preventive maintenance operations of ship 
power plant devices, has been undertaken. 


One of the possibilities may be to elaborate models of fai- 
lure predicting and scheduling maintenance operations, in which 
wear phenomena of device elements would be accounted for. 
A way of transforming the wear rate distributions of device 
elements into the distributions of time of their correct-to-limit- 
-state operation, was demonstrated. It makes it possible to esti- 
mate service life indices for devices, e.g. their working time 
resources, and to schedule their maintenance operations ap- 
propriately. 


POSSIBILITES AND LIMITATIONS 
OF THE PREDICTING OF STATE 
OF DEVICES 


Physical ageing of devices usually occurs due to destructi- 
ve processes taking place in them, associated with friction, 
corrosion, errosion, cavitation and fatigue of structural mate- 
rials of their elements. Along with time of their operation it 
leads to often and often occurring failures and in consequence 
to lowering their reliability. 

Service reliability investigations are aimed at determina- 
tion of reliability function of an object, changeable with time 
of its operation, on the basis of data obtained in service. 

Any decrease of reliability level cannot exceed a permissi- 
ble value, hence it should be assumed that 1f such event hap- 
pens the object in question will be taken out of service. 

Below has been undertaken an attempt to present some es- 
sential possibilities and limitations associated with the predic- 
ting of reliability of devices , in order to show that searching 
for new predicting methods is necessary. 

The irreversible wear processes occurring in devices in ser- 
vice cause a monotonically changeable trend of values of diag- 
nostic control parameters to occur. 

A basic condition to succeed in predicting the state of devi- 
ces is to assume a uniform wear rate of their elements, i.e. that 
a trend of measured quantities is known (failure rate - in case 
of reliability considerations). 

Changes of values of diagnostic parameters obtained from 
particular state investigations may differ seriously, and their 
probability density functions are usually not known. 
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Therfore it can be assumed that : 


9 the scheduling ofthe successive diagnostic test of a device 
18 possible as а result ofthe technical state prediction which 
consists in determination of future changes of diagnostic 
parameters and in comparison their instantaneous values 
with those limiting [12,13,15] 


€ assessment of a trend depends on all observations of a va- 
riable, but a weight attributed to relevant realizations de- 
creases in function of the time passing from the instant of 
Observation up to now. It means that the longer time ago an 
event has occurred the smaller the weight, and that its im- 
portance decreases with the time passing from the instant 
of occurence of the event [5]. 


The main tasks for realization of the so described 
prediction process are the following : 


X selection of optimum diagnostic parameters describing 
a current state and its change during operation of a device 


х determination of a predicted value of a diagnostic parame- 
ter behind the time horizon of the diagnosis, by using an 
optimum diagnostic method 


X scheduling the time of the sucessive diagnostic test. 


When considering the real operation processes one has to 
do with both deterministic and random processes. The random 
processes cause that [8] : 


** to formally describe phenomena is not possible or very com- 
plicated 

** theconsidered phenomena are non-measurable, by conven- 
tion 

ж the considered objects and phenomena are subject to conti- 
nuous changes. 


For this reason the following components appear in descrip- 
tion of the real processes which occur during operation of 
a device : 


+ deterministic ones — possible to be calculated by using in- 
ductive methods 

+ stochastic ones — possible to be determined by using de- 
ductive methods on the basis of statistical data, 

5> or not being a basis for any predictions. 


The reliability theory contains only a few theoretical mo- 
dels whose application makes it possible to obtain an exact 
description of influence of different factors on reliability cha- 
racteristics of technical objects; therefore further elaborating 
and testing the mathematical models for reliability investiga- 
tions of such objects, is justified. 

However in any case an empirical mathematical model is 
simplified and it belongs to an a priori determined class. Usa- 
bility and properties of some mathematical models important 
for reliability investigations of ship equipment are presented in 
[6]. In the case of ship equipment traditional methods for veri- 
fication of statistical hypotheses by means of tests may be not 
sufficient to select the most suitable model. Knowledge of 
a kind of device failure (e.g. sudden or progressive) can facili- 
tate the selecting of an appropriate mathematical model [6,17]. 

The systems for reliability prediction of objects, described 
in the subject-matter literature and used in practice, are based 
on investigations of influence of many factors which determi- 
ne failure rate analytically described in an empirical model. 
Accuracy of such models and credibility of reliability predic- 
tions resulting from them, depends mainly on accuracy of em- 
pirical data obtained from laboratory tests or service investiga- 
tions. 
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Wear phenomena in object elements are very complex and 
many factors influence their course therefore to take all of them 
into account is impossible; hence to make choice of an opti- 
mum set of diagnostic parameters is difficult. 

Predicting the technical state of machines on the basis of 
values of diagnostic parameters changeable in service and re- 
lated to a longer time period, is associated with a risk that 
the predicting diagnosis would be based on an out-of-date 
model, i.e. that whose elements (an optimum set of diagnostic 
parameters and an optimum prediction method) do not reflect 
real relationships between the technical state ofthe system and 
the predicting diagnosis, any longer . Hence the optimum pre- 
dicting diagnosis should be stable within the entire period of 
state predicting. If it appears unstable then it may be applied to 
a dynamical operation system of machines. Otherwise should 
be made a decision to modify assumptions and limitations, e.g. 
by deliberate not-accounting for factors causing unstability and 
decreasing versatility of the obtained solution. 

If in a performed predicting action the statistical informa- 
tion required for the used model is unavailable then it becomes 
necessary to apply a testing procedure which makes it possible 
to test in conditions of incomplete data. It is usually associated 
with the necessity to limit the range of predicting analysis to 
a short-range horizon. 

Diagnostic prediction not always provides results expected 
from the point of view of its rightness and accuracy. The pre- 
diction rightness can be determined by probability of fulfil- 
ment of the prediction. The more fully and objectively are pa- 
rametrized realization conditions of a process the greater the 
prediction rightness. 

In the context of the presented possibilities and limitations 
It is necessary to search for ways to improve current methods 
of planning. One of them can be to elaborate such models in 
which wear phenomena of the object's elements would be ta- 
ken into account. To this end may serve the transformation of 
the wear rate distributions of device elements into the distribu- 
tions of their correct-to-limit-state operation, which would 
make predicting their technical state and appropriate schedu- 
ling their maintenance operations, possible. 


TRANSFORMATION OF THE WEAR RATE 
DISTRIBUTION MODELS INTO THE 
MODELS OF TIME OF CORRECT- 
-TO-LIMIT-STATE OPERATION 
OF DEVICE ELEMENTS 


To assess service life of devices at the initial stage of their 
operation is only possible by analyzing their wear processes, 
to this end failure predicting models accounting for wear phe- 
nomena of device elements, are the most suitable [14]. 


The reliability model of element’s wear can be presented 
in the following form [14] : 


t со 
[аса = [64 — 120 @) 
0 Z 


g 


where : 
fi(t;)— density function of the serviceability time T, ог: 
0<6<1 
f,(z) – density function of the wear z for the fixed va- 
luetz 0 
t, = current duration time ofthe process of the wear 2. 


The formula (1) expresses equality of the probability ofthe 
event consisting in that in the instant ¢ or earlier a failure of the 
object occurred, and the probability that in the instant ¢ the 
object attained the wear value z, or greater. The relationship 
(1) is illustrated in Fig.1 [14]. 
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Fig.1. Schematic diagram of service life assessment on the basis of wear value taken as a random variable 


Making use of the equality (1) one can perform service life 
assessment by determining the density function fj(t) on the 
basis of distribution characteristics of the wear described by 
the function f(z). 

If the function f>(z) for the instant t is determined on the 
basis of observations of changes of wear distribution para- 
meters, performed in the instants earlier than ¢, then by ap- 
plying the equality (1) one can predict object’s service life. 
To perform such prediction it is necessary to know the allo- 
wable wear value Zg and stochastic model of wear process 
(see [14]). 

A wear rate analysis of elements of ship devices shows that 
their distributions can be approximated by some theoretical 
distributions : normal, Rayleigh's and Weibull's one. The two- 
-parameter Weibull distribution is the most versatile. In some 
special cases the exponential distribution can be also applied 
[2.4,6,7,15]. 

The wear rate distribution empirically or theoretically de- 
termined can be transformed into the distribution of time of 
correct-to-limit-state operation. By assuming that the wear rate 
distribution g(u) was obtained for the sufficiently large set of 
по monomial elements, the number of elements n(T) which 


reached the limit state within the time t = T can be determined 
by means of the following equation : 


n(T) 2 n,-F(T) Q) 
If G(u) = РОЈ; < u) is the cumulative distribution function 


of wear rate, and g(u) its density, then the unreliability func- 
tion F(t) can be expressed as follows [2,3] : 


F(t) 2 P(T«t) =Р "E 24 |= 
zm (3) 
7, 
=P es = [2 du 
z, [t 
Hence : 
z, /t 
R(t) =1-F() = fgu)du (4) 


Tab.1 presents the solutions of the equation (3) for the den- 
sity functions of wear rate g(u) according to the normal, Ray- 
leigh, Weibull and exponential distributions. 


Tab.1. Transformation of wear rate distributions into the distributions of time to correct- to-limit-state operation 


Wear rate distributions of elements 


Distributions of time to correct- to-limit-state operation 


Distribution density 
function 


g(u) 


Probability of 
reaching the limit 
state, F(t) 


Distribution | Parameters 


Normal 


Failure rate 
f(t) 
1- F(t) 


Probability density function 
dF(t) 


Mt) = беја ЕЕ 


Rayleigh's 


Weibull's 


Exponential 
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For instance the way of transforming the density function 
of wear rate g(u) according to Weibull distribution is the follo- 
wing [2,3] : 


© ЕГІ k-1 ü k 
F(t) = feudu = ЕР ехр 4% ди = 


а 
2, [t 2/ 


п (5) 
= –ех EN > =eXp| – = 
p а z, ft p at 
Making use of the relationships : 
dF(t) 
f(t) = —— 6 
(t) E (6) 
and : 
f(t 
A(t) = кя (7) 
1- ЕЦ) 


one obtained the expressions for the distribution density func- 
tion of time to correct-to-limit-state operation, f(t), and the 
failure rate function A(t), presented in Tab.1. This way the wear 
rate distributions were transformed into the distributions of 
time to correct-to-limit-state operation of elements. The para- 
meters of the wear rate distributions and the allowable limit 
wear value of element, z; , become the parameters of the new 
distributions. 


If the parameter W : 
22 
Ұ = = (8) 
u 


is introduced and the expressions for the expected wear rate 
values u for particular distributions are assumed, the following 
forms of the function f(t) are obtained : 


for normal distribution of wear rate, 
where: z, = W.u (9) 


б СЕ. ехр Eu At | 
J/2x ot? 2! ot 


(10) 
for exponential distribution, 
where: z, = W-u п =! (11) 
rose- | (12) 
t t 
for Rayleigh's distribution, 
where: U = mo (13) 
ые 
л ХУ? n(WY 
f(t) = ——exp|-—| — 14 
(t) 21 | А | (14) 
for Weibull’s distribution, 
where : 0 (15) 
крк 1 
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Because the Weibull’s distribution is the most versatile, the 
equation (16) may be taken as the general form of f(t) which — 
for selected values of the parameter k — yields particular solu- 
tions. 

Hence : 


> for the parameter value k = 1 the equation (16) is transfor- 
med into the equation (12) 
> for k=2 the equation (14) is obtained. 


The equation (16) is the probability density function f(t) 
for the two-parameter (W, k) distribution of time of correct-to- 
-limit-state operation , shortly called the W distribution [2]. 

In this case the reliability function R(t) which expresses the 
probability that a device will not reach its limit state, can be 


represented as follows : " 
] 
WI| —+1 
k 
t 


The function R(t) may be used to determine, a.o., y-percent 
resources of working time of devices, hence for scheduling 
their maintenance operations. A good conformity of the course 
of the function R(t) and results of statistical investigations on 
diesel engine cylinder liners is shown in [3]. Therefore, depen- 
ding on information at one’s disposal, service data or theoreti- 
cal relationships can be used for assessing the process of run- 
ning down the working time resources of device elements. 


ASSESSING THE PROCESS OF RUNNING 
DOWN THE WORKING TIME RESOURCES 
OF DEVICE ELEMENTS 


t 
R(t) -1- ff (dt =1-exp4- (17) 
0 


Working time of a device, after which its working time reso- 
urce is exhausted, can be determined on the basis of the working 
time resource specified by the device’s producer or that deter- 
mined with the use of statistical data obtained from service. 

A course of running down the working time resource of 
any device element is exemplified in Fig.2, at accounting for 
the following assumptions [11] : 


> the worn-out element (that of exhausted time resource) is 
replaced with a new one 

> at the beginning of the service process or just after replace- 
ment of an element the ratio of the working time resource f; 
and the time resource f, (specified by the device’s produ- 
cer) equals 1, and it equals 0 when the working time reso- 
urce is run out, i.e the limit state is attained 

2 the running-down process of working time resource of an 
element in service, is linear. 


The diagram shown in Fig.2 presents the way of determi- 
ning the time after which an element is to be replaced with 
a new one or subjected to major repair. 

In the course of operation one does not usually let the limit 
state of device elements to occur, therefore their replacement 
is made in advance. In Fig.2 such procedure is represented by 
dotted lines. Technical state of an entirely used up (worn out) 
element of the device whose values of control parameters are 
contained within the interval between allowable and limit ones 
[2], is so low that further use of the device in which such ele- 
ment is installed is associated with greater and greater expen- 
ditures. Hence the in-advance replacement of an element with 
a new one as a matter of fact shortens its working time, how- 
ever it also makes it possible to restore the demanded technical 
state of the device and to reduce its operational cost. 
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Fig.2. Running-down process of working time resource of an object 


In the case of replacement of the entire device with a new 
one the system (e.g. power plant) in which it has been instal- 
led, becomes restored , however if only one of its elements is 
replaced then the device (or power plant) becomes restored 
partly. Therefore it can be assumed that in the first case the 
device's technical state and working time resource become re- 
stored to their maximum values, whereas 1f the object is sub- 
jected to overhaul its technical state and working time resour- 
ce can be maintained at a given level or restored , however 
they will no longer reach their initial values (those at the be- 
ginning of the device's use) and the device in question will be 
no longer considered as „а new one". It means that in the case 
of repaired objects the line of the maximum relative working 
time resource shown in Fig. 2 will monotonically descend. The 
entire problem in question is presented a.o. in [12]. 

By applying the way of the determining of working time 
resource, presented in Fig. 2, great advantages are obtained by 
means of analysis of running-down process of working time 

Я 


Injector 


Injection 
pump 


resources of elements of complex technical objects, e.g. diesel 
engine, in order to answer the question : which the engine's 
elements require to be serviced, and after which time of 
operation (see Fig.3). Similar analysis can be performed for 
e.g. ship power plant. As a result of it specification of the ob- 
jects which require maintenance operations in particular years 
of ship service would be obtained. 

Various possible schedules of servicing can be assumed, 
e.g. after every voyage, each year etc; in each case the scope, 
time and cost of servicing will be different. It will also depend 
on service conditions of particular objects, which greatly in- 
fluence working time resources of the objects. 

For instance, Fig.4 presents the way of determining the re- 
sources of working time to preventive maintenance, routine 
repair and major repair, counted from a given service year £4 , 
at different service conditions. 

In the example shown in Fig.4 the object's serviceability 
time was divided into equal intervals in such a way that the 
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Fig.3. Running-down process of working time resources of combustion engine elements (an example) 
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Fig.4. Influence of service conditions on working time resource of an object and on scheduling its servicing operations 
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time intervals to preventive maintenance, routine repair and 
major repair were equal (usually given by device producers or 
shipowner). Any change of service conditions would result in 
changing the working time resource of the device in question , 
manifested by shifting-up the line of running-down the wor- 
king time resource of the device, shown in Fig.4. 

For instance, an improvment of service conditions stands 
for the passing from the point A to the point B and for an ap- 
propriate increase of the distances x, y and z determined relati- 
ve to a new line of running-down the working time resource. 
As aresult it would be necessary to modify the initial schedule 
of servicing the device. 


FINAL REMARKS 


% Inthe subject-matter literature there are no established pro- 
cedures for the determining of working time resources of 
technical objects in changeable service conditions. 


** The presented concept of the predicting of working time 
resources may facilitate the scheduling of preventive ma- 
intenance of objects and the choosing of optimum servi- 
cing plan. 


** Usefulness ofa given prediction method increases if a com- 
prehensive physical analysis of development of wear pro- 
cesses and their features are performed in advance. 


For all prediction realizations it is important to know history of 
development of the to-be-predicted process. Moreover, an ana- 
lysis of the conditions which determined past course of the 
process, assessment of significance of structural components 
influencing its course, as well as estimation of possible influ- 
ences in the future, is also important. 


** Wear rate distributions of device elements can be used for 
determining distributions of their working time to limit sta- 
te and in consequence to determine service life indices cha- 
racterizing the devices e.g. their working time resources. 


** From theoretical and practical point of view the presented 
distribution W can be useful in describing progressive fai- 
lures. 


NOMENCLATURE 

а - scale parameter of Weibull distribution 

f(t) — probability density function of time to limit state (time to 
failure) 


F(t) — failure probability function (unreliability function), 
cummulative distribution function of serviceability time 

g(u) — wear rate probability density function 

h(z) — wear-out probability density function 

k shape parameter of Weibull distribution 

ny, -— number of objects in question 

n(t) — number of objects failed to the instant 7 

R(t) – reliability function 

t = time 

t, -— time of correct operation of object at y-percent probability 
(Y-percent time resource) 


Т - time to limit state (random variable) 

u — wear rate (realization); u = z/t 

Up = limit wear rate; = zy 

U  — expected wear rate value 

U - wear rate (random variable); U = Z/t 

U, = limit wear rate (random variable); U, = z,/T 

W -— expected service life value; W = z,/u 

z  — wear (realization); increment of control parameter value 
7а - initial (lower) wear value 

Zg - limit (upper) wear value 
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2 - wear (random value) 

à  — failure rate, exponential distribution parameter 
A(t) — failure rate function 

с - normal distribution parameter; standard deviation. 
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ABSTRACT 


The paper concerns a propulsion system of merchant ships intended for sailing in the Baltic Sea zone. Such 
system is to satisfy the ecological requirements determined by relevant international conventions for special 
zones to which the Baltic Sea also belongs. The paper draws attention to gas turbine used as a prime mover 
for such ships, because it satisfies the ecological requirements and has also other advantages. Application 
of gas turbine for ship powering does not require exhaust gas to be purified, however it requires fuel oils 
ofa low sulphur content to be used. If the ecological rules impose the using of the fuel oils of similar quality 
for diesel engines then gas turbine propulsion system will be comparable — also economically (regarding 
specific fuel oil consumption cost) — with that of diesel engine. It would be even more favourable in 
a combine gas turbine /steam turbine system, especially at compound production of electric and heat energy 
(i.e. COGES systems). In the Baltic Sea zone gas turbines will find application to powering a. o. such ships 
as : fast car-passenger ferries, fast cargo ships, special vehicles (hydrofoils, hovercraft, motor yachts). 


Keywords : ship propulsion system, gas turbines, environment protection 


INTRODUCTION 


In the world's merchant fleet diesel engines still have pre- 
vailed as the prime movers of ships. And, about 7096 of ship 
main diesel engines are low-speed ones. The general reason is 
their high efficiency even at partial load, and additionally — 
— a low cost of fuel oil used for low-speed engines (so-called 
heavy oil). 

Gas turbines do not substantially contribute either in total 
number or total output of ship main engines. For recent years 
neither number of gas turbines nor their total output has incre- 
ased significantly [6]. It's worth adding, for instance, that the 
total output of the low-speed diesel engines adjusted for using 
heavy oil, ordered in 2002, dropped by 48% relative to the 
records of the precending year. This is an important circum- 
stance from the point of view of development of gas turbine 
application as ship prime mover because mainly the above 
mentioned advantages of diesel engine over gas turbine are 
brought up in economical comparisons. 

Low speed diesel engines working on heavy oil are first of 
all used to drive large ships. Such engines are disadvantageous 
for ecological reasons (large SO, pollution of the environment). 
However it were not those reasons which made the number of 
orders for them in 2002, lower. Rather general safety reasons 
and economical situation in different regions of the world be- 
longed to them [6]. 

This paper 15 devoted to a propulsion system of merchant 
ship intended for sailing in the Baltic sea zone. Hence such 
system is to satisfy the ecological requirements contained in 
the relevant international conventions for special zones [1, 12, 
14] to which the Baltic Sea belongs too. It draws attention to 
gas turbine as a main engine for such ship, because it satisfies 


the ecological requirements, moreover it has another advanta- 
ges (discussed in detail below). 


ANALYSIS OF ECOLOGICAL 
REQUIREMENTS 
FOR PROPULSION SYSTEMS 
OF SHIPS OPERATING 
IN THE BALTIC SEA ZONE 


In accordance with the provision [14] of the International 
Convention for the Prevention of Pollution from Ships (adop- 
ted in 1973, modified in 1978 and amended in 1997) [12] su- 
phur content in marine fuel oils is not to exceed 4,596. And, in 
the zones subject to control of sulphur dioxide emission (e.g. 
the Baltic Sea, ports, costal waters up to 200 NM from land, 
and other zones specified by IMO) one of the two following 
conditions is to be fulfilled : 


» sulphur content in fuel oil is not to exceed 1,596 
> sulphur dioxide emission is not to exceed 6,0 g/kWh. 


A reason for such control and limitation of the environ- 
menal pollution by sulphur oxides is the protection of human 
health and the environment. Sulphur dioxide penetrates the envi- 
ronment as a product of combustion of natural fuels used for 
energy production purposes. 80% of natural fuels 15 combu- 
sted on land, and 2096 on ships. Most of sulphur dioxide rele- 
ased from ships pollutes the marine environment. Sea pollu- 
tion due to sulphur dioxide emitted from ships does not parti- 
cipate much in its total content in the marine environment if to 
compare the total amount of the sulphur contained in sea water 
(over 1099) with that contained in natural fuel resources 
(101%), However despite that, in some zones of the marine envi- 
ronment, including the Baltic Sea, the amount of sulphur oxi- 
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des emitted from ship propulsion systems is a threat to the envi- 
ronment. 

Among other, the Baltic Sea environment is polluted to an 
extent greater than average. The reason of the susceptibility of 
the Baltic Sea is its rather small water depth and salinity. 

Some reduction of the environmental pollution by sulphur 
oxides emitted from ships may be obtained in the following 
ways [4] : 


> by lowering sulphur content in fuel 
> by increasing propulsion efficiency 
> by separating sulphur oxides from exhaust gases 
2 by applying new energy sources, 
other than those have been using so far. 


However until now, application of an energy source other 
than fuel oil for ship propulsion has appeared either irrational 
(e.g. nuclear energy) for economical reasons, or impractical 
(e.g. wind). 

A source of sulphur dioxide in exhaust gas is the sulphur 
contained in fuel oil. Substitution of a fuel oil of high sulphur 
content for that of lower one or natural gas would result in 
a lower sulphur dioxide content in exhaust gas. 

However in this case some limitations associated with ava- 
ilability and cost of a fuel oil of low sulphur content, could 
appear. In present, the sulphur weight content in marine fuel 
oil amounts to 2.8%--3.0%. The situtation has not changed for 
the last ten years. The reducing of the sulphur dioxide content 
in exhaust gas by 50% , postulated by IMO, would require to 
lower halfas much the present sulphur content in fuel oil. Cost 
of such fuel oil would reach even 50 USD/t. 

For the zones of a more limited sulphur dioxide content in 
exhaust gas, such as the Baltic Sea, the limitation of the sul- 
phur content in fuel oil down to 196 1s considered. Cost of such 
fuel oil would increase by 58-85 USD/t. 

Due to efforts resulting from competition between ship 
engine producers diesel engine efficiency of both two-stroke 
and four-stroke engines has reached a high value (i.e. a low 
value of specific fuel oil consumption). To attain any higher 
efficiency would be very difficult. Hence the increasing of en- 
gine efficiency has a small influence on the reducing of sul- 
phur dioxide content in exhaust gas. 

In marine conditions the washing purification method of 
exhaust gas from sulphur dioxide by making use of absorbing 
features of sea water (respective to sulphur oxides) is conside- 
red. However this is a sophisticated and expensive system. 
Moreover, in such case fuel oil consumption increases and – in 
consequence – also the environmental pollution due to carbon 
dioxide. Moreover, its effectiveness depends on salt content in 
sea water, which is low in the case of the Baltic Sea. 

In association with the above given comments it seems that 
some lowering of sulphur content in marine fuel oil may be 
expected. Trends to lower sulphur dioxide content in exhaust 
gas emitted from ships would probably lead to elimination of 
heavy oil application for ship propulsion systems. This seems 
unavoidable especially in such zones as the Baltic Sea. This 
way a higher share of diesel oil in total fuel oil consumption 
for ship propulsion would cause that the gas turbine would 
become competitive against the high-pressure combustion pis- 
ton engine, also from the point of view of specific fuel oil con- 
sumption cost, especially if applied in a compound gas- 
-steam system (see Tab.1). 

The presented (Tab.1) specific fuel consumption of a gas 
turbine working in a compound gas-steam system is a little 
higher than that of a diesel engine. However as distinct from 
the gas turbine, the diesel engine consumes a rather large amo- 
unt of lubricating oil, which results in that the specific power 
cost of both kinds of propulsion becomes comparable. 


74 POLISH MARITIME RESEARCH, Special issue, 2004 


Tab.1. Sulphur content in marine fuel oil 


Kind of ship Sulphur 
propulsion Fuel content 
system 
Low-speed 
diesel engine 


Specific fuel 
consumption 
[kg/MWh] 


Sulphur 
emission 


liquid : 
heavy oil 
diesel oil 1. 
liquid : 
Gas turbine | diesel oil 
natural gas 


Compound liquid : 
gas-steam diesel oil 
system natural gas 0 


COMPARISON 
OF TWO SHIP PROPULSION SYSTEMS : 
BASED ON DIESEL ENGINE 
AND GAS TURBINE, RESPECTIVE 
TO TECHNICAL, ECONOMICAL 
AND ECOLOGICAL ASPECTS 


Gas turbine used as a ship main engine 
has many advantageous features, namely : 


e 


% 


low specific initial cost 

small dimensions 

low weight 

modular structure 

low level of vibrations and noise 

long time between successive overhauls. 


e 


ЫХ 


e 


= 


ЫХ 


e 


= 


X 


% 


Despite that its share as а ship prime mover іп total output 
and number of applied engines on ships is still low, because 
gas turbine has some significant drawbacks in comparison with 
high-pressure combustion engine; these are : 


$ low efficiency when working in a simple 
propulsion system 

€ fast dropping efficiency at decreasing load 

it requires light fuel oil to be applied, which 
is more expensive than heavy oil used for low- 
-speed diesel engines (covering more than 
70% of total power used for ship propulsion). 


For this reason the intensive research ordered by US Navy 
is carried out , aimed at development of a new type of marine 
gas turbine working on low-cost fuel oil [26, 27, 28]. 

The gas turbine applied on ships are mainly of an aircraft 
type (about 95%). Their output is in the order of 3000+27000 
kW, initial gas temperature : 850+1100°C, compression ratio: 
10+26, efficiency: 27+40%, specific fuel oil consumption: 
220+260 g/kWh [13]. Most of them is delivered by General 
Electric, among which LM 2500 turbine is most popular. 

The investigations in question are aimed first of all at in- 
creasing initial temperature of combusted gas in order to make 
combusting a lower quality fuel oil possible. Moreover, waste 
heat recovery from exhaust gas is also planned to lower fuel 
oil consumption down to 170+180 g/kWh. The application of 
exhaust gas recuperation at outlet from turbine, and of intersta- 
ge compressor cooling makes it possible to decrease specific 
fuel oil consumption by about 30%. The investigations are also 
aimed at obtaining a lower fuel oil consumption at partial load 
levels of gas turbine. 

TheRolls-Royce, a British firm in cooperation with the 
Westinghouse, a US firm has built WR-21 ship gas turbine of 
a new type in which the interstage cooling and recuperation 
system (ICR) have been applied. It is a triple-shaft turbine of 
25000 kW output and a modular structure. 
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General Electric has modernized its known LM2500 turbi- 
ne (now marked : LM2500R) by adding an interstage cooling 
and recuperation system [26,28]. 

At outlet either from a turbine or a diesel engine a steam 
generator supplying a steam turbine can be installed to utilize 
exhaust gas heat. In the case of gas turbine, exhaust gas trans- 
fers to environment more heat than in the case of diesel engine 
in which a significant part of the heat is absorbed by the water 
cooling the engine . Therefore in a combined gas turbine/ste- 
am turbine system, the steam turbine's output is much greater 
than that in a combined diesel engine/steam turbine system. 


For ship propulsion two kinds 
of combined gas/steam systems are used : 


© in (һе COGAS system (COmbined Gas And Steam turbi- 
ne) ship's screw propeller is driven through a mechanical 
transmission gear 

9 in the COGES system (COmbined Gas and Steam turbine 
Electric) the so-called electrical transmission is used. The 
turbines drive electric generators which supply a common 
electrical power system (electrical network), and ship screw 
propeller(s) is (are) driven by electrical motors. 


The Deltamarin, a Finnish advising firm has compared the 
COGES propulsion systems and diesel-electric propulsion sys- 
tems (fitted with an electric transmission) installed on diffe- 
rent merchant ships. On this basis it has specified the follo- 
wing advantages of the COGES propulsion system [17] : 


* low initial cost 

ж high thermodynamical efficency: specific fuel oil consump- 
tion equal to or lower than that of a diesel engine of a com- 
parable output (working through an electrical transmission) 
low operational cost (servicing, maintenance, lubricating oil) 
low sulphur content in exhaust gas (exhaust gas purifica- 
tion devices are not necessary) 

low level of vibrations and noise (e.g. passenger cabins can 
be located in the vicinity of engine room) 

high availability (high reliability, and possible fast opera- 
tion of repair or replacement of engine) 

much smaller number of auxiliary devices, which makes 
initial and operational costs lower 

small engine room dimensions resulting in a greater space 
for cargo. 


+ ж жа жа HH 


Technical characteristics 


Out of the technical factors influencing choice of a kind 
of ship propulsion system the following are considered : 


Жж operational reliability 
The gas turbines show low failure frequency and high avai- 
lability greater than that of the diesel engines and steam 
turbines. Reliability of gas turbines used as ship prime mo- 
vers amounts to about 98% [9]; the result is based on multi- 
-year service experience on aircraft, ship and industrial tur- 
bine engines. 


æ overhaul labour consumption 
Due to a modular structure of gas turbine, access to its par- 
ticular units and replacement of failed parts is very fast. It 
is possible to overhaul the engine without taking the ship 
out of operation. Gas turbine engine may be replaced wi- 
thin less than 24 h. Its small weight and gabarites make it 
possible to handle the whole engine at once, if necessary. 


** power transmission 
In the case of fully mechanical powering, the gas turbine 
propulsion system requires a mechanical transmission gear 
to be applied. It is also required in the case of using high- 


and medium-speed diesel engines as ship prime movers, 
however low-speed diesel engines do not require it. 

As applications of the electrical transmission are still gro- 
wing the way of power transmission from both gas turbine 
and diesel engine becomes analogous. 


ж manoeuvrability 

Diesel engines are reversible; gas turbines are irreversible 
as a rule. The stage of controllable guide vanes which ma- 
kes reversible action of gas turbine possible, is rarely ap- 
plied. The disadvantage of gas turbine can be overcome by 
using reversible transmission gears or controllable pitch 
propellers. Small gabarites of gas turbine rotor result in its 
very small inertia. Due to this, manoeuvrability of gas tur- 
bine ships is very high. Starting operation of gas turbine is 
very fast; starting from cold state, it is capable of develo- 
ping full output within several minutes. 


з automation of operation 
In comparison with diesel engines, gas turbines can be au- 
tomatically controlled in a more easy way. It leads to possi- 
ble reduction of a number of operators. 


** vibrations, noise 

Vibrations generated by gas turbine in operation are much 
smaller than those due to diesel engine. They are easier 
damped due to their higher frequencies than those genera- 
ted by diesel engines. Noise intensity due to gas turbine 
operation is much lower than that generated by diesel engi- 
ne. Moreover, gas turbines are usually covered by an aco- 
ustic insulation. 


Economical characteristics 


Economical characteristics of the considered kinds of ship 
propusion system were compared by using a selected ship as 
an example. 

It was the fast Ro2000 ship of 10.000 DWT load-carrying 
capacity, built by German Flender shipyard [2,3]. Its two con- 
trollable pitch propellers driven by two medium-speed, four- 
-stroke diesel engines of 50 MW total output , propell the ship 
up to 28 kn speed. The ship's propulsion system was compared 
with a system ofthe same power but driven by a gas turbine, as 
well as with the COGES, a combined gas/steam system. 


As far as the diesel engine propulsion system 
is concerned the following information was assumed valid : 


* Two MAN-B&W 12V 48/60 medium-speed diesel engi- 
nes, each of 12600 kW output at 600 rpm, were assumed to 
be the ship's prime mover 

* Electrical power was produced by four generating sets of 
1000 kW output each, driven by four HFC5 632/14K- 
-T260L-EX Yanmar diesel engines 

* It was assumed that four waste heat boilers delivered each 
of 1000 t/h of steam 

* The amounts of consumed electric and heat energy were 
estimated on the basis of available data dealing with a ro-ro 
ship of a similar load-carrying capacity. 


As far as the gas turbine propulsion system 
Is concerned it was assumed that : 


* the system consisted of two WR-21Northrop Grunman / 
Rollce-Royce gas turbines of 25000 kW output each, dri- 
ving two ship screw propellers 

* the turbines are equipped with an interstage cooling and 
recuperation systems due to which the turbine's efficiency 
amounts to 4096 

* Exhaust gas from the recuperator supplies a waste heat 
boiler. 


POLISH MARITIME RESEARCH, Special issue, 2004 75 


An analysis of characteristics of ship gas turbine propulsion system 


* Specific fuel oil consumption of the turbines not much va- 
ries within 50+100% interval of the rated output. 


As far as the COGES system is concerned it was assumed 
that two LM2500 General Electric gas turbines of 22400 kW 
output each, as well as a steam turbine of 10000 kW output, 
was applied. 

In order to compare the above specified propulsion 
systems the following measures were used : 


propulsion efficiency : 
Р, =, -3600 
| Во Wa 
energy efficiency : 
‚_ (P, N +Р,):3600+0 
ДЕ Bh Wa 
specific fuel oil consumption : 
B, 
Ens 


where : 


bs 


Ву - total fuel oil consumption per hour 
(by main engine, auxiliary engines 
and boilers) 

P, - main engine output 

Ра - generated electric power 

Q - hourly flow rate of produced heat 

W,- net calorific value of fuel oil 

Ng - transmission gear efficiency. 


Appropriate values of the economical characteristics ofthe 
ship propulsion systems in question are compared in Tab.2. 


Tab. 2. Comparison of the economical characteristics 
of the considered ship propulsion systems 


WR-21 COGES 
gas turbine} system 


и MAN-B&W 
elements 
diesel engine 


propulsion efficency, 
[20] 


specific diesel fuel 
oil consumption, 194 
[g/kWh] 
energy efficency, [%] 
specific lubricating 
oil consumption, 
[g/kWh] 


specific cost 
of marine diesel oil, 
USCent/kWh 
specific cost 
of lubricating oil, 
[USCent/k Wh] 


price of heavy oil, 
[USD/t] 
specific cost of heavy oil, 
USCent/kWh 
It was assumed that both the diesel engine and gas turbine 


is supplied with the same kind of fuel oil (marine diesel oil). 
In such conditions the specific fuel oil consumption of the 
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COGES system appears the smallest, somewhat smaller than 
that of the medium-speed diesel engine in question. However, 
the diesel engine also consumes lubricating oil (much more 
than the combined system COGES), which makes that the spe- 
cific power cost of the COGES system amounts to about 0.88 
of that of the diesel engine propulsion system (see Tab.2). 

The specific power cost of the gas turbine propulsion sys- 
tem (of a simple configuration) is only a little greater and it 
amounts to about 1.03 of that of the diesel engine propulsion 
system. 

The specific cost of heavy oil is given in the last row of 
Tab.2. On its basis it can be stated that the specific power cost 
of the COGES system amounts to about 1.4 of that of the low- 
-speed diesel engine propulsion system. 

As a matter of fact at partial load diesel engine's efficiency 
decreases much slower than that of gas turbine. At the half- 
-load, diesel engine's efficiency drops e.g. by 5% relative to its 
value at the rated load, whereas the efficiency of LM2500 gas 
turbine drops by 10%. However, in the COGES system 50% 
of the rated load can be obtained by keeping in operation at its 
rated load only one of the two gas turbines. 


Ecological characteristics 


Exhaust gas from diesel engines contains more not only 
sulphur oxides but also nitrogen oxides due to a greater effec- 
tiveness of fuel combustion processes, in consequence — a gre- 
ater temperature. The exhaust gas also contains more carbon 
oxide and hydrocarbons. Tab.3 enables to compare in this re- 
spect the considered ship propulsion systems. 


Tab. 3. Comparison of the ecological characteristics 


of the considered ship propulsion systems 


propulsion 


diesel engine 


gas turbine — fuel oil 
simple natural gas 
configuration 
gas turbine — | fuel oil 1. Fi 
COGES system natural gas 


The necessity of complying with the environmental pro- 
tection rules entails undertakings dealing with construction of 
the diesel engine itself, as well as with the exhaust gas purifi- 
cation devices [31]. Such undertakings are costly. However the 
gas turbine propulsion system can satisfy the environmental 
protection rules without any additional devices. 


COMPARISON OF SHIP POWER 
PLANT CHARACTERISTICS 
RESPECTIVE TO A KIND 
OF SHIP PROPULSION SYSTEM 


A kind of ship propulsion system (a kind of prime mo- 
ver) influences such characteristics of the ship power plant 
as : its weight, volume of occupied space , number of opera- 
tors. They play an important role from the economical point 
of view. By decreasing the volume of its space it is possible 
to increase amount of cargo or number of passangers to be 
shipped. The lower the weight of the power plant, the smal- 
ler the power demanded for ship propulsion (the smaller fuel 
consumption). A smaller number of operators leads to a lo- 
wer personnel cost. 
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Comparison ofthe ship power plant respective to its weight 
and volume with taking into account a kind of its propulsion 
system, can be performed by using two indices : the so-called 
specific mass and power density (concentration). The specific 
mass 15 defined as the ratio of the total power plant mass (of 
main engine, auxiliary devices, and empty piping systems) and 
the engine's rated output. The power density is defined as the 
ratio of the engine's rated output and the engine room volume. 
Values of the indices, presented in Tab.4, make it possible to 
compare the example power plant with the respect to the three 
considered kinds of propulsion system. 


Tab. 4. Specific mass and power density of the example ship power plant 
in three considered versions of propulsion system 


Considered Пена propulsion system 
pose on: 


12V 48/60 WR-21 COGES 
diesel engine | gas turbine system 


Index 
че c mass 
B ru a | је 


(concentration) 
[kW/m'] 

Out of the considered kinds of ship propulsion system, the 
diesel engine power plant is ofthe greatest specific mass, abo- 
ut six times greater than the gas turbine one, and over three 
times greater than that of COGES type. The power density (con- 
centration) of the gas turbine power plant is over five times 
greater than that of the diesel engine power plant; the greatest 
power density has the COGES version of the power plant - 
over sixteen times greater than that of the diesel engine power 
plant, and over three times greater than that of the gas turbine 
power plant (see Tab.4.). 

As it has been already mentioned , the gas turbine power 
plants (including their COGES versions) are suitable for auto- 
mation and control. This way, decreasing the number of power 
plant operators and simultaneous decreasing volume and mass 
of power plant as well as lowering its operational cost, can be 
possible. 

Growing interest for application of electrical drive to ship 
propellers (by using the electrical transmission) also directs grea- 
ter and greater attention to gas turbines. In this respect it is worth 
adding that in such version of propulsion system gas turbines 
can be located on the deck, which eliminates a troublesome 
arrangement of compressor's air supply pipes delivering air from 
the deck to the power plant. 


SPECIFIC ASPECTS OF GAS TURBINE 
APPLICATION ON SELECTED SHIPS 


Choice of a main engine location depends on a type of ship 
and applied engine, and a power transmission way from the 
engine to propeller. For instance application of an electrical 
drive to ship propeller provides the designer with a great fre- 
edom in locating the main engine relative to a mechanical trans- 
mission of power. Any propelling engine due to its mass, thrust, 
torque, transient state loads and sea state must have a suitable 
foundation firmly connected with ship's hull structure. Such 
propelling engine as gas turbine is a source of high temperatu- 
re and noise, hence protective means should be applied to lo- 
wer noise and high temperature. The features have an influen- 
ce on location of main engine within engine room as well as on 
construction of its foundation. In ship power plants based on 
gas turbines the main engine is usually installed in an individu- 
al casing seperated from ship's hull structure. The gas turbine 
module is connected to air intake and exhaust gas systems. 


Such solution of the engine's casing (container-like) enhances 
also fire-safety of the whole power plant. 

Location of engine's seating depends on many aspects. One 
ofthem is troublesome centering the engine and propeller sha- 
fts. To avoid such operation a non-mechanical way of power 
transmission from the engine to propeller, usually an electrical 
transmission system, 1s applied. In this case different variants 
ofthe engine's location within ship's engine room are possible. 
The main engine can be located either aft or fore; if just aft 
then it is possible to make the shaft between the engine and 
propeller as short as possible , which reduces engine room's 
volume. The electrical transmission makes it possible to locate 
the main engine fore, without applying any long shaft. In most 
cases ship's power plant is usually located within ship hull's 
space. 

Along with application of gas turbines new possibilities to 
locate the main engine even on the ship's deck, have appeared. 
Small mass of gas turbine provides the designer with a greater 
freedom in locating the power plant. For instance it is allowed 
to locate ship power plant above the main deck both on large 
and small ships, also in the case of application of heavy indu- 
strial gas turbines. 

In Fig.1 is presented a variant of location of a power plant 
equipped with a heavy gas turbine installed on the main deck 
ofa tanker. In the power plant a gas turbine with heat recupera- 
tion system, as well as electrical drive of propeller, are applied. 
Mass of propeller, shaft, electrical motor and gas turbine deter- 
mines mass centre location below the main deck level (engine 
- generator - electrical motor - turbine shaft). 

Dimensions of air and gas ducts of gas turbine influence 
location of the gas turbine itself. If the gas turbine is located 
high then lengths of inlet and outlet pipes are shorter, hence 
mass and dimensions of air and gas ducts are also changed. 


Fig.l. An example arrangement of tanker's power plant fitted 
with gas turbine and electrical drive of propeller 


In Fig.2 a gas turbine power plant installed on a container- 
-ship is presented. Location of the gas turbine with compressor 
aft makes it possible to shorten inlet piping as much as possi- 
ble.In such location of the gas turbine the reduction gear is pla- 
ced before the engine, and the propeller shaft is led aft below 
the engine. This way the length ofthe power plant appears shor- 
ter in comparison with those of other types of heat engines. 

Application of gas turbines to ship propulsion makes free 
access to all devices possible, both during service and over- 
hauls. Moreover, lengths of air and gas ducts can be minimi- 
zed. The number of devices associated with auxiliary mecha- 
nisms is reduced, and the engines themselves do not require to 
be additionally cooled with sea water (hence sea water pumps, 
their piping systems etc are not necessary). 


POLISH MARITIME RESEARCH, Special issue, 2004 77 


An analysis of characteristics of ship gas turbine propulsion system 


Fig. 2. An example solution of gas turbine power plant of a containership 


SUMMARY 


О Application of gas turbines on merchant ships is still limi- 
ted [6] due to their high operational cost (specific fuel oil 
consumption cost) especially in comparison with that of 
low-speed diesel engines working on heavy oil. However 
the ecological requirements established by the Internatio- 
nal Maritime Organization (IMO), limiting content of sul- 
phur oxides in exhaust gas from ship propulsion engines 
[1,12,14], can make it possible — or even force — a broader 
application of gas turbines as ship prime movers. For spe- 
cial zones to which the Baltic Sea belongs too, the require- 
ments are more stringent. 


O Fulfilment of the requirements by diesel engines could be 
achieved by : 
c» using a fuel oil of an appropriately low sulphur content 
(it would mean to stop using heavy oil) 
c» installing exhaust gas purifying devices. 
Both the ways of reducing content of sulphur oxides 
in exhaust gas are expensive. 


O Application of gas turbine to ship propulsion does not re- 
quire exhaust gas to be purified, however it requires a fuel 
of low sulphur content to be used. In the case if the ecologi- 
cal requirements force application of fuel oil of a similar 
quality for diesel engines , then the gas turbine propulsion 
system will be comparable with the diesel engine one also 
respective to economical aspects (specific fuel oil consump- 
tion cost). In this respect it would be even more favourable 
if used in a combined gas turbine /steam turbine system, 
especially at common production of electrical and heat ener- 
gy (COGES system). 


O Gas turbine propulsion system has many technical advan- 
tages which make gas turbine superior , in this respect, over 
diesel engine. Diesel engines, especially low-speed ones 
working on heavy oil, are superior over gas turbines only 
in the respect of specific fuel oil consumption cost. This is 
why their share in the total propulsion power installed in 
merchant ships amounts to about 70%. 
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O In the Baltic Sea zone the gas turbines are expected to be 
more and more used for propulsion systems of such ships as : 
c» fast car-passenger ferries 
c» fast cargo ships 
= special vehicles (hydrofoils, hovercraft, motor yachts). 
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A proposal of ship gas-turbine driven waterjet 
propulsion — preliminary considerations 


Wiestaw Prochnicki, Assoc.Prof.,D.Sc. 


ABSTRACT 


In the paper are presented preliminary considerations concerning the efficiency of water- 
Jet ship propulsion system, as well as the calculation of main dimensions of waterjet chan- 
nel. The friction and momentum losses of the flow channel have been roughly estimated by 
using Fliegner 5 equations. An important conclusion is confirmed that the summary losses 
are inversely proportional to square of ship velocity (~1/17). On the other hand the ship 
propulsion power is directly proportional to third power of ship velocity (~из). Therefore 
to minimize ships hull resistance, hulls of waterjet-driven ships ought to be of a great 


slenderness — e.g. L/B — 15, stabilized by sponsons, or of semi-swath hydrofoil-supported construction. 


Keywords : ship propulsion, hydromechanics, waterjet 


INTRODUCTION 


Attempts to apply waterjet propulsion to ships have already 
appeared in the ancient ages?. However the Hamilton's low- 
-power waterjets (1954) intended for propulsion of fast river 
boats can be considered as the first developed construction), 
Nowadays many firms build waterjets of a wide range output : 
from a few hundred kW to a few dozen MW. 

Rolls-Royce, a leader of modern waterjet constructions, 
offers complete propulsion systems for very fast ferry ships, 
cargo vessels, warships, cruise ships and yachts. 

Up to now about 1400 assemblies ofthe output from 220 to 
70 800 kW have been built. And, a ship propulsion system of 
250 MW output is under development 


An example of Rolls-Royce waterjet propulsion system 
for a passenger ferryship is shown in Fig.1. 


Fig. 1. ROLLS-ROYCE waterjet propulsion system for a passenger ferryship 


CONSTRUCTIONAL ELEMENTS 
OF WATERJET PROPULSOR 


In Fig.2 the main functional elements 
of waterjet propulsor are diagramatically presented. 


oil ond Guide vanes 
Pump (diagonal, helicoidal, propeller) 


1. $ 
u 
===> 


№ un 
m 


Diffuser 


Fig. 2. Main elements of waterjet propulsor. 
Notation : w, — angular velocity of pump; wj, w2, Wp — inflow, outflow, and 
pump flow velocity, respectively; u — ship speed; H — water elevation 


Labyrinth glands’ Flexible coupling 
Em uu = 


The flow channel consists of a diffuser, rotordynamic 
pump, guide vanes (stationary) and confusor. The channel 
surface should be of high smoothness (low coefficient of rela- 
tive roughness of surface). 

The ship manoeuvring devices 
are shown diagramatically in Fig.3. 


b) 


ES 


STEERING 
STARBOARD 


NO SPEED 


Fig. 3. Schematic diagram of ship manoeuvring elements: 
a) for ship reversing; b) for ship turning 
Notation : a - steering blades, Q - blade turn angle 


The effectiveness of ship steering is very high because the 
waterjet force Ғ, for u—0 is Ғ.(о)->Е, max, and the manoeu- 
vring operations are usually executed at low ship's speed : ио. 


ANALYSIS OF WATERJET 
PROPULSION EFFICIENCY 


By taking into account that the relation : 
w =u (1) 
is satisfied the absolute values of water velocities 
at channel inflow and outflow, are : 
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c = № и = 0 Q) 


The force of water stream reaction, is : 


Е = m(c; – с) = m(w, -u) (3) 
where : m [kg/s] — rate of water outflow. 
Foru-o 
E, (0) = Е, max € m: w5 (4) 
Therefore the ship propulsion power is : 
N, =F, u=m(w,-u)-u (5) 


and carry — over loss output : 
2 2 
с W^5,-—Uu 
N, =т- = та и 
2 


Applying the notation : 


(6) 


М, = pump propulsion power 
Np — overall pump efficiency 
Мр – power of summary flow 
channel losses 
one obtains : 
Мр пре М +N. +N, (7) 


The propulsion efficiency ny, is defined as follows : 


n ERE (8) 
Ма 
N, 
hence : 
N,:m 
u Тр 
c (9) 
қ. NAN, +N; 
If it is assumed that the efficiency 1} = 1, Шеп: 
=n = Nu 10 
Nn, m= =n, N, +N, +N, (10) 
Inserting (5) and (6) in (10) one obtains 
2 - 20 -u):u _ 21—-и)-п (1D 
Ке =. 


+——; 1- 1-2 М. 
ту т.о 


0 rus 
where : и = — - velocity index. 
Wo 
The expression of relative (nondimensional) losses power : 


NL 


2 
ти 


2 


which may be represented as a sum of n particular losses 
of the waterjet flow in channel, yields the following [5] : 


N n 
2—1-- 5L (12) 
Fig.4 shows the diagram of efficiency TIN. in function 

of the velocity index : 


u =— 
Wo 
and of the summary losses УС. 
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oK .91 02 03 04 ^w; 07 08 09 
Fig. 4. Propulsion efficiency NN, (for hp = 1) 
as a function of = = апа Уб, see [3] 
2 


The maximum (optimum) efficiency is the value : 


о — до /— 
Nn, ymax ` NN, (Woot) 


where : 
-1 


Уб | = ла — 09 


i-l 


In Fig.5 the theoretically determined summary losses 
(of friction and momentum) are shown : 


N n 
E = 
2—L- v 
ти i=l 
as a function of : 


n 
У С = НА) [5] (14) 
= where : 

Х(Ее)- the friction number (dependent on the Rey- 
nolds number Re, and surface roughness co- 
efficient of the flow channel), are shown for 
0 = 0.8 and Re-10/. 


LU Е c] LÀ 
ж E 
10 u-0.8 
0, Е ‘a 
? id 
> = 
= = ] 
ei 3 
3 
2 
‚Г [ 1.0 
0.5 z 
H=0 
0.1 E = 
= === 5 
TI [E 
6 pen 
== еы ы 
0.02 шн 


| 
2 4 6 8 10 12 14 16 18 20 22 24 26 28 30 
u [m/s] 


Fig. 5. Coefficient of summary losses (5245 ) for ü —0.8 and Ке=107 
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In Fig.6 the efficiency my max is presented in function 
of the summary losses | Yi, 6) 


0.9 
0.8 T4 
o H ЫШ 
n Nu , тах + | 
0.6. | 
T + 

0.5 HHH 

0.4 an Mi 4 BED 

0.3 | 
001] 2 3456789] 2 3 4567810 2 3 4567810 
Fig. 6. The optimum efficiency TIN, max in function ој (X. Ci) 


MAIN DIMENSIONS 
OF WATERJET PROPULSION SYSTEM 


The power of ship propulsion is : 


N, =m-(w,-u)-u = 5402.16:0-0? (15) 


where : 
40° [HP 
UN T Fa - power coefficient 


Cy - Admiralty coefficient 

Dy [t] = рУ; - ship mass 

Pw = 1,021.03 [t/m' ] - water density 
У, [m'] - ship displacement 

u [m/s] - ship speed. 


Inserting the equation of waterjet flow rate : 


m-2-:A,:w;:p, P (16) 
into (15) one obtains : 
2 
каа Rs p к= 

2А, | W3 W3 

hence : 

u? 

А, 22.650. —— (7) 


The graph of the waterjet outflow area of one of both nozzles 
А, =A,(U, о) is shown in Fig.7. 


Dimensional proportions of the outflow nozzles 
are shown in Fig.8 [6]. 


AN EXAMPLE ARRANGEMENT 
OF THE SHIP WATERJET PROPULSION 
SYSTEM WITH GAS-TURBINE DRIVE 


The considered waterjet system is intended 
for a passenger ship of the following parameters : 


ship length L = 80 m, L/B = 6.3 , ship mass Dw = 4000 t, speed 

u — 6 m/s (11.7 kt.), Admiralty coefficient С, = 300, power 
output № = 980 kW, having main engine : one-cylinder gas 
turbine with heat exchange regeneration (air preheater), of ef- 
fective efficiency : 


no! = 0.33 +0.34 


A» [m?] 


| | 
0.5 06 u=u/w2 0.7 0.8 
Fig.7. Waterjet outflow area of one of both nozzles A» in function 


u 
of the power coefficient ct and velocity index u = — 


€) – racing boat of u = 25 m/s 


Assuming the velocity index u = = 2:00:75 


2 
and taking into account the power coefficient 


a= ES = 0.84 


one gets from the diagram in Fig.7 : Аз = 5 m? 
hence, the outflow diameter is: d? = 1.8 m. 


i=l 


If H - 0 , = 0.9 and Р: | 0.085 


is assumed (Fig.5), then : 


d EL ефт 


n 


1+ Уб 


i=l 


hence : 


NW, =NN, ‘Np = 0.696 
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Fig.8. The stern area of various types of ships (maintaining the dimensional 
proportions): a) – passenger liner of L = 280 т, и 
b) — passenger vessel ој L = 88 m, и = 6 m/s, а = 0.907 HP/kt2; 
c) – passenger ship of L = 150 m, и = 9 m/s, а = 2.23 ИРМА; 
d) – tug of L = 40 m, и = 5.6 m/s, а = 0.207 HP/kt.^; 


12 m/s, а: = 5.01 НР/М; 


8l 
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The effective power of the gas turbine is : 


N9" - Nu 21400 kW 


1х 


u 


The waterjet propulsion system with gas turbine, 
located in the stern part of the ship, is shown in Fig.9. 


Fig.9. The waterjet propulsion system located in the ship s stern part: 
1 - diffuser; 2 - confusor (outflow nozzle); 3 - guide vanes; 4 - thrust & 
journal bearing; 5 — diagonal-helical pump; 6 - labyrinth gland/journal 
bearing; 7 - highly flexible coupling; 8 - gas turbine; 9 - combustion 


chamber; 10 - gear box; 11 - axial-radial compressor; 12 - recovery heat 


exchanger (air preheater); 13 - exhaust gas duct; 14 - air inflow duct 


The driving gas turbine with regeneration heat exchanger, 
is shown diagramatically in Fig.10. 


і 5m »- 


Fig.10. Longitudinal cross section of a gas turbine 
with regeneration heat exchanger (suitably shaped and located) 


Ifthe ship's hull is very slender, e.g. of L/B = 15, and stabi- 
lized by two pairs of the slim sponsons (s) — see Fig.11, the 


ship will develop the speed up to 8 m/s (15.55 kt.), at the same 
output of main engine. 


0 | 20 40 60 


Fig.11. Profile of a very slender ship (L/B = 15) 
propelled by a gas turbine waterjet system 


FINAL REMARKS 


2 The waterjet propulsion system has some essential advan- 
tages in comparison with that conventional based on screw 
propellers : 

absence of any rotational elements behind the hull 

smaller losses due to wake current (of about 13% = 

= 19% [2]) 

favourable efficiency index at partial loads 

torsional vibration diminution 

effective ship turning and reversing 

possible elimination of rudder blade, resulting in a lo- 

wer hull resistance. 


995 22 
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> Summary losses are inversely proportional to из, and mo- 
notonically decreasing with the ship speed u increasing. 
They are very weakly dependent on the velocity index u. 
> The great dimensions of flow channel are a serious pro- 
Мет of waterjet arrangement. The outflow area A, is func- 
tionally dependent as follows : 
—2 
А, со a( ы —) 
]—H 
where : о, acc. to eq.(15), is a coefficient dependent on dy- 
namic ship hull properties, whereas the expression ar 
is directly proportional to the propulsion efficiency, which 
means that the greater the propulsion efficiency, the greater 
the nozzle dimensions (!). 
2 The following problems are important for improving con- 
struction of waterjets : 
c» smoothness maintenance 
c» the pump and flow channel protection against acciden- 
tally flowing-in hard pieces 
c» development of automatic control elements for ship tur- 
ning and reversing 
c» development of construction of journal & thrust bearing 
situated inside pump, e.g. a ceramic bearing with wa- 
ter-film lubrication. 
NOMENCLATURE 
B - ship breadth 
ср, еә - relative water speed at inflow and outflow, respectively 
d, - diameter of pump outflow 
4, - diameter of confusor outflow 
Fy - propulsion force 
H - waterjet height in relation to sea level 
HP - horsepower 
kt - knots 
L - ship length 
m - flow rate 
№ - carry-over loss output 
N, - ship propulsion power 
u - ship speed 
W,,W>2 - inflow and outflow water velocity, respectively 
& - nondimensional loss of i-th component 
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a) Archimedes axial water pump (287-212 AChN), Leonardo 
da Vinci pump (1452-1519), Toogwood's & Hayese's patent 
(1661), Beniamin Franklin's pulsators (1706-1790) 

b) So-called „First Hamiltonians” 

€) The pump may be of various type: diagonal, helicoidal 
or propeller pump 

4) Water flow rate of two waterjet nozzles 


Automatic control systems for ships fitted 
with podded propulsion drive (POD) 


Czestaw Dymarski, Assoc.Prof., D.Sc. 
Karol Dabrowski, M.Sc., Eng. 


ABSTRACT 


The paper shows purposefulness and possiblity of automation of ship propulsion systems , especially those 
of POD type, intended to be used on four Baltic navigation ships : a containership, ro-ro ship, product 
tanker, as well as river-sea-going ship, designed within the frame of the Eureka ,, Baltecologicalship " 
project. Requirements of ship classification societies, and general ones for modern ship control systems 
were presented, as well as a review and analysis of currently applied power plant automation systems for 
diesel- electric ships fitted with podded propulsion drive (POD), were performed. Moreover real ways and 
possibilities of controlling ship s course and speed were indicated. Information contained in this paper may 
facilitate selecting appropriate design assumptions for a propulsion control system to be applied on the 
above mantioned ships. 


Keywords : ship propulsion, azimuthal propulsion, ship control systems 


INTRODUCTION 


The outfitting of ships with automatic control systems has 
already been a standard for many years. Such systems can be 
found practically in any shipboard device. However develop- 
ment of computer and data transmission technology has ope- 
ned new possibilities for control systems. They are associated 
with joining control systems into large ones, which is a suces- 
sive phase of development of automatic control for ships. 

In a further prospect the development should lead to buil- 
ding ship power plants practically unmanned regarding their 
automation level. Unfortunately the goal is still very distant at 
the present development level, mainly due to too low reliabili- 
ty of mechanical devices. 

The podded propulsion drive (POD) is a step forward in 
this respect as instead of a single high-output main engine a set 
of relatively low-output engines driving electric generating sets 
can be applied. A failure of one of the engines does not cause 
the blackout of the entire propulsion system to occur. 

The other advantage of the POD system is a large number 
of electrical elements which can be easily automated.The de- 
vices are often fitted with its own microprocessor system and 
communication interface. This way vast possibilities have been 
opened for designers of automatic control systems. 

The presented paper concerns the automatic control sys- 
tems for four types of ships : a containership, ro-ro ship, pro- 
duct carrier, driven by azimuth propellers, as well as a river- 
-sea-going ship driven by screw-rudder propellers. It can be 
assumed that a ship type only to a small degree influences ove- 
rall structure of automation systems, therefore in this work no 
separate analysis for each of the ship types was performed. 


REQUIREMENTS 
OF SHIP CLASSIFICATION SOCIETIES 


The ship control systems are to provide high level of safe- 
ty. The assumption has to be realized by satisfying a number of 
requirements established in the rules of ship classification so- 
cieties. Below, the following requirements for ship control sys- 


tems are formulated on the basis of the rules and publications 
of Polish Register of Shipping [11,12] : 


2 ensuring control of functioning and technical state of all 
important devices, i.e. electric motors installed in PODs, 
frequency converters, diesel engines etc, as well as genera- 
ting appropriate signals and alarms to inform operators on 
occurrence of a malfunction 


2 making it possible to control ship's motion by influencing 
its propulsion system both from the bridge and other places 
on the ship, at which such need could appear 


У ensuring an automatic response of the system in the case of 
hazardous situations. 


ADDITIONAL REQUIREMENTS 
FOR CONTROL SYSTEMS 
OF POD-FITTED SHIPS 


The PODs differ substantially from the classical ship pro- 
pulsion system of a fixed or controllable pitch propeller and 
passive rudder system. Their essence consists in joining toge- 
ther propeller and rudder systems into the whole. In this case 
the propeller simultaneously fulfils the role of active rudder. 
Hence the automatic control system has to integrate both the 
function. Even if it is composed of two subsystems : steering 
and propelling, they will be expected to be functionally inte- 
grated into an entity. First of all they should be felt by the user 
as a single system. For this reason it is very important to appro- 
priately design and locate control heads. 

However the POD does not substitute all ship manoeuvring 
devices. For manoeuvring in ports bow thrusters are used. Ther- 
fore they are also provided on the designed ships in question. 
Because of high availability of electric energy on the ships it is 
recommended the thrusters to be electrically driven as it ma- 
kes it possible to connect their control systems with the POD 
control system. This way ship's crew can have at its disposal 
a complete ship control system both during sea voyage and 
port manoeuvring, integrated in a single panel which makes it 
possible to control all propulsion and manoeuvring devices. 
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Such solution makes it necessary to design a superior sys- 
tem responsible for control of information flow within the en- 
tire system. It may be a separate computer which, on the one 
hand, could control communication between control systems 
of particular devices : i.e. frequency converters, main engines, 
switching devices etc., on the other hand it could serve as 
a module making it possible to hand-over control from one to 
another control heads. It is assumed that direct control over 
a given device can be executed from one control head only, 
hence a system giving assent for such control is necessary in 
this case. 

The remaining requirements for a control system one can 
describe by dividing it into the following constructional ele- 
ments : 


» System's wiring — should ensure fast and safe data trans- 
mission. Simultaneously, number of communication cables 
should be limited. The first postulate could be best fulfilled 
by lightguide connections having rotatable optical connec- 
tors installed in the places where the cables pass from PODs 
to the deck.The best way to limit number of cables is to 
apply data buses (e.g. profi-BUS, an up-to-date and still 
modernized solution). This solution has been used for many 
years and it has proved to be correct. 


» Drive control and supervison system — should prevent 
against exceedance of limiting values in frequency conver- 
ters. In the case when the control system covers a greater 
number of devices (e.g. thrusters) several drive systems are 
involved. 


> Drive system — contains logic and drive supervision sys- 
tem. It analyzes operation of supply network, motor, trans- 
formers and remaining drive control elements. Its work is 
aimed at : 


A automatic switching-off in the case of exceedance of 
appropriate values in power supply circuit 

A safeguarding against overloads and shortings 

A checking the system before starting, generating availa- 

bility signals and stopping the drive in case of incorrect 

values of operational parameters 

control of starting and stopping the drive 

rotational speed control of driving motor 

communication with the remaining automation elements 

switching on/off of activity of control panels on ship. 


>>» 


Practically this is a management system 
which plays a decisive role in drive operating. 


> Data transmission control — ensures distribution and re- 
cording of the data obtained from transducers and measu- 
ring instruments, moreover it makes communication be- 
tween control heads and the drive control & supervision 
system, possible. 


> Operator panels, control heads — serve as basic commu- 
nication means between ship’s crew and automation sys- 
tem. It should be remembered that to apply a POD system 
the ship is be fitted with an electric power plant , associated 
switching devices and a transforming station. They are also 
equipped with an automatic control system of electric po- 
wer production and distribution. All the systems should be 
coupled together, and they should be made accessible from 
control heads grouped in one place. It is also crucial to de- 
sign interfaces of operator panels in such a way as to ensu- 
re easy operation ofthe system, without any additional cost 
of ship’s crew training. 


A separate problem is to ensure cooperation with an auto- 
pilot. It can be realized by making use of a data transmission 
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control system. The autopilot would be considered by such 
system as an additional control panel, without any possibility 
to control drive of thrusters. However the problem of a way of 
ship course-keeping remains to be solved. In the case of apply- 
ing the POD system it can be attained in two ways : 


1* — inthe same way as it is realized in the case of the conven- 
tional propulsion system ( screw propeller and passive 
rudder), i.e. by controlling direction of propeller thrust 
vector. The solution requires the entire POD to be rota- 
ted, but rotational speed of propellers can be maintained 
constant. 

274 _ the other way is to change rotational speed, hence by chan- 
ging the thrust delivered from the propeller of one POD 
only. Such action would produce a non-uniform propul- 
sion generating a ship course changing moment. Choice 
of one of the possibilities should be based on an effecti- 
veness analysis of both solutions and their energy effi- 
ciency. From the point of view of automatic autopilot 
system both the versions are rather easily realizable. 


ANALYSIS OF POSSIBLE APPLICATION 
OF DIFFERENT CONTROL SYSTEMS 


A system approach is necessary for designing a modern 
ship propulsion system fitted with electric power plant and 
switching station . An analysis of a ship as the whole , and not 
as a sum of its particular equipment elements, may help in lo- 
wering risk of malfunctioning, incorrect communication be- 
tween devices etc. Also, users of such system may be easier 
and faster trained in operating it. In the case in question conti- 
nuous development of user's interfaces (operator panels) is also 
of a great importance. It makes faster and more effective reac- 
ting by ship crew during ship operation, possible, as well as 
that special qualifications for operating automation devices by 
ship crew, are not required. 

The ship electric power plant supplies all driving devices 
and delivers electric energy to the ship electric network. As all 
supply systems are mutually connected the intergrating of au- 
tomation systems into a global system whose elements are ca- 
pable of adapting their work to an energy situation ofthe entire 
ship, becomes rational. 


Propeller thrust control 


Two kinds of ship propulsion control are possible: maintai- 
ning a set value of ship speed, and maintaining a set value of 
power. In the first case a speed setting will be translated into an 
appropriate number of propeller rotations. It will cause demand 
for engine output changeable. In the case of the engine output 
control a reference value is the main shaft power, and the shaft 
rotational speed should be so changed as to maintain the po- 
wer value constant. 

Setting an appropriate speed value on the bridge should 
initiate a power changing program which determines slope 
values of power rising and trailing edges in function of the set 
value and current load. It has to prevent the propulsion system 
from an overload at too sudden speed change. An example dia- 
gram showing increase of main shaft power, speed and torque 
Is presented in Fig.1. 


Speed change versus load | Power change versus load 


Slope of Slope of 
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Speed based propulsion ramp 


— Rotations per minute (rpm) 
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0 5 10 15 20 25 35 40 45 50 55 60 65 
Time [s] 
Fig. 1. An example diagram showing increasing course of main shaft power, speed and torque in function of time 
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Fig. 2. Schematic diagram of the ship speed control procedure 


As soon as a power level becomes steady a limiting value 
of propeller thrust is checked. The maximum set value of ship 
speed will be hence dependent on limitation of propeller thrust 
force. Next the system checks a value of engine shaft torque 
and adjusts its rotations in such a way as not to exceed the 
maximum shaft torque. After performing the procedure an ap- 
propriate value is sent to the engine. 

A sudden drop of rotational speed ofthe engine should cause 
the engine passing onto generating work mode and transfer- 
ring a part of its output to an inverter and network, and a power 
surplus to a resistor. 


Propeller thrust direction control 


Rotation of propeller's POD makes it possible to produce 
large forces transversally directed to ship's longitudinal axis. 
It may be hazardous as a large heeling moment appears in such 
situation. Apart from a large heel of ship, which itself 1s hazar- 
dous (possible shipping of water on deck) it generates large ac- 
celerations which may cause a shift of cargo and – in consequen- 
ce — trigger other hazards including loss of stability and ship’s 
capsize. Hence the control system is required not only to control 
a propeller thrust value depending on a value of POD’s rotation 
angle, but also to limit speed of its rotation. The speed limitation 
should not be dependent on an actual thrust value. At its low 
values the POD’s rotation can be executed rather fast. Greater 
thrust values should cause POD’s rotation speed lowering. 


Limitation of power level in the system 


A number of generators running in the system depends on 
an actual power input. Appropriate limits for power input after 
exceedance of which the number of running generators is chan- 
ged by one, are established. This way the energy production 
becomes adjusted to the power demand from the side of ship's 
network. 


Limitation of engine output 


In the case of sudden manoeuvres it may occur that engin's 
output or torque reaches its respective maximum value earlier 
than set speed value. In the situation the system should be ca- 
pable of stopping the rotational speed increasing. 

Besides, in the case of automatic switching - off the gene- 
rators the control system reduces engine output to such an extent 
as to adjust its value to the power available in the network. 


Requirements for safeguarding system 


It is necessary to control the following : 


+ 


2 
Ф 


pump operation parameters of cooling system 
temperature of stator's winding of electric motor 
temperature of lubricating oil or bearing 

temperature of transformer's winding 

overloads of frequency converter 

temperature of frequency converter's cooling medium. 
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Automatic control systems for ships fitted with podded propulsion drive (POD) 


If any ofthe above enumerated parameters is not contained 
within its permissible range the system has to start procedure 
of propulsion power lowering. 

The control system should be capable of auto-diagnosing. 
Due to the function it is possible to monitor : all programs 
executed by the system, values of voltage on auxiliary devices, 
supervision over communication system , I/O panel of frequ- 
ency converter, as well as analogue input signals. Besides it is 
important to ensure possiblity of recording course of the sys- 
tem's operation in order to perform later its analysis in the case 
if ship's crew is not able to cope with a failure. 


Starting and stopping main drive 


Starting the engine is usually executed from Central 
Control Station (CCS). АП auxiliary devices such as : 
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cooling water pump of engine cooling system 
cooling water pump of inverter cooling system 
fan of engine cooling system 

auxiliary voltage source of inverter 

pumps and fans of cooling system 


2% 


= 


2% 
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X 


% 


X 
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are to be started separately before starting the main engine. 


The following sequence of the main drive starting 
operation should be realized : 


X 


% 


switching-on an appropriate number of generating sets 
starting auxiliary devices 

releasing the engine's breaking system 

starting the main engine. 


X1 


% 


2%” 


% 


2% 
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Stopping the engine should consist іп lowering 
its rotational speed up to zero and next in the following : 


* switching-on the breaking system 
** switching-off auxiliary devices 
** successive switching-off the generating sets. 


Synchronization of drives 


In the case of applying a greater number of drives on a ship 
itis necessary to synchronize their work to avoid beating which 
happens if their rotational speeds are incompatible. To syn- 
chronize their work the following conditions are to be com- 
plied with : 


= synchronization function is to be activated by means of 
a graphical operator panel 

c» difference between the reference rotational speed and actu- 
al speeds of drives cannot be greater than the value establi- 
shed by their producer 

c» value ofthe reference rotational speed of each drive cannot 
differ by more than the number of rotations per minute (rpm) 
established by its producer 

c» actual values of rotational speed of drives cannot differ by 
more than the number of rpm established by their producer 

= communication between drives is to be realized through 
a data bus. 


However if synchronization of drives were necessary then 
it would be rather impossible to keep the ship on a set course 
by changing the rotational speed of one of the drives only. 


POSSIBLE CHOICE 
OF A CONTROL SYSTEM 


During the analysis of control and automation systems for 
a POD-fitted ship the following products of two leading pro- 
ducers were considered : 
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€ PodCon and Marintronics™ systems of ABB firm 
€ Alstom - Mermaid™ system of Rolls-Royce Group. 


All the products totally cover the automation system of ship 
propulsion and manoeuvring, and they are intended for POD 
drives. Additionally the products can be shaped into many con- 
figurations and matched up to specific features of a ship. More 
information on the systems can be found in [14] and the produ- 
cer technical pamphlets referred to in the below attached Bi- 
bliography. 

From the information contained in the sources it results 
that both firms : ABB and Rolls-Royce-Alstom intend to in- 
troduce the global supervising systems. It makes it possible to 
greatly improve efficiency of a power system, and to reduce 
number of ship's personnel in the future. Even a not large de- 
crease of shipping costs may improve competitiveness of a gi- 
ven shipowner. 

Itis very important to correctly manage the power in chan- 
geable load conditions (frequent speed changes, switching on/ 
off the thrusters etc), as then profits associated with lowering 
fuel oil consumption will be manifested the most. Correct ad- 
justing a number of generating sets which supply the network, 
depending on instantaneous power demand 15 the area of the 
greatest possiblities in improving the power system effective- 
ness. 

In such waters as the Baltic Sea, time of sea voyage of 
a ship 15 rather short, hence the share of duration time of port 
manoeuvers is substantial. Therefore application of an integra- 
ted ship management system may bring in measurable profits 
in the form of lowering operational costs to such an extent as 
to cover initial cost of the system. 


Choice of a definite system is not easy; 
however the problem goes beyond the scope of this paper. 
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FOUNDATION FOR SAFETY OF NAVIGATION 
AND ENVIRONMENT PROTECTION 
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CENTRE 


SHIP HANDLING RESEARCH 
AND TRAINING CENTRE AT ILAWA 


The safe handling of ships depends on many factors such as 
ship's manoeuvring characteristics, human factor (operator expe- 
rience and skill, his behaviour in stress situations etc.), as well as 
environmental conditions, and possible water area restrictions. The 
interface between ship operator and the available control devices 
is also a very important factor. 

Results of analysis of CRG (collisions, rammings and groun- 
dings) casualties show that human errors are involved in one third 
of them, and the same amount of CRG casualties is attributed to 
the poor controllability of ships. Training on ship handling is large- 
ly recommended by IMO as one ofthe most effective methods for 
improving the safety at sea. The goal of the above training is to 
gain theoretical and practical knowledge on ship handling in 
a wide number of different situations met in practice at sea. 

The main advantage of the method of training on ship han- 
dling based on manned models is realistic representation of all 
complex hydrodynamic phenomena affecting ship manoeuvrabi- 
lity. The models which satisfy similitude criteria have hydrody- 
namic characteristics very similar to those of real ships. This is 
very important e.g. for manoeuvres in restricted waters and ship- 
to-ship interaction, where bow cushion and suction areas conside- 
rably influence ship behaviour. Also manoeuvres where hull-pro- 
peller-rudder interaction plays a significant role can be reproduced 
on models in a realistic way. 

Training on manned models makes it possible to account for 
psychological aspects of training in comparison with that on an 
electronic simulator, by better feeling of effects of groundings, 
rammings and collisions, as well as environmental effects such as 
wind and current. 

Since 1980 more than 2000 ship masters and pilots from 27 
countries have been trained at Ilawa Ship Handling Research and 
Training Centre in weekly courses. The Foundation for Safety of 
Navigation and Environment Protection, being a non-profit organi- 
sation, is re-investing all spare funds in new facilities and each year 
new models and new training areas have been added to the existing 
ones. Also, the existing training models are modernised each year. 


Currently six models representing a wide spectrum of ship 
types are available at the Centre for training purposes. The ship 
models are presented in the table and figure. The models are 
equipped with all necessary devices simulating ship's systems, 
and the basic navigation aids (gyro, log, GPS, navigational lights, 
wind velocity and direction sensor, etc.). 

Last year a new system for simulation of tug-ship co-opera- 
tion entered into service. The system was installed on VLCC, LCC 
and Large Container models, and similar system will be soon in- 
stalled onboard other ship models. 

In summer 2000 the Training Centre at Паза was equipped 
with a high precision tracking system based on GPS technology. 
At present the GPS receivers working in RTK mode and allowing 
to measure the model's position with accuracy of 0.01m are in- 
stalled on two models, and in the next years similar receivers will 
be also fitted on other models. The new installation improved very 
much quality of training. 


The following training areas are currently available : 


€ Mock-ups of harbours with different configuration of entran- 
ces, turning areas, piers and quays, ferry terminals, docks and 
locks. The main exercises realised by using these training areas 
are : approaching harbours and terminals, berthing and unber- 
thing under various conditions, with wind and current pre- 
sent, entering the lock . 

€ Training areas including straight line and curvilinear water- 
ways and sharp turns. All passages are marked by buoys and 
leading marks for daylight exercises and a set of navigational 
lights for night-time exercises. Main current generator is also 
installed in this area. 

€ Open sea training area, which includes curvilinear deep water 
canal, sea berths, and a SPM (single point mooring buoy). 
This area is also used for realisation of standard manoeuvres, 
rescue and search manoeuvres with taking into account MER- 
SAR recommendations. Ship — to — ship manoeuvre can also 
be trained there. 


Main particulars of the models 


Model's 
type 


Length/Breadth Lgp/B 
Length/Draft Lgp/d 
Breadth/Draft B/d 

Block coefficient Св 
Number and type of 
propellers 
Dprop/ d 
Rudder area ratio 
Apg/(Lxd) [%] 
Scale 


5.79 


3.13 2.52 


1.76 1.39 


LCC tanker Киш Ro-ro ship Passenger VLCC tanker Large container 
bulk carrier car ferry 
|590 | 


|. — 1813 | 16 | 96 | 94 | 1573 | _ 2098 | 
| зз | 252 | 36 | 4% | 27 |) 29 | 


82 2x1.33 1.75 


1.85 
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WARTA 


Model of a LCC tanker 


SZCZECIN 


Model of a Panamax bulk carrier 
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Model of a VLCC tanker 


Model of a large container ship 


Ship training models used at Ilawa Centre 


€ Shallow water curvilinear canal. It consists of restricted cross- 
-section canal with a bend and straight sections, 140 metres in 
length. The canal is mainly used for passing and overtaking 
manoeuvres. 

€ River training area consisting of curvilinear stream restricted 
from both sides where current is created. Piers for mooring 
ship models in current are installed there. 

€ Restricted water areas for demonstration of bank and shallow 
water effects. 


On the shore site there are two hangars for storing models, as 
well as slipways and the building which contains lecture, recrea- 
tion and computer rooms, staff accommodations, stores, small 
workshop, and other facilities. 


The scope of lectures and programme of practical exercises is 
flexible. Each year new items are including into the course pro- 
gramme, so that the current programme differs much from the 
training realised some years ago. The new training areas together 
with the new models of modern ships (the seventh model repre- 
senting a LNG carrier is under construction now) and the mo- 
dernized and new equipment of training models make now possi- 
ble the realisation of training on specific types of ships. It is also 
possible to arrange special courses focused on chosen problems 
of manoeuvrability as well as to realize training on specific types 
of ships. 


Jacek Nowicki 
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